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Laboratory  (  formerly  Antenna  Laboratory)  ,  Department  of  Electrical 
Engineering,  The  Ohio  State  University  at  Columbus,  Ohio.  Research 
was  conducted  under  Project  5635  of  the  Air  Force  Cambridge  Research 
Laboratories,  Office  of  Aerospace  Research,  United  States  Air  Force, 
Bedford,  Massachusetts.  Dr.  John  K.  Schindler,  CRDG,  was  the 
Microwave  Physics  Laboratory  Program  Monitor  and  the  research  was 
conducted  under  Contract  Fl 9628-67-C -0239. 

It  is  a  pleasure  to  acknowledge  the  assistance  and  cooperation 
of  Professor  E.  M.  Kennaugh  in  the  preparation  of  this  report.  The 
author  is  also  indebted  to  the  staff  of  the  ElectroScience  Laboratory. 
Helpful  discussions  with  Professor  Jack  H.  Richmond,  Professor  Leon 
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The  material  contained  in  this  report  is  also  used  as  a  dissertation 
submitted  to  the  Department  of  Electrical  Engineering,  The  Ohio  State 
University  as  partial  fulfillment  for  the  degree  Doctor  of  Philosophy. 


ABSTRACT 


\ 

A  linear  system  analysis  is  applied  to  electromagnetic  scattering 
problems.  The  relationship  between  the  incident  and  scattered  fields  is 
reduced  to  a  one-dimensional  linear  system  by  suitable  restrictions,  and 
modeled  as  a  linear  two-port  with  time -invariant  parameters.  The  two- 
port  is  characterized  by  a  real,  time -dependent  impulse  response  function 
and  a  transform-related  complex,  frequency-dependent  phasor  response. 

Two  fundamental  approaches  for  obtaining  approximate  solutions  to 
electromagnetic  scattering  problems  are  seen.  Estimates  of  either  the 
frequency-dependent  phasor  response  or  the  time-dependent  impulse  re¬ 
sponse  can  be  attempted.  The  material  in  this  study  is  concerned  with 
the  latter  approach. 

\ 

The  concept  of  irhpulse  and  transient  response  approximations  in 
electromagnetic  scattering  problems  and  the  general  properties  of  such 
waveforms  are  reviewed  in  Chapter  I.  It  is  demonstrated  that  the  impulse 
and  transient  response  waveforms  provide  a  primary  conceptual  model  for 
electromagnetic  scattering  which  permits  the  integration  of  all  existing 
frequency  data,  calculated  or  measured,  and  is  simple  in  form.  The  model 
is  directly  related  to  the  geometrical  and  physical  properties  of  the  scatter¬ 
ing  object;  does  not  change  drastically  for  small  perturbations  of  parameters; 
and  can  be  used  to  predict  the  response  of  the  object  to  any  type  of  incident 
plane  wave.  Analytical  and  experimental  methods  for  obtaining  estimates 
of  the  response  waveforms  are  discussed,  and  a  systematic  procedure  for 
securing  such  estimates  is  suggested. 

Estimates  of  the  impulse  and  transient  response  waveforms  of  a 
perfectly  conducting  prolate  spheroid  are  obtained  in  Chapter  IL  Two 
results  are  derived.  The  first  is  restricted  to  axial  incidence  and  utilizes 
the  geometrical  similarity  of  the  sphere  and  spheroid  to  obtain  the  response 
waveforms  of  the  spheroid.  The  second  result  is  for  arbitrary  orientation 
of  the  spheroid  and  arbitrary  linear  polarization  of  the  incident  field. 
Calculations  from  both  results  are  compared  with  measured  data  on  a  2:1 
axial  ratio  spheroid. 

Derivation  of  the  response  waveforms  of  the  prolate  spheroid 
target  requires,  as  input  data,  the  Rayleigh  scattering  coefficient  of  the 
target.  The  Rayleigh  coefficients  of  prolate  and  oblate  spheroids  for 
arbitrary  orientation,  axial  ratio,  and  material  properties;  and  for 
arbitrary  linear  polarization  of  the  incident  field  are  given  in  Appendix  I. 

Time  domain  interpretations  of  high  frequency  asymptotic  estimates  of 
the  creeping  wave  on  conducting  spheres  and  cylinders  are  considered  in 
Appendix  II.  The  non-causal  nature  of  creeping  wave  estimates  for  the 
sphere  is  noted  and  suggestions  for  obtaining  a  causal  response  waveform 
are  made. 
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CHAPTER  I 
INTRODUCTION 


Theoretical  and  experimental  studies  of  the  interaction  of  electro¬ 
magnetic  waves  and  finite  material  objects  have  received  the  attention 
of  numerous  investigators  for  several  decades.  One  therefore  finds  in 
the  literature[  l]  an  imposing  array  of  sophisticated  mathematical 
techniques  for  calculating  such  interactions  and  an  equally  impressive 
collection  of  calculated  and  measured  data  on  the  interactions  for  spe¬ 
cific  objects.  These  data  are  usually  presented  in  the  form  of  curves 

or  patterns  of  the  radar  cross  section  of  the  object.  With  the  exception 

« 

of  spherical  objects,  however,  the  theoretical  contributions  are  charac¬ 
terized  by  restrictions  on  the  orientation  and  material  properties  of  the 
object  and  on  the  frequency,  polarization  and  character  of  the  incident 
electromagnetic  field.  In  the  same  vein,  experimental  studies  reported 
in  the  literature  cannot  hope  to  be  all  inclusive  and  inevitably  involve 
some  ot  the  same  limitations. 

*  Contributions  prior  to  1957  are  covered  by  earlier  Georgia  Institute 
of  Technology  bibliographies  also  given  in  Reference  1. 
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One  might  summarize  the  present  state  of  the  art  in  prediction 
of  the  radar  cross  section  of  an  object  as  follows;  if  the  object  is  very 
large  in  terms  of  the  wavelength,  then  several  approximate  or  asymp¬ 
totic  theories  of  varying  degree  of  complexity  and  accuracy  are  avail¬ 
able  for  cross  section  calculations.  The  major  difficulty  encountered 
in  this  portion  of  the  spectrum  is  in  estimating  how  large  a  particular 
object  must  be  in  terms  of  the  wavelength  in  order  to  achieve  a  given 
accuracy  in  the  cross  section  calculation  using  a  specific  approxi¬ 
mation.  If  the  object  is  very  small  in  terms  of  wavelength  then  again 
several  approximate  theories  are  available.  In  this  case  no  particular 
difficulty  is  encountered  in  the  accuracy  of  the  calculations  but  the 
maximum  object  size  which  can  be  handled  is  rather  strictly  limited. 
Such  limitations  are  imposed  by  either  the  geometrical  shape  and 
constitutive  parameters  of  the  object  or  by  the  storage  capacity  of 
present  digital  computers. 

Exceptions  exist  of  course,  but  in  general  one  can  achieve  a 
reasonable  estimate  of  the  radar  cross  section  of  an  object  at  the 
extremities  of  the  spectrum.  In  what  is  somewhat  loosely  referred 
to  as  the  resonance  region[  2]  of  an  object,  however,  no  such  reason¬ 
able  estimate  is  readily  obtainable.  In  certain  cases,  i.  e.  ,  for  a 
specific  object  geometry,  material  properties  and  orientation,  the 
extension  of  a  particular  high  or  low  frequency  approximate  or 


asymptotic  theory  into  the  resonance  region  has  been  achieved.  No 
general  extension  can  be  implied  from  these  results  however,  and  in 
fact  would  appear  to  fail  by  the  very  nature  of  the  assumed  fields. 

It  would  serve  no  purpose  here  to  detail  all  of  the  approximate 
theories  referred  to  above.  The  point  to  be  made  is  that  for  the 
general  finite,  non- spherical  object  a  significant  portion  of  the  spec¬ 
trum  usually  cannot  be  analyzed  by  existing  asymptotic  or  approxi¬ 
mate  scattering  theories.  Unfortunately,  this  gap  in  the  spectrum 
is  often  of  most  practical  interest  for  cross  section  calculations. 
Three  additional  comments  are  in  order.  First,  it  is  customary  with 
most  approximate  methods  to  assume  a  plane  monochromatic  incident 
field.  Treatment  of  an  arbitrary  type  signal  often  becomes  a  much 
more  complex  problem.  Second,  the  approximate  methods  in  vogue 
at  present  are  not  capable  of  systematically  integrating  calculated  and 
measured  data  from  various  portions  of  the  spectrum  into  a  single 
model  with  a  wider  applicability.  Thus,  for  example,  good  estimates 
of  high  and  low  frequency  scattering  properties  of  an  object  arc  not 
combined  to  improve  the  estimate  of  its  resonance  region  behavior. 
Finally,  the  geometrical  shape  of  an  object  is  not  simply  related  to 
the  form  of  the  predicted  scattered  field.  Consequently,  minor 
perturbations  of  the  object  geometry  often  require  a  tedious  new 

* 

solution  rather  than  a  minor  modification  of  an  existing  solution. 
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It  is  the  opinion  of  the  author  that  no  measurement  or  calculation 
program,  regardless  of  its  scope,  can  hope  to  satisfy  the  present  and 
future  requirements  in  radar  reflectivity  without  the  establishment  of 
a  primary  conceptual  model  whereby  all  the  existing  theories  and 
measurements  on  a  given  object  and  on  objects  of  similar  geometry 
can  be  integrated.  Such  a  conceptual  model  should  also  incorporate  an 
applicability  to  all  types  of  signaling  waveform  including  the  ultra-short 
pulse  of  current  interest. 

The  subject  of  this  study  is  the  interpretation  and  application  of 
such  a  primary  conceptual  model  for  electromagnetic  scattering.  This 
model  permits  the  integration  of  all  existing  frequency  data,  calculated 
or  measured,  and  is  simple  in  form.  It  is  directly  related  to  the  geo¬ 
metrical  shape  and  constitutive  parameters  of  the  object;  does  not 
change  drastically  for  small  perturbations  of  parameters;  and  can  be 
used  to  predict  the  response  of  the  object  to  any  type  of  incident  plane 
electromagnetic  wave.  This  primary  conceptual  model,  the  impulse 
and  related  transient  response  waveforms  of  a  scattering  object,  was 
first  proposed  by  Kennaugh  and  Cosgriff[  3]  in  1957.  In  their  paper, 
one  facet  of  the  impulse  response  concept,  namely  the  use  of  high  and 
low  frequency  approximations  to  secure  approximate  resonance  region 
data,  was  exploited.  Since  that  time,  Kennaugh  and  his  colleagues  at 
The  Ohio  State  University  ElectroScience  Laboratory  have  published 


5 


results  dealing  with  various  aspects  and  applications  of  the  theory[4-l6] . 
A  review  of  these  accomplishments;  a  description  of  the  general  proper¬ 
ties  of  the  solutions  to  transient  electromagnetic  scattering  problems; 
and  a  discussion  and  illustration  of  various  methods  for  approximating 
the  appropriate  waveforms  was  published  by  Kennaugh  and  Moffatt[  13] 
in  19^5.  To  the  best  of  the  author’s  knowledge,  only  three  other  publi¬ 
cations  dealing  specifically  with  the  impulse  response  waveform  of  a 
finite  three  dimensional  scattering  object  have  appeared.  Betten[  17] 
in  1962  discussed  the  experimental  determination  of  the  impulse  re¬ 
sponse  waveform  of  an  object  and  the  synthesis  of  an  equivalent  net¬ 
work.  He  also  suggested  a  rather  elaborate  identification  scheme 
utilizing  cross  correlation  of  object  waveforms.  Brown[  18]  in  1962 
presented  an  approach  to  the  solution  of  pulse  scattering  by  finite 
obstacles  which  utilized  the  identification  and  separate  treatment  of 
the  individual  terms  in  a  wavefront  expansion  of  the  transforms  of  the 
field  vectors.  Thus  latter  work,  while  mathematically  quite  elegant, 
adds  little  if  any  insight  to  the  estimation  of  impulse  response  wave¬ 
forms  for  finite  objects.  Finally,  in  1963  Barabanenkov,  et  al.  [  19] 
derived,  apparently  independently,  the  impulse  response  waveform 
produced  by  the  physical  optics  approximation  as  obtained  by  Kennaugh 
and  Cosgriff  in  1957. 
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In  this  study,  the  first  chapter  develops  the  necessary  groundwork 
for  investigating  solutions  to  transient  electromagnetic  scattering  prob¬ 
lems.  The  general  properties  of  such  solutions  are  reviewed  and  defi¬ 
nitions  of  the  pertinent  transient  waveforms  are  given.  Analytical  and 
experimental  techniques  for  obtaining  approximate  estimates  of  the 
impulse  and  transient  response  waveforms  are  discussed. 

In  Chapter  II,  the  electromagnetic  backscattering  by  a  perfectly 
conducting  prolate  spheroid  is  analyzed,  using  estimates  of  the  impulse 
response  and  related  transient  response  waveforms.  Two  results  are 
derived.  The  first  is  restricted  to  axial  incidence  and  utilizes  the  geo¬ 
metrical  similarity  of  the  sphere  and  spheroid  to  obtain  an  impulse 
response  for  the  spheroid  in  terms  of  known  approximate  waveforms 
for  the  sphere.  The  second  result,  based  in  part  on  the  results  of  the 
axial  solution,  is  for  arbitrary  axial  ratio  and  orientation  of  the 
spheroid  and  for  arbitrary  linear  polarization  of  the  incident  electro¬ 
magnetic  field.  Calculations  from  both  results  are  compared  with 
measured  data  for  a  2:1  axial  ratio  spheroid  and  previously  calculated 
data  for  a  10:1  axial  ratio  spheroid. 

Two  Appendices  are  included  in  this  study.  In  the  first  Appendix, 
the  Rayleigh  scattering  coefficients  of  prolate  and  oblate  spheroids  are 
obtained.  The  results  hold  for  arbitrary  orientation,  axial  ratio,  and 


material  properties  of  the  spheroid  and  for  arbitrary  linear  polari¬ 
zation  of  the  incident  wave.  These  results  are  utilized  in  Chapter  II 
in  deriving  response  waveforms  of  the  prolate  spheroid. 
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The  second  Appendix  contains  a  time  domain  intepretation  of 
certain  high  frequency  asymptotic  estimates  of  the  creeping  wave  on 
cylinders  and  spheres.  The  non-causal  nature  of  particular  asymp¬ 
totic  estimates  is  demonstrated  and  suggestions  for  obtaining  a 
causal  response  waveform  are  made. 

A.  The  Concept  of  Impulse  and 

Transient  Response  Waveforms 
in  Electromagnetic  Scattering 
Problems 

Two  basic  assumptions  are  made  which  hold  throughout  this 
investigation.  It  is  assumed  that  between  the  input  terminals  of  the 
radiating  system  producing  an  incident  field  and  the  output  terminals 
of  the  receiving  system  detecting  the  scattered  field  only  linear 
processes  with  time -invariant  parameters  are  involved.  The  scat¬ 
tering  process  is  viewed  as  a  cause  and  effect  relationship  whereby 
a  transducer  operates  on  the  input  or  cause  at  the  input  terminals  to 
produce  the  output  or  effect  at  the  output  terminals.  With  the  as¬ 
sumptions  given  above,  the  transducer  is  a  linear,  time -invariant 
operator  and  is  uniquely  determined  from  the  knowledge  of  a  single 
function;  the  impulse  response.  Specifically,  the  cause  is  the  incident 
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electromagnetic  field  in  which  the  scatterer  is  immersed  and  the  effect 
is  the  scattered  electromagnetic  field  defined  conventionally  i  the 
difference  between  the  total  and  incident  fields.  It  is  further  assumed 
that  the  input  and  output  terminals  are  both  in  the  far-field  of  the  scat¬ 
tering  object. 

Thus,  as  shown  in  Fig.  1,  the  scattering  process  is  modeled  by 
a  passive  linear  two-port  with  time -invariant  parameters.  The  input 
is  E(  t')  ,  the  output  F(  t')  ,  and  the  two-port  has  an  impulse  response 
function  Fj(  t ' )  .  The  input  and  output  are  related  through  the  corvo  - 
lution  integral  as 

00 

(1)  F(t')  Fj(t)  E(f-T)dT 

ti 


F(t') 


Fig.  l--Two  port  representation  of  scattering  problem. 


In  order  to  define  the  various  quantities  in  Eq.  (  1)  and  to  explain  the 
scalar  treatment  of  the  electromagnetic  fields,  it  is  necessary  to 
introduce  a  coordinate  reference  frame.  This  coordinate  frame  is 
shown  in  Fig.  2.  The  incident  plane  electromagnetic  field  polarized 
in  the  x  direction 

(2)  El  =  x  E(t  -  z/c) 


Fig.  2--Coordinates  of  scattering  problem. 
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evaluated  at  z  =  0  has  an  intensity  E(  t) .  With  this  plane  wave  incident 
on  the  scattering  object  as  shown  in  Fig.  2,  a  scattered  field  is  pro¬ 
duced.  At  an  arbitrary  location  P  in  the  far-field  of  the  scattering 
object,  a  normalized  transverse  component  of  the  scattered  field  has 
an  intensity  F(  t) .  In  order  to  remove  the  time  delay  between  scatterer 
and  observer,  a  new  time  scale  t1  =  t-r/c  is  introduced.  The  input 
E(  t')  is  simply  E(  t)  with  t  replaced  by  t'.  For  the  output  F(  t')  , 
t-r/c  is  replaced  by  t'  in  F(  t) .  The  impulse  response  waveform  Fj(t') 
is  the  response  when  the  input  E(t')  is  impulsive,  i.  e.  ,  E(t')  =  6(t'). 
The  lower  integration  limit  t\  in  Eq.  (1)  is  the  initial  value  of  t'  at 
which  the  impulse  response  waveform  Fj(t')  departs  from  zero.  This 
limit  is,  in  general,  not  zero  since  the  initial  contribution  need  not 
arrive  at  a  time  t  =  r/c.  If  L  is  the  maximum  object  dimension  then 
tj  >  -L/c  in  the  half-space  z  >  0  and  t[  >  0  in  the  half  space  z  <  0. 

The  two-port  has  a  frequency-dependent  phasor  response  G(  jo>)  which 
is  related  to  the  radar  cross-section  of  the  scatterer  as 

(3)  <r(  w)  -  it  |G(  ju  )  | 

The  frequency-dependent  phasor  response  G(  jw)  and  the  time- 
dependent  impulse  response  waveform  Fj(t')  form  a  Laplace  trans¬ 


form  pair. 
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(4)  G(ju)  =  c  ^  Fjtt')  e'->wt'  df 

0 

CO 

F,(t')=  _L_  V  G(j“)  etlj“  du 

2"jc 

or  using  s  notation  (  s  =  jw ) 

CO 

(  5)  G(  s)  =  cj  Fj(t')  e"st'  df 

0 

1 

Ft(  t')  =  -  \  G(  s)  eVs  ds 

1  2trjc  J 

-j°o 

Equations  (  3)  and  (  4)  define  the  normalization  mentioned  above.  It 
should  be  noted  that  the  scattering  object  in  Fig.  2  is  located  in  the 
half- space  z  >  0  of  the  coordinate  frame.  A  more  detailed  develop¬ 
ment  of  the  normalized  impulse  response  waveform  Fj(t')  is  found 
in  Kennaugh  and  Moffatt[  13]  . 

Certain  points  should  be  stressed  concerning  the  two-port  in 
Fig.  1.  Its  time-dependent  impulse  response  Fj(  t')  is  dependent 
on  the  orientation  of  the  scattering  object,  the  observation  angles 
(  but  not  range)  and  the  particular  transverse  component  of  the 
scattered  field  selected.  The  electromagnetic  scatterer  properties 
are  uniquely  determined  by  Fj(  t')  ,  but  one  cannot  extrapolate  this 
statement  without  proof  to  imply  that  Fj(  t')  uniquely  determines  the 
physical  properties  of  the  scatterer.  It  seems  doubtful  that  two 


12 


geometrically  or  materially  different  scatterers  could  (  for  fixed  source 
and  receiver  locations  and  polarizations)  ,  have  exactly  the  same  phasor 
response  over  the  entire  spectrum,  but  a  unique  determination  of  the 
physical  properties  of  the  scatterer  has  not  been  proven. 

The  impulse  response  waveform  Fj(t')  is  the  ime-dependent 
electromagnetic  field  intensity  produced  at  the  output  terminals  when 
the  input  E(t')  is  an  impulsive  plane  electromagnetic  wave,  i.  e.  , 

E(t')  =6(t').  OnceFi(t')  is  known,  the  response  waveform  for  any 
incident  waveform  is  determined  by  Eq.  (  1).  Two  other  particular 
response  waveforms  are  of  interest  in  this  study  and  need  to  be  defined. 
The  step  response 

00 

(6)  Fu(t')  =  ^  Fx(  t)  u(  t'-T)  dT 

t; 

and  the  ramp  response 

®  t' 

(7)  F r(  t' )  =y  Fj(  t)  ( t'-T)  u(  t'-T)  dT  =  ^  Fy(  T)dT 

ti  t; 

The  Laplace  transform  relations  in  Eqs.  (  4)  and  (  5)  state  that 
Fj(t')  andG(jw)  can  each  be  derived  from  the  other.  ButG(jw)  is 
known  exactly  for  only  one  finite  three  dimensional  shape;  the  sphere, 
and  even  for  this  shape  the  transformation  to  obtain  Fj(  t')  cannot  be 
achieved  exactly.  Thus,  a  study  of  the  scattering  problem  in  the  time 
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domain  consists  essentially  of  the  development  of  a  reasonable  esti¬ 
mate  for  the  impulse  response  waveform  Fj(t').  In  exactly  the  same 
way,  the  study  of  scattering  problems  in  the  frequency  domain  has 
been,  as  discussed  in  the  introduction  to  this  chapter,  the  development 
of  reasonable  estimates  to  the  phasor  response  G(jw),  but  usually  for 
either  the  high  or  low  extremities  of  the  spectrum.  A  number  of 
distinct  advantages  to  the  time  domain  approach,  as  well  as  some  dis¬ 
advantages,  will  be  enumerated  in  the  course  of  this  study.  It  is  hoped 
that  the  time  domain  approach  can  be  presented  with  sufficient  clarity 
and  simplicity,  and  its  utility  demonstrated  on  a  sufficiently  complex 
problem,  that  future  investigators  will  move  freely  between  the 
domains  to  exploit  the  maximum  potential  of  either  approach. 

It  would  be  premature  at  this  point  to  discuss  the  relative  ad¬ 
vantages  and  disadvantages  of  a  frequency  domain  versus  time  domain 
approach  to  an  electromagnetic  scattering  problem.  What  is  hoped  to 
be  demonstrated  is  that  one  can  incorporate  into  the  time  domain  model 
most  of  the  advantages  of  a  frequency  domain  approach  while  at  the 
same  time  exploiting  certain  unique  features  of  the  time-dependent 
waveform  model.  Consideration  at  the  moment  is  to  be  restricted  to 
simply  the  solution  of  a  given  electromagnetic  scattering  problem 
involving  a  single  isolated  object.  In  the  opinion  of  the  author,  the 


following  remarks  justify  consideration  of  the  time  domain  approach 
and  at  the  same  time  indicate  the  basic  premises  upon  which  such  an 
approach  is  based. 

1.  It  is  to  be  generally  expected  that,  compared  to  the 
complicated  complex  phasor  response  G(  s)  of  a  given 
object  as  a  function  of  frequency,  the  real,  time-de¬ 
pendent  waveform  Fj(t')  will  be  relatively  simple.  In 
the  context  of  this  study,  the  scatte^er  can  be  represented 
by  a  distributed  constant,  linear,  time -invariant  two-port 
network.  As  such,  its  impulse  response  waveform  should 
be  relatively  simple  compared  to  the  transcendental 
response  functions  characterizing  distributed  systems  in 
the  frequency  domain. 

2.  The  impulse  and  transient  response  waveforms  of  a 
scattering  object  must  be  related  in  a  rather  direct  way 
to  the  geometry  and  constitutive  parameters  of  the  object. 

It  will  be  shown,  for  example,  that  the  area  beneath  the 
ramp  response  waveform  is  proportional  to  the  Rayleigh 
coefficient;  and  hence,  to  the  volume  of  the  scatterer. 

More  fundamentally,  however,  as  the  impulsive  or 
transient  illumination  moves  across  the  object,  only  that 
portion  of  the  object  which  has  been  illuminated  can  possibly 
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contribute  to  the  scattered  field  waveform.  Therefore*  up 
until  the  time  when  the  transient  illumination  has  passed 
completely  over  the  object,  there  is  a  direct  correlation 
between  the  response  waveform  at  a  given  time  and  a  spe¬ 
cific  portion  of  the  object.  Furthermore,  since  we  are 
dealing  with  a  physically  realizable  system  which  cannot  be 
a  predictor,  two  objects  which  present  initially  identical  geo¬ 
metrical  and  physical  properties  over  a  given  region  must 
yield  identical  response  waveforms  up  to  the  time  corre¬ 
sponding  to  complete  illumination  of  this  region,  regardless 
of  their  geometrical  and  physical  properties  beyond  this  point. 
3.  In  principle  at  least,  it  is  possible  to  incorporate  into  an 
estimate  of  the  waveform  Fj(  t')  all  of  the  best  features  of 
various  approximate  or  asymptotic  estimates  of  G(  s)  while  at 
the  same  time  utilizing  ctrtain  unique  features  of  the  time -de¬ 
pendent  waveform.  If  one  has  certain  estimates  of  G(  s)  whose 
validity  is  restricted  to  particular  portions  of  the  spectrum,  it 
is  far  from  clear  how  consideration  of  these  estimates  can  be 
used  to  approximate  a  G(  s)  corresponding  to  the  remainder 
of  the  spectrum.  In  the  time  domain  however,  the  estimates 
of  G(  s)  become  time -limited  port  .-ns  of  the  waveform,  and  it 
is  known  that  these  pieces  must  combine  with  other  pieces  to 
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produce  a  single  waveform.  Certain  conditions  on  this  total 
waveform  are  known  from  low-frequency  derived  moment 
conditions.  Even  very  crude  estimates  of  how  the  pieces 
are  combined  must  add  some  new  knowledge  concerning  G(  s) . 

4.  The  convolution  integral  in  Eq.  (  1)  relating  the  response 
waveform  to  the  interrogating  signal  provides  an  additional 
understanding  of  the  relationship  between  the  contributions 
of  G(  s)  from  various  portions  of  the  spectrum  and  the  response 
waveform  Fj(  t').  Let  the  input  signal  be  a  monochromatic 
continuous  wave.  The  graphical  interpretation  of  convolution 
is  that  of  reversing  one  of  the  signals  with  respect  to  time  and 
then  sliding  one  over  the  other.  At  any  given  time,  the  response 
is  given  by  the  integral  of  the  product  of  the  two  waveforms  over 
that  time  interval  where  they  coincide.  For  the  monorhromatic 
input,  this  consists  of  sliding  a  sinusoid  of  a  given  period  across 
the  reversed  waveform  Fj(  -t').  It  follows  that  at  relatively  low 
frequencies,  the  response  can  be  influenced  little  by  the  minute 
details  of  the  waveform;  therefore  the  response  at  low  fre¬ 
quencies  is  basically  dependent  on  the  general  size  and  shape 
of  the  waveform.  As  the  input  frequency  increases,  more  and 
more  of  the  waveform  detail  is  important  and  slowly  varying 
portions  of  the  waveform  less  important  since  the  contributions 
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from  these  are  effectively  cancelled  by  the  positive  and  negative 
portions  of  the  sinusoid.  Two  conclusions  can  be  drawn  from 
these  remarks[  13]  :  (1)  The  general  shape  and  size  of  the 
impulse  response  waveform  Fj(  t1)  is  dictated  by  the  low-fre¬ 
quency  response  of  the  object  and  (  2)  The  fine  structure  and 
detail  of  the  waveform  is  controlled  by  the  high-frequency 
response  of  the  object.  Thus  if  one  constructs  an  estimate 
of  Fj(t')  using  the  moment  conditions,  it  is  to  be  expected 
that  such  a  model  would  fit  relatively  well  at  low  frequencies, 
but  would  require  modifications  and  corrections  to  predict 
the  response  well  at  high  frequencies. 

5.  The  use  of  experimental  data  in  deriving  estimates  of  G(  s) 
or  Fj(t')  is  an  important  consideration.  The  most  straight¬ 
forward  but  probably  least  effective  method  for  utilizing  such 
data  is  simply  the  comparison  of  theoretical  results  based  on 
a  model  obtained  wholly  from  analytical  considerations  with 
such  data.  If  only  a  limited  amount  of  experimental  results 
are  available,  such  comparisons  are  not  conclusive  and  often 
not  indicative  of  the  changes  required  in  the  model.  From 
the  discussion  in  (  4)  above,  it  is  apparent  that  one  should  be 
able  to  construct  a  crude  model  of  the  response  waveforms 
from  a  limited  number  of  low-frequency  measurements,  and 
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that  systematic  procedures  for  incorporating  such  data  with 
theoretical  estimates  of  Fj<  t')  should  be  feasible. 

This  study  does  not  fully  develop  the  methods  and  techniques 
mentioned  in  (  5)  above,  and  it  will  be  found  that  the  model  of  Fj(t') 
developed  in  Chapter  II  for  the  prolate  spheroid  target  is  obtained 
without  consideration  of  experimental  data;  such  data  only  being  used 
to  compare  with  calculations  from  the  model.  However,  we  shall 
discuss  how  such  experimental  data  might  be  more  effectively  used 
and  the  possible  directions  of  future  research  in  this  area. 

B.  General  Properties  of  the 
Impulse  and  Transient 
Response  Waveforms 

The  frequency-dependent  phasor  response  G(  jw)  and  the  corre¬ 
sponding  time -dependent  impulse  response  waveform  F j(  t 1 )  are  both 
governed  by  Maxwell's  equations.  Therefore  Fj(  t*)  must  be  a  causal 
function,  i.  e.  ,  a  backscattering  response  cannot  appear  at  the  output 
terminals  before  the  time  t'  =  0.  The  square  integrable  criterion 

CO 

(8)  -  j  |  G(  jw)  |  dw  <  ®  , 

0 

can  be  applied  provided,  as  pointed  out  by  Paopulis[  20 ]  ,  all  constants 
and  positive  powers  of  jw  are  first  removed  from  G(  jw ) .  Note  that  it 
does  not  necessarily  follow  that  approximate  or  asymptotic  estimates 
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of  G{  jw )  will  also  lead  to  causal  time  functions.  It  will  be  shown  in 
Appendix  II  that  certain  high-frequency  asymptotic  estimates  of  G(  jw) 
lead  to  non-causal  time  functions. 

The  definition  in  Eq.  (  3)  shows  that  G(  jw)  must  have  the  di¬ 
mension  of  length.  From  Eqs.  (  4)  ,  (  6)  ,  and  (  7)  respectively,  it  is 
seen  that  the  impulse  response  Fj(t')  is  dimensionless;  the  step  re¬ 
sponse  Fjj(t')  has  the  dimension  time;  and  the  ramp  response  Fj^(  t1) 
the  dimension  time  squared.  It  is  most  convenient  to  choose  a  time 
scale  measured  in  units  of  the  transit  time  for  some  characteristic 
dimension  of  the  scatterer.  If  I  is  such  a  characteristic  length  then 
the  transit  time 

(9)  t0  =  i/c. 

is  set  equal  to  unity.  This  can  be  interpreted  as  choosing  a  time 
scale  such  that  c  is  unity  and  a  distance  scale  such  that  a  free  space 
wave  traverses  the  length  I  in  unit  time. 

The  low  frequency  scattering  properties  of  any  finite,  three  di¬ 
mensional  object  provide  interesting  and  useful  conditions  on  the 
impulse  and  transient  response  waveforms.  At  sufficiently  low  fre¬ 
quencies,  it  is  known  that  the  phasor  response,  G(  s) ,  of  a  scattering 
object  can  be  expanded  in  a  Taylor  series  about  the  origin  s  =  0  as 
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CO 

(  10)  G(s)  ^2  an  *n 

n=0 

where  Rayleigh's  theory  dictates  that  the  coefficients  a0  and  aj  in 
such  an  expansion  are  zero  while  the  coefficient  a2  is  proportional 
to  the  scatterer  volume.  Note  that  the  coefficient  a2  depends  upon  the 
shape,  orientation,  and  constitutive  parameters  of  the  scatterer  as 
well  as  the  polarization  of  the  incident  and  scattered  fields.  Expand¬ 
ing  e"st'  in  the  definition  integral  in  Eq.  (5)  in  a  Taylor  series  about 
the  origin  s  =  0  there  results, 

CO 

(11)  y  F^t'Jdt'  =  a0  =  0 

0 


CO 
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These  are  moment  conditions  on  the  impulse  response  waveform 
Fj(t').  An  integration  by  parts  of  the  second  and  third  expressions 
in  Eq.  (11)  shows  that  the  first  three  moment  conditions  can  be 
interpreted  as  requiring  that  the  net  area  beneath  the  impulse  and 
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step  response  waveforms  be  zero  and  that  beneath  the  ramp  response 
waveform  proportional  to  the  so-called  Rayleigh  coefficient.  a2,  of 
the  scatterer.  A  word  is  in  order  concerning  this  proportionality 
factor  which  has  caused  some  difficulty  in  the  past.  With  the  definition 


in  Eq.  (  11)  ,  the  Rayleigh  term  (for  small  w)  is 


-ai(2ir)Vi*  -as(2ir)2  /i  \* 


(  12) 


G(ju)  = 


i2  X2 


where  t  is  some  characteristic  dimension  of  the  scattering  object.  If 
the  Rayleigh  coefficient,  K,  is  defined  by 


(  13) 

G(j»)  -Kfcj 

then 

-a2(  2 it)  2 

(14) 

K=  t2 

lo 

and 

00 

r* 

(15) 

^  FR(t')dt'  = 

0 

-K  to 
(  2it)  2  c 


Note  that  the  moment  conditions  in  Eq.  (11)  were  not  obtained  from 
the  final  value  theorem  of  the  Laplace  transform  theory 


lim  sG(  s)  =  lim  Fj(  t') 
s—0  t'—  °° 


(16) 
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which  theorem  simply  yields  the  result  that  in  the  limit  of  large  time 
the  response  must  decay  to  zero.  This  latter  result  is  clearly  reason¬ 
able  on  physical  grounds. 

The  expansion  for  G{  s)  in  Eq.  ( 10)  is  not  unique  unless  a  phase 
center  is  defined.  Since  the  first  non-zero  term  is  the  coefficient  of 
s2,  the  choice  of  phase  center  does  not  change  the  first  three  moments. 
However,  to  employ  moments  higher  than  the  third  both  the  choice  of 
phase  center  and  the  origin  of  coordinates  on  the  t'  scale  must  be 
specified.  The  first  two  moment  conditions  are  applicable  to  any 
finite,  three  dimensional  object.  For  a  great  many  targets,  the 
Rayleigh  term  (  coefficient  of  s2)  can  also  be  estimated  making  the 
third  moment  applicable.  The  higher  order  coefficients  can  only  be 
obtained  for  particular  objects,  therefore  the  usefulness  of  the  higher 
order  moment  conditions  (n  >  2)  is  restricted  to  a  small  class  of 
targets.  An  exception  to  this  exists  for  the  case  of  objects  which  have 
a  center  of  symmetry.  ‘  Stevenson[  21]  has  demonstrated  that  for 
this  class  of  targets  the  coefficient  of  the  s3  term  is  identically  zero 
when  the  phase  center  for  the  expansion  is  chosen  at  the  center  of 


*  An  object  has  a  center  of  symmetry,  0,  if  it  is  symmetrical  to  a 
transformation  consisting  of  rotation  by  tt  radians  about  an  axis 
through  0  and  then  reflection  in  a  plane  perpendicular  to  the  axis  of 
rotation  and  containing  0. 
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symmetry.  Therefore  for  this  class  of  objects,  the  fourth  moment 
condition  can  be  utilized.  If  the  phase  center  is  chosen  at  the  center 
of  symmetry  and  the  origin  of  the  t1  axis  at  the  point  where  the  inci¬ 
dent  impulsive  wave  first  encounters  the  body  then 

00 

(*  12  La, 

(  17)  j  f»  Frffjdf  =  c "  ■ 1 
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where  L  is  twice  the  linear  distance  from  the  center  of  symmetry  to 
the  t'  origin  and  is  a  negative  quantity.  Equation  (  17)  may  also  be 
written  as 


(  18) 


j  t'FR(f)dt'  = 


The  moment  conditions  given  in  Eq.  (  11)  implicitly  assume 
that  all  the  moments  of  the  time -dependent  impulse  response  wave¬ 
form  exist[  22]  .  It  is  clear  that  for  the  n^1  moment  to  exist,  Fj(t') 
for  large  t'  must  vanish  at  least  as  fast  as  l/t,n+^.  In  the  latter 
portion  of  this  section  a  discussion  of  the  proper  analytical  models 
for  Fj(t')  and  G(  s)  makes  it  clear  that  insofar  as  the  approximations 
of  interest  for  these  functions  are  concerned,  the  assumption  of  the 


existence  of  the  moments  is  not  unreasonable. 
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A  well  known  high  frequency  estimate  of  G(  s)  is  the  physical 
optics  approximation.  It  can  be  shown£  23]  that  this  approximation 
predicts  the  following  form  for  the  step  response  waveform  Fy(  t') 
for  conducting  objects  if  the  observation  location  is  constrained  to 
the  E  or  H  planes  of  the  incident  field. 

E  plane,  direct  polarized  component 

i  (•  c  f  p  •  ( g  - 1)  l  A  A 

(19)  Full',  B.P)  =  —  j  j  6  t' - - -  i-  n  ds 

8ill  L  c  J 

E  plane,  cross  polarized  component 


Fig.  3 --Coordinates  of  scattering  problem. 

point  on  the  scatter er  surface.  P  is  the  angle  between  the  unit  vectors 
g  and  £  (  P  <  ir)  and  slu  is  the  illuminated  surface  of  the  scatterer. 

For  the  backscattering  case  (  P  =  v)  specializing  to  the  orientation  in 
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The  function  A(  z)  measures  the  area  of  the  scatterer  surface  between 
the  Xf  y  plane  and  a  cutting  plane  at  z  projected  orthogonally  on  the  x,  y 
plane.  A  detailed  derivation  of  Eq.  (  24)  is  given  in  Reference  13. 

The  bistatic  results  in  Eqs.  (19)*  (21),  and  (  22)  have  also  been  in¬ 
terpreted  in  terms  of  area  functions[  23]  .  It  can  be  shown[  13]  that 
the  physical  optics  result  in  Eq.  {  24)  always  satisfies  the  zero 
moment  condition  [  first  of  Eq.  (  H)  J  but  does  not  in  general  satisfy 
the  higher  order  moments. 

Similar  time  domain  interpretations  can  be  made  of  other  high- 
frequency  approximations,  usually  however  only  with  reference  to  a 
particular  scattering  object.  If  the  inverse  transformation  [  Eq.  (  4)] 
can  be  achieved  in  closed  form  then  an  analytical  result  is  obtained. 

If  the  inverse  transformation  cannot  be  achieved,  a  Fourier  synthesis 
procedure  can  be  used  to  obtain  an  estimate  of  the  time-dependent 
waveform.  The  Fourier  synthesis  procedure  is  adequately  described 
in  Reference  13.  In  essence,  given  the  phasor  response  G(  s)  at  a 
sufficient  number  of  harmonically  related  frequencies,  estimates  of 
the  time -dependent  response  waveforms  can  be  generated. 

It  will  be  instructive  at  this  point  to  examine  certain  of  the 
approximate  impulse  and  transient  response  waveforms  which  have 
been  obtained  for  specific,  objects.  The  interest  here  is  in  noting  the 


salient  features  and  characteristics  of  the  waveforms.  It  is  stressed 
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that  these  are  far  or  radiation  field  response  waveforms  and  the  con¬ 
clusions  drawn  from  them  are  restricted  to  that  region.  The  back- 
scattering  impulfo,  step,  and  ramp  response  waveforms  of  a  perfect¬ 
ly  conducting  sphere[  13]  are  shown  in  Fig.  4.  The  sphere  radius  is 
a  and  the  time  scale  is  in  units  of  the  sphere  diameter  transit  time 
( tQ  =  2a/c) .  In  this  figure,  as  well  as  those  presented  in  the  re¬ 
mainder  of  this  study,  the  delta -function  pulse  of  weight  ot  is  denoted 
by  a  vertical  arrow  with  a  by  the  arrow  head.  For  the  impulse  re¬ 
sponse  waveform,  the  incident  waveform  is  a  positive  delta  function 
pulse  of  weight  unity.  The  impulse  response  waveform  consists  of  a 
negative  impulse  of  weight  a/c  followed  by  a  positive  jump  of  0.  5  and 

O 

a  secondary  negative  maximum  at  approximately  t'  equal  2.  57t0. 
Beyond  this  secondary  negative  peak  the  waveform  damps  rapidly, 
becoming  negligible  in  comparison  to  the  earlier  portions  for  times 
greater  than  approximately  6  transit  times  for  the  sphere  diameter. 
The  secondary  negative  peak  at  t1  ~  2.  57t0  corresponds  to  the  arrival 
of  a  diffracted  front.  Both  the  negative  impulse  and  0.  5  jump  at 
t '  =  0  are  correctly  predicted  by  approximate  high  frequency  solutions. 
The  response  waveforms  of  the  conducting  sphere,  particularly  the 
diffracted  or  creeping  wave  peak  at  t'  =  2.  57t0,  will  be  discussed  in 
Appendix  II.  One  should  note  that  the  exact  nature  of  the  creeping 
wave  maximum  is  not  resolved  in  Fig.  4.  It  is  important  to  discern 
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two  main  points  from  the  waveforms  shown  in  Fig.  4.  The  first  is  that 
the  successive  integrations  to  obtain  the  step  and  ramp  response  wave¬ 
forms  from  the  impulse  response  result  in  a  considerable  smoothing 
of  the  latter  waveforms.  Thus  distinct  singularities  and  peaks  present 
in  the  impulse  response  waveform  are  effectively  masked  in  the  corre¬ 
sponding  ramp  response  and  one  can  conclude  that  a  reasonable 
approximation  to  the  ramp  waveform  is  more  simply  obtained  that  an 
impulse  response  approximation.  Clearly,  each  waveform  contains 
all  the  information;  for  example  the  secondary  peak  in  the  sphere 
impulse  response  becomes  at  most  a  subtle  change  in  slope  in  the 
ramp  waveform,  but  an  estimate  of  the  ramp  waveform  which  ignored 
the  change  in  slope  would  still  be  a  reasonable  approximation  over  a 
range  of  lower  frequencies.  The  second  point  of  interest  is  the  ef¬ 
fective  time-limited  nature  of  the  response  waveforms.  For  times 
greater  than  5  or  6  to,  the  response  waveforms  have  decayed  to  a 
negligible  value.  Regardless  of  the  exact  nature  of  this  decay,  it  is 
evident  that  an  approximation  to  the  waveform  which  fits  well  for 
t*  <  6  tQ  and  which  vanishes  for  large  times  would,  for  all  practical 
purposes,  provide  a  satisfactory  model. 
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The  response  waveforms  shown  in  Fig.  4  were  obtained  by  a 
Fourier  synthesis  procedure  using  calculated  values  of  G(  s)  from  the 
Mie  series.  In  Fig.  5,  the  approximate  impulse  response  waveforms 
of  several  other  objects  are  shown.  The  waveforms  in  Fig.  5  were 
obtained  from  a  combination  of  high-frequency  estimates  of  G(  s)  and 
the  lower  order  moment  conditions.  In  the  case  of  the  cone-sphere[  8] 
(Fig.  5b)  and  the  prolate  spheroid[  1 0 ]  (Fig.  5d)  ,  the  sphere  results 
in  Fig.  4  were  also  utilized.  The  reader's  attention  is  called  to  the 
correction  in  the  response  waveform  shown  for  the  flat-based  right 
circular  cone  in  Fig.  5a.  The  response  waveform  for  this  object 
given  in  Reference  13  inadvertently  reversed  the  time  behavior  of  the 
doubly-diffracted  term  from  the  base. 

Some  additional  characteristics  of  the  impulse  response  wave¬ 
form  can  be  deduced  from  the  waveforms  shown  in  Fig.  5.  There 
is  clearly  a  general  relationship  between  the  object  geometry  and  the 
response  waveform;  i.  e.  ,  the  various  singularities  and  peaks  in  the 
response  waveform  occurring  at  given  times  can  be  directly  associated 
with  a  corresponding  geometrical  feature  of  the  object.  Furthermore, 
the  estimated  character  of  the  singularity  or  peak  can  also  be  deduced 
from  the  geometrical  feature.  For  example,  the  response  waveform 
associated  with  the  diffraction  from  smooth  convex  surfaces  ( the 
secondary  peak  in  the  sphere  impulse  response  of  Fig.  4)  has  the 
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same  general  form  for  the  sphere,  spheroid  and  cone- sphere.  In  a 
similar  fashion,  it  is  seen  from  Fig.  5  that  single  diffraction  from 
edges  has  a  delta  function  character  and  double  diffraction  a 
1  A/t'  -  T  behavior.  It  is  clearly  evident  in  Fig.  5  that  objects  which 
initially  present  identical  geometrical  and  physical  properties  to  the 
incident  transient  excitation  must  also  have  identical  response  wave¬ 
forms  for  an  initial  time  period.  The  impulse  response  waveforms 
for  the  cone,  cone- sphere  and  double  cone  in  Fig.  5  are  identical 
until  a  time  corresponding  to  the  arrival  of  the  incident  impulsive 
wave  at  a  change  in  geometry.  Consider  a  more  explicit  example. 
Compare  two  cone-sphere  targets  (as  in  Fig.  5b)  identical  in  every 
respect  except  for  an  axial  needle  protruding  from  the  sphere  in  one 
case.  For  an  incident  impulsive  wave  along  the  cone  axis  toward  the 
cone  tip,  the  impulse  response  waveforms  of  these  targets  would  be 
identical  for  all  t'  less  than  ^  6.  3t0.  For  times  larger  than  this,  the 
two  response  waveforms  may  differ  markedly.  We  emphasize  the 
fact  that,  as  the  impulsive  wave  moves  over  the  object  in  the  illumi¬ 
nated  region  or  the  transient  currents  excited  by  the  incident  wave 
move  beyond  the  shadow,  a  change  in  the  geometry  or  composition 
of  the  object  cannot  be  anticipated.  The  excitation,  propagating  with 
a  finite  velocity,  must  reach  the  perturbation  before  a  change  in  the 
response  waveform  can  occur.  Thus  if  one  were  attacking  the  scatter¬ 
ing  problem  of  the  perturbed  target,  the  early  portion  of  the  waveform 
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would  be  known  and  the  analysis  could  be  concentrated  on  the  effect  of 
the  perturbing  geometry.  A  question  clearly  exists  as  to  exactly  when 
the  forcing  function,  i.  e.  ,  the  incident  impulsive  wave,  ceases  to 
influence  the  transient  currents  induced  on  the  target.  From  a  geo¬ 
metrical  optics  viewpoint,  the  forcing  function  would  shut  off  at  the 
shadow  boundary.  More  realistically,  the  impulsive  wave  probably 
continues  to  influence  the  induced  current  somewhat  beyond  the 
shadow  boundary  but  confined  to  the  neighborhood  of  the  boundary. 
Therefore  in  the  case  where  a  change  in  the  geometry  or  composition 
occurs  close  to  the  shadow  boundary,  some  uncertainty  would  exist 
in  predicting  exactly  when  the  response  waveform  should  depart  from 
the  unperturbed  case.  This  causality  feature  of  the  response  wave¬ 
form  with  respect  to  a  given  geometrical  feature  of  the  target  is 
extremely  useful  in  devising  estimates  of  the  waveform.  That  is,  in 
most  cases  the  location  in  time  if  not  the  character  of  the  response 
due  to  a  specific  feature  can  be  determined. 

The  most  important  question  in  the  application  of  the  time  do¬ 
main  concepts  discussed  thus  far  in  this  chapter  is:  What  is  a  valid 
general  analytical  model  for  the  time -dependent  response  waveform 
Fj(  t ' )  ,  or  alternatively  for  the  frequency-dependent  phasor  response 
G(  s)  ,  in  a  scattering  problem?  In  the  context  of  this  study,  the 
scattering  object  is  a  two-port  transducer  and  the  two-port  is  a  linear, 
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distributed-constant  system.  Thus  G(  s)  must  be  some  combination 
of  rational  and  transcendental  functions  of  s.  Unfortunately,  the 
literature  on  the  analysis  of  distributed  constant  systems  via  trans¬ 
form  methods  is  quite  limited  compared  to  that  on  lumped  constant 
idealizations.  Two  points  should  be  made;  first,  the  impulse  re¬ 
sponse  of  a  lumped  constant  system  is  only  an  approximation  (at 
sufficiently  high  frequencies  the  system  is  distributed)  ,  but  one 
which  is  sufficiently  good  for  the  frequency  ranges  of  practical  inter¬ 
est,  secondly,  some  circuit  concepts  developed  by  transform  methods 
with  reference  to  lumped  constant  systems  may  have  to  be  modified 
for  application  to  a  distributed  constant  system. 

It  must  be  stressed  that  in  the  context  of  this  study,  what  is 
sought  is  an  approximate  analytical  model  capable  of  producing  a 
reasonable  estimate  of  Fj(t')  or  G(  s)  in  the  range  of  practical  inter¬ 
est.  For  example,  the  dc  and  optical  properties  of  the  scatterer  are 
not  really  of  interest  here;  and  consequently,  the  estimate  of  G(  s) 
need  not  be  valid  at  these  limits  of  the  spectrum.  Such  a  band-limited 
approximation  for  G(  s)  ,  or  alternatively  for  Fj(  t1)  ,  is  further 
justified  when  one  considers  the  effective  bandwidths  of  practical 
radar  systems.  The  question  becomes  therefore:  What  is  a  sufficient¬ 
ly  general  physical  model  of  the  scattering  process  in  such  a  band- 
limited  approximation?  Clearly,  the  necessary  complexity  of  the 
physical  model  depends  to  some  extent  on  the  specific  scattering 
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object  and  one  cannot  rigorously  justify  the  assumption  of  any  particu¬ 
lar  model  for  the  general  object.  However,  the  representative  wave¬ 
forms  shown  in  Figs.  4  and  5  do  illustrate  some  of  the  characteristics 
the  general  model  must  be  able  to  reproduce.  Specifically,  it  can  be 
seen  that  over  limited  time  ranges,  the  model  must  be  capable  of 
generating  an  increasing  rather  than  decreasing  function  with  time 
and  also  admit  to  at  least  approximations  of  various  types  of  singu¬ 
larities.  Finally,  for  moderately  large  times  the  model  must  produce 
a  rapid  damping. 

Consider  first  the  large  time  behavior  of  the  impulse  response 
waveform.  It  is  known[  24]  that  in  any  passive  or  stable  active 
lumped-constant  network,  the  force-free  behavior  decays  exponential- 
ly  so  that  a  short  time  following  the  transient  interval  only  the  forced 
behavior  survives.  In  the  present  case  of  a  distributed  constant  sys¬ 
tem,  this  generalization  cannot  be  made  precisely  because  of  the 
presence  of  the  transcendental  functions  of  s.  However,  given  suffi¬ 
cient  elements,  one  can  synthesize  an  approximation  to  the  distributed 
constant  system  using  lumped  elements.  Clearly  this  lumped  constant 
approximation,  regardless  of  the  number  of  elements  needed  to  secure 
a  reasonable  estimate,  must  decay  exponentially  for  large  time.  For  an 

*  In  the  general  case  one  has  a  summation  of  exponentials,  but  for 
large  time  the  minimum  decay  rate  dominates  and  the  response 
approaches  an  exponential  decay. 
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approximate  model  for  Fj(  t* )  an  exponential  decay  for  large  time  is 
assumed,  i.  e.  ,  beyond  some  finite  value  of  t'  which  will  vary  with 
the  scatterer,  an  exponential  decay  of  the  waveform  is  enforced. 
Physically,  this  exponential  decay  is  associated  with  radiation  and 
dissipation  effects.  If  the  inverse  transformation  of  Eq.  (  5)  is 
written  as 

i  j® 

(  25)  Fj(  t1)  =  -  \  G(  s)  et  8  ds  y  >  -o 

2irjc  J 

V-J® 

then  a  is  the  minimum  exponential  decay  rate  of  the  impulse  response 
waveform  for  large  values  of  t'  .  The  results  in  Figs.  4  and  5  indi¬ 
cate  that  a  is  reasonably  large  for  practical  scatterers. 

It  is  also  known(  24]  that  for  lumped  constant  systems  with 
simple  poles  the  most  general  representation  of  the  impulse  response 

waveform  is  a  linear  combination  of  terms  of  the  type 

N 

(26)  f(t)  =2  Anean‘u(t) 

n=l 

This  in  fact  is  the  most  general  time  function  having  a  rational 
Laplace  transform  when  the  Re(  an)  <  0.  It  is  immediately  clear 
that  one  would  require  N  large  to  approximate  the  response  wave¬ 
forms  in  Figs.  4  and  5  with  a  function  of  this  type,  except  for  large 
t'  .  However,  if  the  real  part  of  the  an  are  unrestricted  and  a  delay 
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factor  is  introduced  then  a  linear  combination  of  terms  of  this  type 
should  be  capable  of  achieving  reasonable  estimates  of  the  time  re¬ 
sponse  waveforms  with  just  a  few  terms.  A  possible  analytical  model 
for  the  impulse  response  waveform  is 


M 


(27) 
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k=l  n=l 
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The  model  in  Eq.  (  27)  consists  of  weighted  delta  function  singularities 
at  selected  instants  plus  sections  of  sums  of  exponential  functions  and 
a  final  exponential  decay.  Note  that  there  are  no  restrictions  on  the 
real  part  of  the  b^  except  in  the  final  time  interval,  i.  e.  ,  t'  >  . 

In  this  final  time  interval 


(28)  Re(bnK)  <  0 

to  ensure  an  exponential  decay  of  the  response  waveform  for  large 
time. 

Excluding  the  delta  function  singularities  and  the  exponential 
decay  for  large  time,  the  analytical  model  in  Eq.  (27)  is  a  super¬ 
position  of  bounded,  time -limited  waveforms  in  sub-intervals  or 
sections.  The  number  of  sub-intervals  into  which  the  response 
waveform  is  separated  is  finite  for  any  scatterer,  and  judging  from 
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the  response  waveforms  in  Figs.  4  and  5  will  usually  be  no  more  than 
3  or  4.  The  point  is  that  the  number  of  sub-intervals  (K-l) ,  their 
relative  location  ( T^)  and  duration  (T^+i  -  T^)  is  selected  from  an 
examination  of  the  geometry  of  the  target.  Locations  on  the  target 
where  changes  in  the  curvature  or  the  derivatives  of  curvature,  or 
the  composition  of  the  scatterer  occur  are  natural  join  points  of  the 
sub-intervals.  With  this  type  of  model  it  is  clear  that  both  experi¬ 
ence  and  a  knowledge  of  the  known  response  waveforms  for  various 
objects  ( a  catalog  of  known  response  waveforms)  is  an  indispensable 
aid  in  making  a  judicious  choice  of  the  model  for  a  given  target. 

It  is  not  intended  to  imply  in  the  compactly  written  analytical 
model  in  Eq.  (  27)  that  the  same  number  of  exponential  terms  neces¬ 
sarily  be  used  in  each  sub-interval.  The  number  of  terms  required 
in  each  sub-interval  will  depend  on  the  complexity  of  the  waveform, 
primarily  at  the  end  points  of  the  interval.  It  is  intended  that  the 
character  of  the  response  at  the  joins  of  the  sub-intervals  be  deter¬ 
mined  from  various  approximate  or  asymptotic  theories  and  the  ge¬ 
ometry  of  the  shape.  In  this  regard  it  is  evident  that  the  analytical 
model  in  Eq.  (  27)  cannot,  in  its  present  form,  faithfully  reconstruct 
singularities  such  as  those  corresponding  to  double  diffraction  (  Fig. 
5a).  However,  such  cusp  type  singularities  are  usually  predicted  by 
an  asymptotic  theory  and  are  associated  with  abrupt  discontinuities  in 
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the  curvature  of  an  object.  Such  discontinuities  are  mathematical 
idealizations  for  real  targets  which  exhibit  less  perfect  curvature 
transitions.  It  is  contended  that  a  practical  interpretation  of  the 
singularities  resulting  from  such  discontinuities  is  a  bounded  vari¬ 
ation  of  the  response  waveform.  Therefore  an  approximation  to  the 
singularity  can  be  obtained  with  the  analytical  model  proposed  in 
Eq.  (27). 

Neglecting  the  delta  function  singularities,  this  analytical 
model  of  the  impulse  response  waveform  is  a  superposition  of  bound¬ 
ed,  time-limited  waveforms  with  various  delays  plus  a  final  expo¬ 
nential  decay.  In  the  s  domain,  the  model  corresponds  to  a  combi¬ 
nation  of  entire  functions[  25]  ,  i.  e.  ,  functions  which  have  no  singu¬ 
larities  in  the  finite  portion  of  the  plane;  and  rational  functions. 

Such  a  function,  i.  e.  ,  one  made  up  of  entire  and  rational  functions, 
is  known  as  meromorphic  in  the  unextended  plane.  The  analytical 

model  for  the  phasor  response  is 
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The  analytical  time  domain  model  in  Eq.  (  27)  and  the  corre¬ 
sponding  phasor  response  in  Eq.  (  29)  have  not  been  based  on  an 
assumed  physical  model  for  the  distributed  constant  system  but  rather 
on  the  known  impulse  response  waveforms  of  certain  representative 
scatterers.  It  is  clearly  not  imperative  that  the  analytical  model  be 
linked  to  a  particular  physical  model;  but  such  a  correlation  is  advan¬ 
tageous  as  a  partial  justification  of  the  model.  Franks[26]  has  shown 

that  the  finite  impulse  sequence 

M 

(30)  f(t)  =2  AmMt-W 

m=l 

given  as  the  first  part  of  Eq.  (  27)  corresponds  to  the  impulse  re¬ 
sponse  of  a  network  consisting  of  delay  lines  lumped  loaded  with  re¬ 
sistors.  He  has  also  shown  that  the  impulse  response  of  delay  lines 
with  general  lump  loading  can  have  the  bounded,  time-limited 
behavior  corresponding  to  entire  functions  in  the  s  domain.  Franks 
suggested  a  segmental  impulse  response  involving  an  exponential  sum 
approximation  in  each  interval  similar  to  the  second  part  of  Eq.  (  27) 
for  such  a  network.  In  his  case  however,  it  was  stipulated  that  the 
real  part  of  the  bj^  be  less  than  zero.  However,  the  purpose  of 
Frank's  approximation  was  for  the  realization  of  a  segmental  impulse 
response  which  consisted  of  segments  of  time  functions  having  rational 
Laplace  transforms  between  critical  break  points  corresponding  to 
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lengths  of  delay  line.  If  one  replaces  the  delay  line  with  finite  lengths 
of  uniform  transmission  line,  then  the  segments  of  time  functions  will 
no  longer  have  rational  transforms.  Th general  representation  in 
Eq.  (  26)  for  time  functions  with  rational  transforms  is  no  longer  ade¬ 
quate  and  a  less  restrictive  representation  is  required.  A  possible 
interpretation  of  the  physical  model  corresponding  to  the  approximate 
analytical  representation  in  Eq.  (  27)  consists  of  a  finite  length  of 
uniform  transmission  line  lump  loaded  at  various  intervals.  As  a 
practical  physical  model  to  approximate  the  scattering  of  an  electro¬ 
magnetic  field  by  a  finite  three  dimensional  object,  such  a  represen¬ 
tation  does  not  seem  unreasonable. 

Some  final  comments  on  the  analytical  model  in  Eq.  (  27)  are 
necessary.  The  impulse  response  waveform  Fj(t')  must  be  real, 
therefore  a  necessary  relationship  between  the  coefficients  in  the 
exponential  sums  exists.  If  the  bnk  are  complex,  they  must  exist 
as  conjugate  pairs  and  the  corresponding  Bnk  are  aal  in  magnitude. 
If  these  Bnk  have  the  same  sign  they  are  real,  :i  they  have  the  op¬ 
posite  sign  they  are  imaginary.  If  the  bnk  are  purely  real  then  the 
Bnk  are  also  purely  real.  In  some  cases  it  may  prove  advantageous 
in  terms  of  determining  the  b,*  and  in  a  particular  interval  to 
write  the  exponential  in  a  shifted  form 
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(31)  Bnkeb»k<t'-Tk> 

in  order  to  isolate  when  t1  equals  T^. 

At  best,  the  preceeding  discussion  has  perhaps  achieved  a  rea¬ 
sonable  rationalization  for  the  use  of  the  analytical  model.  A  rigorous 
demonstration  that  the  analytical  model  has  the  physical  interpretation 
suggested  has  not  been  given,  nor  has  it  been  demonstrated  that  these 
models  are  necessarily  a  valid  approximate  represents  .on  of  an 
actual  scattering  problem.  At  this  point  in  the  development  of  impulse 
response  methods  in  scattering  problems,  however,  it  is  felt  that  such 
rigorous  demonstrations  are  not  necessary.  There  is  evidence[  13] 
that  the  analytical  model  in  Eq.  (  27)  yields  adequate  representations 
of  the  response  waveforms  for  certain  objects.  That  is,  the  published 
estimates  of  the  response  waveforms  for  particular  targets!  13]  , 
which  have  been  verified  by  comparisons  of  the  corresponding  phasor 
responses  with  measured  or  calculated  data,  were  achieved  with  the 
implicit  assumption  of  the  analytical  model  suggested  here.  Thus 
while  a  more  rigorous  basis  for  the  analytical  model  is  desirable, 
it  is  not  a  prerequisite  for  further  applications  of  the  impulse  re¬ 


sponse  approach. 


43 


C.  Analytical  and  Experimental 
Approximations  of  Impulse 
and  Transient  Response 
Waveforms 

The  purpose  of  this  Section  is  to  discuss  in  a  general  manner 
the  development  of  approximate  impulse  and  transient  response  wave¬ 
forms.  To  place  the  discussion  in  proper  perspective,  it  is  perhaps 
advisable  to  begin  with  some  comments  regarding  electromagnetic 
scattering  studies.  If  the  ultimate  goal  of  such  research  is  an 
"understanding"  of  the  scattering  by  various  objects,  then  the  most 
direct  route  to  this  "understanding"  lies  via  the  development  of  the 
time  domain  approach  rather  than  the  development  of  theoretical 
solutions  per-se.  The  point  is  that  the  ultimate  refinement  of  a 
given  theoretical  solution  lends  little  insight  to  the  general  problem. 
Using  the  time  domain  approach,  it  is  possible  to  secure  estimates 
of  the  response  of  a  variety  of  object  shapes  in  a  limited  amount  of 
time.  A  most  important  consideration  is  the  applicability  of  experi¬ 
mental  data.  Despite  the  vast  literature  on  electromagnetic  scatter- 
ing[  1]  ,  one  is  hard  pressed  to  find  examples  wh^re  experimental 
data  have  been  used  in  a  direct  way  to  help  construct  a  theoretical 
estimate.  The  handbook  approach,  followed  almost  without  exception, 
has  been  to  devise  an  estimate  of  the  phasor  response  based  purely  on 
theoretical  considerations  and  then  compare  calculations  from  this 
model  with  experimental  data.  Our  criticism  is  not  of  experimental 
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data,  but  rather  of  the  fact  that  such  data  have  been  generated  haphaz¬ 
ardly  without  any  useful  integration  of  the  results.  With  the  ever  in¬ 
creasing  complexity  and  number  of  objects  of  practical  interest,  the 
utilization  of  experimental  data  must  be  extended  beyond  simply  a 
confirmation  of  theoretical  estimates  or  the  compilation  of  patterns  on 
a  specific  shape. 

The  time  domain  approach  appears  to  offer  the  first  practical 
way  in  which  to  build  one's  knowledge  and  understanding  of  scattering 
by  the  systematic  improvement  of  a  basic  model.  Thus,  experimental 
measurements  as  well  as  computer  derived  solutions[  27]  can  be  used 
in  a  direct  way  to  refine  the  theoretical  model  discussed  in  the  previous 
section.  It  is  to  this  systematic  procedure  that  the  discussion  in  this 
section  is  directed. 

It  is  postulated  that  the  fundamental  building  block  in  the  develop¬ 
ment  of  a  general  approach  for  obtaining  the  response  waveforms  of  an 
object  is  the  ramp  response  waveform  Fj^(t').  There  are  several 
reasons  which  make  an  initial  estimate  of  the  ramp  response  wave¬ 
form  a  logical  starting  point.  It  was  pointed  out  in  the  previous  section 
that  one  could  estimate  the  ramp  response  waveform  using  fewer  pa¬ 
rameters  than  would  be  required  for  estimates  of  the  corresponding 
step  or  impulse  response  waveforms.  If  the  estimated  ramp  response 
waveform  has  a  finite  integral  over  all  t'  then  the  first  two  moment 
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conditions  are  automatically  satisfied.  Forcing  this  finite  integral  to 
be  proportional  to  the  Rayleigh  coefficient  (  Eq.  (15))  then  satisfies 
the  first  three  moment  conditions.  Finally,  it  has  been  shown[  13) 
that  whereas  many  values  of  the  phasor  response  are  required  to 
synthesize  an  estimate  of  the  impulse  response  waveform,  as  few  as 
five  values  of  the  phasor  response  can  be  used  to  generate  an  estimate 
of  the  ramp  response  waveform.  The  required  input  data  are  low  fre¬ 
quency  harmonics  with  the  fundamental  wavelength  approximately  10 
times  the  maximum  linear  dimension  of  the  object.  Note  that  the  neces¬ 
sary  phasor  response  data  could  be  obtained  from  experimental  measure¬ 
ments  as  well  as  theoretical  calculations;  consequently,  an  experimental¬ 
ly  derived  first  estimate  of  the  ramp  response  waveform  is  also  feasi- 
ble[  13]  .  Regardless  of  the  method  used  to  generate  this  first  estimate, 
the  corresponding  estimate  of  the  phasor  response  will  be  a  low-fre¬ 
quency,  i.  e.  ,  Rayleigh  and  low  resonance  region  approximation. 

In  most  previous  applications  of  the  time  response  concepts, 
direct  estimates  of  the  impulse  response  waveform  have  been  attempted 
employing  segmental  exponential  sum  approximations  similar  to  Eq. 

(27).  In  specific  cases,  sinusoids  or  polynomials  in  t'  have  been 
used  in  place  of  an  exponential  sum  in  some  sub-intervals.  In  general, 
the  physical  optics  approximation  is  used  to  predict  the  short  time 
behavior  of  the  response  waveform.  Later  portions  of  the  waveform 
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are  selected  either  from  a  knowledge  of  the  character  of  the  response 
corresponding  to  specific  geometrical  features  of  the  object  or  simply 
from  a  rough  guess  as  to  its  probable  form.  These  various  pieces  of 
the  waveform  are  then  joined  with  sufficient  undetermined  parameters 
to  permit  the  known  moment  equations  to  be  satisfied.  Exactly  this 
type  of  approach  is  used  in  Chapter  II  to  obtain  an  estimate  of  the 
impulse  response  waveform  of  a  prolate  spheroid  target;  the  extension 
being  to  a  general  orientation  of  the  spheroid  and  an  arbitrary  linear 
polarization  of  the  incident  plane  wave. 

The  approach  outlined  above  for  obtaining  estimates  of  the 
impulse  response  waveform  of  an  object  is  clearly  inadequate.  The 
fact  that  such  a  procedure  has  led  to  significant  contributions  to 
electromagnetic  scattering  theory  is  probably  the  strongest  possible 
recommendation  for  further  study  in  the  time  domain.  The  most 
obvious  difficulty  with  the  method  above  is  that  a  successful  estimate 
of  the  impulse  response  waveform  is  almost  wholly  dependent  on  the 
ingenuity  of  the  analysist  in  deducing  the  character  of  the  response 
waveform  for  other  than  short  times.  The  moment  conditions  can 
also  lead  to  problems  in  that  enforcing  the  moments  results  in  a 
system  of  simultaneous  nonlinear  equations.  Even  a  very  simple 
model  for  the  response  waveform  can  lead  to  an  extremely  difficult 
system  of  equations.  Thus  in  many  cases  it  is  necessary  to  devise 


47 


a  model  such  that  a  soluble  system  of  equations  is  obtained,  and  this 
usually  requires  some  sacrifice  in  the  desired  character  of  the  wave¬ 
form.  The  approach  also  lacks  the  property  of  finality  in  that  if  an 
additional  exponential  term  is  added  in  any  one  sub-interval,  then  all 
of  the  coefficients  must  be  redetermined.  The  final  criticism  is  that 
there  is  no  established  procedure  for  selecting  the  distribution  of  pa¬ 
rameters  in  the  various  sub-intervals  which  are  to  be  determined 
from  known  conditions  or  the  waveform.  It  seems  apparent  that  a 
need  exists  for  studies  on  techniques  for  securing  estimates  of  the 
response  waveforms  of  distributed  constant  systems.  It  is  not  in¬ 
tended  to  pursue  such  an  investigation  in  this  research  beyond  the 
point  of  indicating  some  possible  lines  of  attack.  The  reader  should 
note  that  one  feature  of  the  segmental  impulse  response  waveform 
will  probably  be  retained,  ultimately,  in  any  derived  estimate  of  the 
response  waveform  of  a  given  object.  The  singularities,  disconti¬ 
nuities,  and  peaks  in  the  impulse  response  waveform  are  directly 
associated  with  the  arrival  of  reflected  and  diffracted  fronts  from 
specific  geometrical  and  physical  features  of  the  object.  Since  the 
arrival  time  of  these  fronts  can  usually  be  accurately  predicted, 
procedures  for  estimating  the  response  waveforms  of  the  object 
must,  at  one  point,  successfully  exploit  this  knowledge. 
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In  considering  alternative  methods  for  estimating  the  response 
waveforms,  one  is  naturally  led  to  a  consideration  of  the  techniques 
developed  for  time  domain  synthesis  and  time  domain  identification 
in  network  and  control  system  studies.  The  literature,  as  was  noted 
in  the  case  of  linear  system  analysis  via  transform  methods,  is  pre¬ 
dominantly  devoted  to  a  consideration  of  lumped  constant  idealizations. 
Horowitz[  28]  has  discussed  a  number  of  these  techniques.  Two  main 
points  are  noted;  emphasis  is  on  expansion  of  the  impulse  response 
waveform  in  a  finite  series  of  orthogonal  functions,  and  an  error 
criterion,  usually  the  integral  square -error  or  the  weighted  integral 
square-error  is  used  to  judge  the  accuracy  of  the  approximation. 

Both  time  domain  synthesis  and  time  domain  identification  are  con¬ 
cerned  with  the  following  problem.  How  does  one  find  a  frequency- 
dependent  response  Ha(  s)  whose  inverse  transform  approximates  a 
given  time  function  f(  t)  ?  As  previously  stated,  Ha(  s)  is  usually 
assumed  to  be  a  rational  function  of  s.  Clearly,  the  problem  at  hand 
is  somewhat  different.  Simply  stated  it  is;  how  does  one  find  a  time- 
dependent  impulse  response  waveform  whose  transform  approximates 
the  phasor  response  of  a  system  over  the  range  of  practical  frequen¬ 
cies,  given  certain  conditions  on  the  waveform  and  estimates  of  its 
value  at  particular  times?  There  are  common  areas  of  the  problems. 
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however,  which  enable  one  to  consider  application  of  some  of  the 
techniques  described  in  the  literature. 

It  was  concluded  previously  that  the  initial  estimate  of  the 
time-dependent  waveform  should  be  made  on  the  ramp  response 
waveform  of  the  target.  From  a  mathematical  viewpoint  it  would 
seem  that  the  ramp  response  waveform  should  be  expanded  in  a 
set  of  functions,  complete  and  orthogonal  on  the  range  0  to  ®  . 

The  reason  for  this  is  that  if  it  is  assumed  that  the  ramp  response 
waveform  Fp(  t* )  is  square  integrable 

GO 

(32)  ^  |  Fr(  f)  |2  df  <  ® 

0 


(35) 
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In  a  great  many  cases  it  is  reasonable  to  expect  Fp(  t')  to  be  continu¬ 
ous,  since  a  jump  discontinuity  in  Fp(t')  would  mean  a  doublet  in  the 
impulse  response  waveform.  The  coefficients  in  such  an  expansion 
have  the  property  of  finality,  an  obvious  advantage.  The  problem  of 
course  is  in  obtaining  the  coefficients  An  in  the  expansion.  Formally, 
if  the  4>n  are  assumed  to  be  an  orthonormal  set  then 

CO 

(36)  An  =  )  Fr<  f)  4>n(t')dt' 

0 

Note  that  one  could  therefore  express  estimates  of  the  ramp  response 
waveform  obtained  by  any  method  in  the  form  of  Eq.  (  33). 

Kautz[  29]  has  described  the  conversion  of  polynomial  sets 
which  are  orthogonal  over  a  finite  interval  into  exponential  sets  or¬ 
thogonal  over  the  semi-infinite  interval  0  <  x  <  ®.  An  advantage  of 
this  technique  is  that  whereas  simply  transforming  sets  orthogonal 
over  a  finite  interval  into  sets  orthogonal  over  the  semi-infinite  inter¬ 
val  leads  to  poles  in  the  s  plane  equally  spaced  along  the  negative  real 

Jr 

axis,  the  exponential  set  permits  an  arbitrary  location  of  the  poles. 
Assuming  a  set  4>n,  orthonormal  over  the  semi-infinite  interval  Kautz 
forms 


*  An  exception  is  the  Laguerre  polynomials,  orthogonal  over  the 
semi-infinite  interval  with  weight  e"^at  ,  which  lead  to  a  simple 
multiple  order  pole  at  s  =  -a[  29]  . 
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{  37)  <Mf)  =  ReCn  e'8,t' 

4>z(t')  =  ReCzi  e“Sl 

4>3  ( t')  =  ReC3i  e_Slt'  +  ReC32  e’®2*' 

4>4  (  f)  =  ReC4i  e"8lt'  +  ReC52  e”®2*'  +  ReC53  e~8*V 
etc.  , 

where  the  are  arbitrary.  From  the  orthonormal  con¬ 

dition  on  the  <j>n  it  is  shown  that  the  coefficients  CVJj  can  be  determined 
and  relatively  simple  computational  techniques  are  given  involving 
partial  fraction  expansions  of  the  Laplace  transform  of  the  <|>n.  From 

Eq.  (  36)  ,  taking  as  an  example  the  third  coefficient 

00  00 

(38)  A3  =  Re ^  C31  e’Slt'  FR(t')dt'  +  Re^  C32  e"®2*'  FR(  t')  dt'  . 

0  0 

But  from  Eq.  (  5)  , 

00 

r  _st.  1  G(S) 

(39)  )  e  FR(t')dt'=  -  — 

0 

Thus  Eq.  (  38)  becomes 

C3i  G(  si )  C32  G(  s2) 

(40)  A3  =  Re  -  — : —  +  Re  - - - 

c  sj  c  sf 


or  in  general[  29]  ,  An  is  replaced  by  the  sum 
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(41) 


n 


i  =  l 


G  <■!> 


Now  if  ail  the  poles  are  chosen  to  lie  equally  spaced  along  the 
jw  axis  (  Sfe  =  jPk)  then  this  configuration  corresponds  to  a  Fourier 
series  in  the  time  domain,  and  each  orthonormal  function  of  the  set 
consists  of  a  single  sine  or  cosine  term.  Note  that  G(  jPjJ  is  the 
phasor  response  of  the  scatterer  at  the  angular  frequency  P^,  a 
measurable  quantity.  For  the  special  pole  configuration  leading  to 
a  Fourier  series,  the  result  is  simply  an  alternative  Fourier  synthe¬ 
sis  procedure.  However,  the  pole  locations  (s^)  are  arbitrary  and 
one  could  choose  unequally  spaced  imaginary  poles  or  any  desired 
combination  of  real  and  complex  poles.  For  example,  the  above 
procedure  could  be  used  when  calculated  or  measured  values  of  the 
phasor  response  were  available  at  nonharmonically  related  frequen- 
cies.  When  only  imaginary  poles  are  utilized  Kautz[  29]  (for 
equally  spaced  poles)  suggests  a  predistortion  of  the  response  wave¬ 
form 

+at' 

(42)  FR(  t' )  =  e  F R(  t ' ) 


*  A  procedure  for  this  situation  was  also  suggested  in  Reference  30. 
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That  is,  the  ramp  response  waveform  is  divided  by  e  *  before  ex- 

at1 

pension  and  then  multiplied  by  e  afterwards.  This  has  the  effect 
of  multiplying  the  error  bound  as  a  function  of  time  by  e”at  [  29]. 

Thus  far  it  is  clear  that  an  analytical  estimate  of  the  ramp  response 
waveform  can  be  put  in  the  form  of  an  orthonormal  expansion  and  that 
an  initial  orthonormal  expansion  estimate  of  the  ramp  response  wave¬ 
form  can  be  obtained  from  calculated  or  measured  data  on  the  phasor 
response.  It  is  also  known[  13]  that  such  data  should  be  predominant¬ 
ly  low-frequency. 

An  attractive  feature  of  an  orthonormal  expansion  is  that  it 

automatically  minimizes  the  rms  error[  29]  .  That  is,  if  Fp(t')  is 

the  actual  ramp  response  waveform  and  F^(  t')  an  orthonormal 

approximation  such  as  Eq.  (  33)  then 

00 

(43)  ‘2  =  \[FR(t')  -  FRA(t')]*dt' 
becomes 

00  N 

(44)  «2  =  y  Fr(  t1)  dt'  *2  An 

0  n=0 
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where  G(  •)  corresponds  to  FR(  t')  and  G^(  s)  to  Fp^(t').  Equation 
(  45)  may  also  be  written  as 


where 


A(jw)  =  4r(jw)  +jAi(jw)  =  G(jw)  -  Ga(  jw ) 


Regardless  of  how  the  approximation  of  the  ramp  response  is  obtained, 
an  estimate  of  the  rms  error  can  be  obtained  from  phasor  response 
data.  How  good  an  estimate  of  the  rms  error  one  can  obtain  clearly 
depends  on  the  extent  and  nature  of  the  phasor  response  data.  In  any 
practical  situation,  at  most  a  band  limited  knowledge  of  the  phasor 
response  may  be  known  and  Eq.  (  46)  becomes 
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where  u>f  is  the  cut-off  frequency  of  the  band-limited  phasor  response. 
Note  that  if  in  the  range  0  <L  w  <  wf,  the  maximum  value  of  A(  jwj/w2 
occurs  at  u>p  then 
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Both  Fp(t')  and  FRA(  t')  must  be  causal  time  functions,  therefore 

)|g 

Eq.  (  47)  may  be  written  as[  20] 


(49) 
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and  in  the  notation  of  Eq.  (  48) , 
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In  a  great  many  practical  cases,  only  data  on  the  radar  cross  section 
are  available.  One  possible  method  for  securing  an  error  estimate, 
not  of  the  waveforms  themselves  but  of  the  waveforms  convolved  with 
themselves,  can  be  given. 

It  can  be  shown  that  if  a  ramp  response  waveform  is  convolved 
with  itself,  i.  e.  ,  the  autocorrelation  of  the  ramp  waveform,  the  re¬ 
sultant  waveform  is  given  by 


(  51)  Frr(  t')  = 


m 

G(  jw)  I  •' 

-  1 

jw) 

U)2  1 

i 
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*  Equation  (  49)  actually  only  holds  for  Wf  *  ®,  but  for  the  present 
case,  i.  e.  ,  G(  s)/s2,  it  is  felt  to  be  a  reasonable  approximation 
for  finite  wf. 


Thus  forming  the  error  measure 
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Data  on  the  phasor  response  or  radar  cross  section  are  usually  given 
as  a  function  of  an  electrical  length,  say  ki  .  The  equivalent  expres¬ 


sion  for  Eq.  (  53)  is 


(54) 
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f  r<r(Zki)  -  <rA(2ki) 

0  L  (2k i)4 


d(kf  ) 


where  tQ  =  2f/c  is  set  equal  to  unity.  Note  that  the  error  estimates 
given  above  are  applicable  to  any  approximation  of  the  ramp  response 
waveform.  They  are  not  confined  to  orthogonal  expansions  of  the  re¬ 
sponse.  Also,  similar  expressions  involving  the  step  and  impulse  re¬ 
sponse  waveforms  of  the  target  can  be  obtained  with  the  restriction 
that  these  waveforms  be  real,  causal,  and  square  integrable.  It  is 
possible  therefore  to  obtain  an  estimate  of  the  rms  error  for  a  given 
approximation  and  deduce  from  this  when  further  modifications  of  the 
model  are  not  justified.  Such  a  gauge  is  important  in  that  there  are 
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really  three  approximation  problems  to  be  solved,  sequentially  the 
ramp,  step  and  impulse  response  waveforms  of  the  object,  and  one 

must  economize  the  treatment  of  the  preliminary  steps. 

>1 

Simply  stated,  the  suggested  approach  to  the  approximation 
problem  is  as  follows:  From  a  combination  of  theoretical  estimates 
and  experimental  data,  a  basic  model  of  the  ramp  response  waveform 
is  obtained.  This  is  a  low-frequency  model  and  the  number  of  pa¬ 
rameters  involved  should  be  minimized  consistent  with  obtaining  an 
rms  error  of  possibly  5  to  10  percent.  The  model  is  then  interpreted 
as  a  step  response  waveform  (  differentiation)  ,  additional  parameters 
added  and  an  rms  error  of  again  5  to  10  percent  achieved.  Note  at 
this  point  that  two  simultaneous  conditions  must  be  satisfied  by  the 
added  parameters.  A  specified  minimum  error  in  the  step  waveform 
is  sought  but  at  the  same  time  the  ramp  response  error  must  be 
maintained.  Proceeding  to  the  impulse  response  waveform,  addition¬ 
al  parameters  are  once  again  included.  At  this  point  it  is  desired  to 
build  into  the  model  all  available  information.  That  is,  approximate 
or  asymptotic  high  frequency  estimates  of  the  reflected  and  diffracted 
fields,  natural  points  of  discontinuity  or  rapid  change  of  the  response 
corresponding  to  geometrical  or  physical  features  of  the  object,  high 
frequency  experimental  data,  etc.  There  are  now  three  conditions  to 
satisfy,  the  ramp,  step  and  impulse  response  error. 
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Realistically,  it  must  be  recognized  that  how  meaningful  these 
errors  are  depends  upon  the  extent  of  the  phasor  response  data  avail¬ 
able.  A  band- limited  estimate  of  the  area  beneath  the  real  and  im¬ 
aginary  parts  of  the  phasor  response  as  a  function  of  frequency  clearly 
requires  much  more  data  than  to  estimate  the  area  beneath  these 
same  curves  divided  by  the  frequency  squared.  The  point  is  that  a 
reasonable  estimate  of  the  integral  in  Eq.  (  49)  might  be  obtained 
from  a  few  properly  distributed  values  of  Ar(  ju>)  or  Aj(  ju ) ,  but  the 
corresponding  integrals  for  the  step  and  impulse  error  would  require 
an  increasing  density  of  known  values.  The  upper  bound  on  the  error 
given  in  Eq.  (  50)  is  somewhat  conservative.  Lacking  sufficient  data 
to  estimate  Eq.  (  49)  directly,  a  more  reasonable  estimate  than  Eq. 

(  50)  would  be  to  assume  the  average  Ar(  jw)/  w2  or  A*(  jwj/u2  to  be 
independent  of  frequency.  That  is,  given  Q  values  of  Ar(jw)/w2,  iet 
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Since  the  ramp  waveform  is  a  low  frequency  estimate,  a  weighted 
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might  also  be  used.  The  choice  depends  upon  both  the  cut-off  frequen¬ 
cy  and  the  distribution  of  the  data  within  this  bandwidth. 

If  one  is  to  use  an  orthonormal  expansion  to  estimate  the  re¬ 
sponse  waveform,  the  major  problem  is  to  select  the  proper  set  of 
orthogonal  functions.  Strictly  speaking,  this  involves  the  pole  lo¬ 
cations  of  the  orthonormal  expansion  and  hence  the  natural  resonances 
of  the  transient.  That  is,  the  optimum  location  of  the  poles  is  at  the 
natural  resonances  of  the  scattering  object.  The  flexibility  of  the 
orthogonal  exponentials  of  Kautz[  29]  ,  described  earlier,  have  an 
advantage  in  this  regard  pince  the  pole  locations  may  be  arbitrarily 
specified.  However,  unless  all  the  specified  poles  have  negative, 
non-zero  real  parts,  the  orthogonal  functions  do  not  have  the  neces¬ 
sary  decaying  behavior  with  large  time  and  a  predistortion  would  be 
required.  The  set  of  orthonormal  Laguerre  polynomials  have  a 
decaying  behavior  for  large  time,  but  correspond  to  a  single  multiple 
order  real  pole[  29]  .  It  can  be  shown  that  if  a  given  time  function  is 
expanded  in  a  set  of  orthonormal  Laguerre  polynomials,  the  coef¬ 
ficients  in  the  expansion  are  a  linear  combination  of  the  coefficients 
of  a  power  series  expansion  of  the  corresponding  frequency  domain 
function  about  the  real  pole.  Thus  such  an  expansion  eliminates  the 
possibility  of  directly  determining  the  expansion  coefficients  in  terms 
of  measured  phasor  response  data.  It  is  not  possible,  at  this  point, 
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to  recommend  one  particular  set  of  orthogonal  functions  for  expansion 
of  the  response  waveforms.  There  are  numerous  possibilities!  29] , 
and  perhaps  the  most  expedient  procedure  would  be  to  compare  ex¬ 
pansions  for  several  known  ramp  response  waveforms!  13]  .  Note 
that  the  objective  is  not  a  selected  set  for  each  individual  scatterer 
but  rather  one  particular  set  which  has  a  reasonable  applicability  to 
numerous  objects.  Assuming  that  a  particular  orthonormal  set  has 
been  selected,  logical  procedures  for  estimating  the  ramp  response 
waveforms  of  an  object  can  be  postulated. 

An  initial  estimate  of  the  ramp  response  waveform  can  be  se¬ 
cured  in  several  ways;  the  Fourier  synthesis  procedure  using  experi¬ 
mental  data  and  an  analytical  approach  combining  the  physical  optics 
approximation  and  the  second  moment  condition  as  mentioned  previ¬ 
ously.  It  is  not  expedient  to  adopt  the  latter  approach,  precisely 
because  of  the  necessity  for  selecting  a  specific  functional  form  for 
the  model.  There  is  no  reason,  however,  why  one  cannot  simply 
sketch  a  curve  for  the  ramp  waveform  utilizing  exactly  the  same 
information,  i.  e.  ,  physical  optics,  the  geometry  of  the  object,  and 
the  second  moment.  Such  a  possibility  has  been  suggested  previ¬ 
ously!  13]  .  Envisoned  is  some  type  of  analog  interplay  with  a  digital 
computer  whereby  one  could  very  rapidly,  by  trial  and  error,  secure 
a  ramp  waveform  similar  to  that  predicted  by  physical  optics  over 
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most  of  the  illuminated  region  of  the  target*  which  had  the  proper  net 
area*  and  which  agreed  with  experimental  data.  It  must  be  stressed 
that  the  anticipated  estimate  of  the  ramp  waveform  is  a  simple*  con- 

$ 

tinuous  curve  with  at  most  possibly  three  significant  zero  crossings. 
Thus  either  from  a  Fourier  synthesis  of  experimental  data  or  the 
above  described  sketched  estimate*  a  first  model  for  the  ramp  re¬ 
sponse  waveform  is  secured.  The  model  is  simply  a  curve  in  graphi¬ 
cal  form,  and  is  not  described  analytically.  This  model  can  now  be 
approximated  by  an  orthonormal  expansion*  the  integrations  in  L'q. 

(  36)  being  approximated  using  digital  computer  techniques.  The 
fact  that  a  time-limited  model  for  the  ramp  response  waveform  is 
used  makes  this  approach  feasible.  One  also  sees  the  advantage  of 
using  one  orthonormal  set  for  at  least  a  class  of  objects;  the  ortho¬ 
normal  functions  need  only  to  be  calculated  once*  over  an  initial  time 
range,  at  suitable  increments. 

The  second  moment  condition  introduces  a  constraint  on  the 
coefficients.  Let  $n(  s)  be  the  Laplace  transform  of  4>n(  t)  *  from 
Eqs.  (  1  5)  and  (  36) 


*  A  complicated  (  in  terms  of  geometrical  or  physical  properties) 
target  may  have  a  ramp  waveform  with  more  than  three  significant 
zero  crossings,  but  a  first  estimate  of  this  waveform  would  proba¬ 
bly  not  recognize  this  fact. 
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whereby  the  ramp  response  or  its  derivatives  is  constrained  to  a 
specific  value  at  a  fixed  time.  It  is  not  clear  at  this  point  which 
constraints  should  be  utilized,  aside  from  that  in  Eq.  (  58)  ,  or  at 
what  point  in  the  development  they  should  be  introduced.  Note  that 
the  type  of  constraint  in  Eq.  (  60)  is  an  obvious  way  to  introduce 
into  the  model  both  a  knowledge  of  the  arrival  time  of  diffracted 
fronts  and  asymptotic  estimates  of  their  value.  When  constraints 
on  the  orthonormal  expansion  are  introduced,  minimization  of  the 
rms  error  in  Eq.  (  43)  requires  the  use  of  a  Lagrange  multiplier. 

A  generalized  method  for  the  case  of  several  constraints  of  the 
form  of  Eq.  (  60 )  ,  yielding  the  Lagrange  multipliers  and  the  new 
expansion  coefficients  has  been  given[  29]  . 


A  systematic  procedure  for  progressing  beyond  the  ramp  re¬ 
sponse  waveform,  i.  e.  ,  to  the  step  and  impulse  response,  is  more 
difficult  to  formulate.  The  first  estimate  of  the  step  response  (  adding 

M  terms  to  the  ramp  model) 

N+M 

(61)  Fu(t')  An4>;<t') 

n=0 

has  M  undetermined  coefficients.  The  constraint  in  Eq.  (  58)  must 
now  be  imposed  on  the  N+M  coefficients.  It  need  not  have  been  im¬ 
posed  explicitly  on  the  N  coefficients  of  the  initial  ramp  waveform 
since  the  sketched  or  synthesized  ramp  response  satisfied  the  second 
moment.  The  form  of  solution  with  the  constraint  requires  that  the 
M  added  coefficients  be  calculated  without  the  constraint  and  then  a 
new  set  of  N+M  coefficients  be  computed.  Therefore  M  adjustments 
on  the  step  response  waveform  can  be  made.  That  is,  the  corrected 
step  waveform  is 

N+M 
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N+M 

A  *a<°> 

n=0 

\  - - 

N+M 

£  *n<°> 
n=0 

The  problem  is  in  calculating  the  M  new  coefficients  without  the  con¬ 
straint.  It  is  these  coefficients  which  are  to  add  the  necessary  ad¬ 
ditional  detail  to  the  step  waveform.  Given  sufficient  phasor  response 
data  to  estimate  the  rms  error  in  the  step  response  waveform  in  Eq. 

(  61)  ,  the  new  parameters  could  be  adjusted  to  reduce  this  error. 
Lacking  such  data,  two  alternatives  are  seen.  The  first  is  to  in¬ 
corporate  certain  theoretical  estimates  of  the  step  response  detail 
into  a  revised  sketch  of  the  waveform  and  then  proceed  as  with  the 
ramp  model  to  calculate  the  coefficients.  New  coefficients  would 
then  be  computed  as  in  Eq.  (  62).  The  second  alternative  is  to  in¬ 
corporate  additional  constraints  of  the  form  in  Eq.  (  60)  ,  calculate 
the  additional  unconstrained  coefficients  from  the  ramp  model,  and 
then  solve  for  the  new  coefficients  subject  to  the  constraints.  When 
more  than  one  constraint  is  introduced,  a  system  of  linear  equations 
must  be  solved  to  determine  the  Lagrangian  multipliers.  Essentially 
the  same  problems  and  alternatives  exist  in  going  to  the  impulse  re¬ 
sponse  waveform.  In  this  case  one  should  not  attempt  to  incorporate 
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into  the  model  delta  function  or  higher  order  singularities  in  the  wave¬ 
form,  since  these  can  be  added  independently. 

An  attempt  has  been  made,  in  the  preceding  pages,  to  outline 
a  systematic  general  procedure  for  estimating  the  canonical  response 
waveforms  in  a  transient  electromagnetic  scattering  problem.  The 
effective  integration  of  analytical  and  experimental  data  into  a  single 
model  is  stressed,  along  with  an  orderly,  predictable  method  for 
determining  the  parameters  of  the  model.  The  use  of  orthonormal 
expansions  of  the  response  waveform  estimates  has  been  suggested, 
and  certain  advantages  of  such  an  approach  discussed.  The  feasibility 
of  the  approach  has  not  been  demonstrated,  and  hinges  primarily  on 
the  availability  of  a  set  of  functions,  orthogonal  and  complete  on  the 
semi-infinite  interval  {  0,®)  ,  which  will  permit  a  reasonable  estimate 
of  the  ramp  waveform  with  just  a  few  terms.  Two  sets,  orthogonal 
exponentials  and  Laguerre  polynomials,  satisfying  the  orthogonality 
and  completeness  criteria  have  been  suggested,  but  their  applica¬ 
bility  has  not  been  shown.  No  claim  is  made  that  the  suggested  attack 
on  the  problem  represents  an  optimum  procedure,  and  it  is  evident 
that  future  research  should  investigate  not  only  orthonormal  expan¬ 
sions  but  other  possible  approaches. 
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One  rather  important  conclusion  can  be  drawn  from  the  discus¬ 
sion  in  this  section.  For  any  given  combination  of  theoretical  esti¬ 
mates  and  experimental  data,  the  best  hope  for  achieving  a  waveform 
estimate  which  closely  approximates  the  true  object  response  lies  with 
the  ramp  response  waveform.  Corresponding  estimates  of  the  step 
and  impulse  response  waveforms  will  be  less  exact,  if  derived  from 
limited  theoretical  and  experimental  data. 


i  • 


CHAPTER  II 

THE  ELECTROMAGNETIC  BACKSCATTERING 
BY  A  PERFECTLY  CONDUCTING 
PROLATE  SPHEROID 


The  investigation  of  scalar  and  vector  boundary  value  problems 
in  general  ellipsoidal  coordinates,  and  more  particularly  in  the  de¬ 
generate  cases  of  prolate  and  oblate  spheroidal  coordinate  systems, 
has  been  a  favorite  of  the  mathematical  analysist  for  almost  nine 
decades.  Flammer[  32]  in  his  "Spheroidal  Wave  Functions"  gives 
an  excellent  historical  survey  of  the  main  contributions  to  the  develop¬ 
ment  of  the  spheroidal  wave  functions  and  the  applications  of  the3e 
functions  to  various  boundary  value  problems.  More  recently, 

Sleator[  33]  has  reviewed  the  analytical  and  experimental  contributions 
to  the  acoustical  and  electromagnetic  scattering  by  prolate  and  oblate 
spheroids.  It  is  vividly  evident  in  this  latter  work  that  despite  the 
long  history  of  interest  in  the  subject,  very  few  numerical  results 
have  been  obtained  for  the  vector  problems.  The  reasons  for  this  are 
discussed  by  both  Flammer  and  Sleator,  and  are  briefly  reviewed 
below.  The  point,  however,  is  that  a  good  approximate  general  solu¬ 
tion  to  the  plane  wave  scattering  by  spheroids  has  not  been  obtained 
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even  in  a  formal  sense.  Those  approximate  solutions  which  are 
amenable  to  computation,  even  with  modern  high  speed  digital  com¬ 
puters,  are  limited  either  to  a  particular  range  of  the  spectrum  or 
to  specific  orientations  of  the  spheroid. 

Recently,  solutions  for  the  scattering  from  conducting  spheroids 
have  been  reported  by  Andreasenf  34]  ,  Oshiro[  3  5]  ,  and  by 
Waterman[  36] .  These  results,  while  differing  somewhat  in  specifics, 
are  essentially  all  applications  of  the  "point  matching"  approach  to 
electromagnetic  scattering  problems.  As  such,  there  is  a  practical 
limitation  imposed  by  the  size  of  present  electronic  computers  and  the 
costs  involved  on  the  maximum  spheroid  size  which  can  be  handled 
successfully.  Estimates  of  the  maximum  size  which  can  be  handled 
vary,  and  much  depends  upon  what  one  considers  reasonable  in  terms 
of  computer  running  time.  For  example,  Waterman[  36]  reports  a 
running  time  of  two  hours  on  an  IBM  7030  computer  to  compute  the 
scattering  pattern  of  a  2:1  prolate  spheroid  with  a  minor  circumference 
of  10  wavelengths.  In  any  event,  such  solutions  must  be  supplemented 
by  approximate  high  frequency  methods  to  cover  the  entire  spectrum. 

A  most  useful  application  of  these  "point  matching"  solutions  for 
sphe  roid  targeVs  would  be  the  computation  of  sufficient  phasor  response 
data  at  harmonically  related  low  frequencies  to  generate  the  ramp  re¬ 
sponse  waveform  for  specific  orientations  via  Fourier  synthesis. 
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In  this  Chapter,  an  approximate  general  solution  for  the  electro¬ 
magnetic  backscattering  by  a  perfectly  conducting  prolate  spheroid  is 
derived  using  time  domain  concepts.  The  formulas  obtained  hold  for 
arbitrary  orientation  and  axial  ratio  of  the  spheroid  and  for  arbitrary 
linear  polarization  of  the  incident  plane  electromagnetic  wave.  These 
results  have  the  extremely  attractive  feature  that  calculations  using 
the  derived  formulas  are  quite  simple,  involving  nothing  more  compli¬ 
cated  than  exponential  functions. 

The  major  difficulty  in  obtaining  a  rigorous  solution  of  the  plane 
wave  scattering  by  spheroids  is  that  the  vector  wave  equation  is  not 
separable  in  either  oblate  spheroidal  or  prolate  spheroidal  coordi¬ 
nates^].  Specifically,  what  goes  wrong  for  the  prolate  spheroid  can 
be  seen  from  a  consideration  of  the  simplest  case,  that  of  a  plane 
electromagnetic  wave  incident  along  the  major  axis  of  a  perfectly  con¬ 
ducting  prolate  spheroid.  This  is  the  problem  which  was  considered 
by  Schultz[  37]  .  As  shown  by  Schultzf  37]  ,  the  incident  field  can  be 
expanded  in  an  infinite  series  of  vector  spheroidal  wave  functions. 

An  expansion  for  the  scattered  field  as  the  sum  of  two  infinite  series 
of  vector  spheroidal  wave  functions  can  also  be  devised,  and  such  an 
expansion  properly  is  a  solution  of  the  vector  Helmholtz  equation, 
has  zero  divergence,  and  satisfies  the  radiation  condition  at  infinity. 
However,  in  attempting  to  satisfy  the  boundary  conditions  at  the 
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spheroid  surface  it  is  found  that  the  field  cannot  be  resolved  into  com¬ 
ponents  such  that  individual  terms  in  the  series  can  satisfy  the  bounda¬ 
ry  conditions.  This  difficulty  arises  for  two  reasons:  (  1)  the  angle 
functions  have  different  indices  in  the  expansion  of  the  incident  and 
scattered  fields,  and  (  2)  the  scale  factors  involved  in  the  definition 
of  the  vector  spheroidal  functions  appear  explicitly  in  the  expansions. 

In  the  classical  solution  of  the  sphere,  the  angular  dependence  is  the 
same  term  by  term  in  the  expansions  of  the  incident  and  scattered 
fields.  Therefore  one  can  utilize  the  orthogonality  of  the  angular 
functions  to  obtain  equations  relating  the  coefficients  of  the  expansions 
of  the  incident  and  scattered  fields.  For  the  spheroid,  this  cannot  be 
done,  i.  e.  ,  the  known  and  unknown  coefficients  cannot  be  related  for 
a  finite  number  of  summation  indices.  It  is  at  this  point  that  the  rigor 
of  Schultz's  solution  is  lost  since  an  infinite  system  of  equations  for 
the  infinite  set  of  coefficients  must  be  truncated  to  obtain  a  solution 
for  the  expansion  coefficients  of  the  scattered  field. 

Although,  as  related  by  Sleator[  33]  ,  certain  advances  have 
been  made  in  the  approximate  or  asymptotic  theories  for  high  and  low 
frequencies,  the  returns  from  the  numerous  studies  of  the  problem 
have  been  woefully  meager.  The  situation  is  summarized  (prior  to 
the  "point  matching"  results)  by  Sleator[  33]  as  follows: 


"Since  the  work  of  Schultz[  37]  and  the  computations  based 
on  this  by  Siegel,  et  al.  [  38]  ,  virtually  no  progress  has 
been  made  in  the  solution  of  the  vector  problem  in  the  reso¬ 
nance  region.  All  existing  techniques  either  break  down 
completely  or  become  prohibitively  difficult  or  tedious  in 
this  region,  and  the  need  for  a  totally  new  approach  becomes 
more  and  more  apparent.  " 

It  is  to  this  new  approach  that  the  work  in  this  Chapter  is  directed. 
The  new  approach  is  an  application  of  the  time  domain  concepts  dis¬ 
cussed  in  Chapter  I.  The  application  of  these  concepts  is  in  its 
infancy,  i.  e.  ,  a  new  approach  to  electromagnetic  scattering  problems 
is  being  developed,  and  the  approach  is  via  a  "model"  rather  than  a 
sophisticated  extension  and  application  of  a  well-known  technique. 
Consequently  the  treatment  will  be  mathematically  simple  and  the 
philosophy  of  an  engineering  approximation  liberally  applied.  The 
goal  of  this  study  is  the  development  of  a  general  model  for  the 
electromagnetic  backscattering  by  a  perfectly  conducting  prolate 
spheroid  for  arbitrary  orientation  and  axial  ratio  and  arbitrary  linear 
polarization  of  the  incident  field.  What  is  sought  is  an  estimate  of  the 
response  waveform  of  the  target  capable  of  yielding  a  reasonable 
approximation  of  the  scattered  field  for  a  variety  of  orientations, 
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polarizations,  and  frequencies.  A  precise  determination  of  the  model 
for  one  fixed  combination  of  these  parameters  is  not  intended. 

The  first  section  of  this  Chapter  briefly  reproduces  a  solution 
for  the  axial  impulse  response  waveform  of  the  spheroid  published  by 
Moffatt  and  Kennaugh[  1 1]  in  1965.  The  reason  for  the  inclusion  of 
this  result  is  two-fold;  first,  much  of  the  study  for  the  general  case  is 
predicated  on  retaining  the  essential  character  of  the  axial  waveform, 
and  second,  the  axial  result  is  intimately  related  to  the  response  wave¬ 
form  of  a  sphere,  an  understanding  of  which  is  essential  for  later  dis¬ 
cussions  of  the  spheroid  waveforms. 

With  this  background  a  basic  general  model  for  the  impulse  re¬ 
sponse  waveform  of  the  spheroid  is  then  derived  and  calculated  results 
from  this  model  compared  with  measured  data  on  a  2:1  axial  ratio 
spheroid  and  calculated  data  (axial)  for  a  10:1  spheroid.  The  Chapter 
concludes  with  a  discussion  of  these  results  and  the  applicability  of 
the  basic  model  for  spheroids  with  other  axial  ratios. 

In  Appendix  II,  that  portion  of  the  sphere  impulse  response 
waveform  corresponding  to  the  creeping  wave  is  examined  in  some 
detail.  This  digression  to  the  sphere  is  used  to  illustrate  the  basic 
nature  of  the  creeping  wave  contribution  in  the  time-domain.  The 
creeping  wave  contributions  predicted  by  high-frequency  asymptotic 
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solutions  are  analyzed  and  certain  inconsistencies  resulting  from  their 
non-causal  nature  discussed. 

A.  The  Axial  Impulse  Response 
Waveform 

Consider  the  perfectly  conducting  sphere  and  the  related  per¬ 
fectly  conducting  prolate  spheroid  with  an  n  to  1  axial  ratio  shown  in 
Fig.  6.  The  spheroid  has  a  semiminor  axis  equal  to  the  sphere  radius 


Fig.  6--Sphere  and  prolate  spheroid  geometries. 
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which  is  taken  as  unity.  Let  the  sphere  and  spheroid  be  illuminated 
by  a  ramp  type  excitation  propagating  in  the  z  and  z  directions  re¬ 
spectively. 

(  63)  E(t')  =  t'  u(t’) 

E(F)  =  t1  u(T') 

where  t '  =  t  -  z/c  and  t'  =  t  -  z/c 

Now  if  the  axial  ratio  of  the  spheroid  is  unity,  clearly  the  response 
waveforms  will  be  identical  for  the  sphere  and  spheroid.  This,  of 
course,  is  true  for  any  excitatioa  As  the  axial  ratio  of  the  spheroid 
is  increased  (n  >  1)  differences  in  the  sphere  and  spheroid  response 
waveforms  will  occur.  The  basis  for  the  approximate  spheroid  re¬ 
sponse  waveform  derived  in  this  section  is  the  assumption  that  the 
ramp  response  waveform  of  the  spheroid  can  be  obtained  by  simply  a 
change  in  the  time  scale  of  the  sphere  ramp  response  waveform. 
Since  the  step  and  impulse  response  waveforms  may  be  obtained  by 
successive  differentiations  of  the  ramp  response  waveform,  a 
change  in  time  scale  for  the  ramp  response  waveform  results  in  a 
change  in  amplitude  as  well  as  time  scale  for  the  step  and  impulse 
response  waveforms.  To  a  first  approximation  at  least,  the  ramp 
response  waveform  of  the  spheroid  undergoes  only  a  time  scale 
change  as  can  be  shown  in  two  ways. 
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For  the  sphere,  using  Eq.  (  24)  of  Chapter  I,  the  response  wave¬ 
forms  predicted  by  the  physical  optics  approximation  are 

(64)  Fr(  f)  =  T-  +  llj  [U(t')  -u(t'-t0>] 

Fuft’)  =  ["  7+  7]  [“(I')  -“( t'-t0)  ] 

Fjtt')  «  -  k  8(F)  +  -  [  u(  t')  -u(  t'-t0)  ] 

2  2 

and  those  for  the  spheroid  are 
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where  tQ  =  2a/c,  and  a  bar  has  been  used  to  denote  the  spheroid  wave¬ 
forms.  If,  in  the  expression  for  the  ramp  response  waveform  of  the 
sphere  the  substitution  t'  =  t'/n  is  made  then  the  ramp  response  wave¬ 
form  of  the  spheroid  is  obtained.  Obviously  successive  differentiation 
leads  to  the  correct  physical  optics  step  and  impulse  responses.  Note 
that  the  same  change  in  time  scale  in  the  step  or  impulse  response  of 
the  sphere  does  not  yield  the  spheroid  physical  optics  result.  It  is 
concluded  therefore  that  insofar  as  the  physical  optics  approximation 
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is  concerned  it  is  the  ramp  response  waveform  which  undergoes  only  a 
time  scale  change. 

A  second  argument  for  only  a  time  scale  change  in  the  ramp 
waveform  is  provided  by  the  moment  conditions;  specifically  the  second 
moment  which  relates  the  area  beneath  the  ramp  response  waveform 
and  the  Rayleigh  coefficient  of  the  target.  The  approximate  axial  ramp 
response  waveform  of  the  spheroid  according  to  this  assumption  is 
obtained  by  a  linear  deformation  of  the  time  axis  in  the  illuminated 
region  and  a  nonlinear  deformation  in  the  shadowed  region.  Corre¬ 
sponding  ordinates  on  the  ramp  response  waveforms  for  the  sphere 
and  spheroid  are  defined  by  equal  cross  sections  in  a  plane  perpen¬ 
dicular  to  the  line  of  sight.  Corresponding  abscissa  are  defined  by 
the  shortest  ray  path  to  the  given  cross  section  following  the  line  of 
sight  to  the  shadow  boundary  and  the  scatterer  surface  beyond. 

Thus  setting 


(66)  Fr(  t ' )  =  Fr(  t')  , 


and  differentiating 


(67) 


F0(t') 


FU<*’> 


dt1 

dF 


=  Fj(t') 


Fj(t') 


+  Fu(t') 
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(68) 


(69) 

and 


(70) 


where 


0  <LV  <£t0 

t0  <  t'  <  tQ  +  S  , 


An  example  of  the  spheroid  ramp  response  waveform  obtained  in  this 
way  is  shown  in  Fig.  7b  for  a  prolate  spheroid  with  a  2:1  axial  ratio. 
The  sphe  re  ramp  response  waveform  is  shown  in  Fig.  7a.  A  graphical 
integration  of  the  spheroid  ramp  response  waveform  shown  in  Fig.  7b 


Fig.  7a- -Normalized  ramp  response  (axial)  and  associated  geometry.  Sphere. 
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Fig.  7b- -Normalized  ramp  response  (axial)  and  associated  geometry.  Spheroid. 
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demonstrates  that  the  second  moment  condition  is  approximately  satis¬ 
fied.  A  comparison  of  the  Rayleigh  coefficient  obtained  by  graphical 
integration  of  the  curve  in  Fig.  7b  and  that  calculated  in  Appendix  I 
shows  an  error  of  approximately  3  percent.  Similar  comparisons  for 
axial  ratios  of  4,6,8,  and  10  gave  a  minimum  error  of  1  percent  (8:1) 
and  a  maximum  error  of  5  percent  (4:1).  No  discernable  trend  of 
error  versus  axial  ratio  was  observed.  The  ramp  response  waveform 
of  the  sphere  (Fig.  7a)  was  obtained  from  the  exact  Mie  solution  using 
the  Fourier  synthesis  procedure.  The  synthesis  procedure  is  described 
in  Reference  13.  Thus  both  physical  optics  and  the  second  moment  con¬ 
dition  (  as  do  the  two  lower  orders)  support  the  assumption  that  the 
ramp  response  waveform  of  the  spheroid  is  obtained  from  that  of  the 
sphere  by  a  change  in  the  time  scale  alone. 

It  is  important  to  emphasize  one  point  in  the  above  results.  The 
Fourier  synthesis  procedure  used  to  obtain  the  approximate  impulse, 
step  and  ramp  response  waveforms  of  the  sphere  does  not  yield  ana¬ 
lytical  expressions  for  these  waveforms,  but  only  their  magnitudes  at 
selected  values  of  time.  The  time  increments  are  arbitrary  so  that 
the  waveforms  may  be  obtained  to  any  accuracy  desired  (within  the 
limitation  of  a  finite  pulse  width)  ,  but  only  in  the  sense  of  a  point  by 
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point  plot.  Consequently  the  ramp,  step  and  impulse  response  wave¬ 
forms  of  the  spheroid  obtained  by  the  above  described  method  are 
also  only  point  by  point  plots  at  arbitrary  values  of  time. 

A  detailed  development  of  the  axial  response  waveforms  of  the 
spheroid  is  given  in  the  1965  publication  by  Moffatt  and  Kennaugh[  11]  . 
Shown  in  Fig.  8  (  solid  curve)  is  the  impulse  response  waveform  of 
a  2;1  axial  ratio  spheroid  as  obtained  by  a  nonlinear  transformation 
of  the  time  scale  for  the  ramp  response  waveform  of  the  sphere.  Also 
shown  in  Fig.  8  (dashed  curve)  is  an  analytical  exponential  sum 
approximation  for  the  waveform.  The  Laplace  transform  of  this 
exponential  sum  approximation  was  found  to  yield  a  satisfactory  fre¬ 
quency-domain  approximation  for  minor  axis  circumferences  of  the 
spheroid  in  wavelengths  greater  than  1.  5.  For  minor  axis  circumfer¬ 
ences  in  wavelengths  less  than  1.  5  it  was  necessary  to  use  convolution 
and  an  approximate  integration  of  the  derived  ramp  response  wave¬ 
form  to  obtain  frequency  domain  data.  A  comparison  of  calculated 
and  measured  echo  area  for  a  2:1  axial  ratio  spheroid  is  shown  in 
Fig.  9.  Data  denoted  as  target  motion  are  explained  in  Reference  11. 

It  has  been  shown  that  an  approximate  solution  for  the  axial 
backscattering  from  a  prolate  spheroid  may  be  obtained  through  a 
modification  of  the  impulse  response  waveform  of  a  spherical  scatterer. 
Experimental  verification  of  the  results  predicted  was  obtained  for  a 
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Fig.  8- -Approximate  impulse  response  waveform  (axial)  of  2:1  axial  ratio 
prolate  spheroid.  Dashed  curve  is  exponential  sum  approximation. 
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prolate  spheroid  of  2:1  axial  ratio.  The  method  used  employs  a  non¬ 
linear  transformation  of  the  time  scale  to  obtain  the  ramp  function  re¬ 
sponse  of  the  spheroid  from  that  of  the  sphere.  The  method  has  not 
been  applied  to  extremely  large  axial  ratios  and  may  require  further 
corrections  in  such  cases. 

The  model  has  a  number  of  deficiencies  which  need  to  be  studied 
further.  First,  for  times  after  the  impulsive  excitation  reaches  the 
rear  of  the  spheroid  it  was  assumed  that  the  time  relations  (  Eqs.  (  69) 
and  (  70)  )  were  the  same  as  the  first  portion,  corresponding  to  the 
excitation  progressing  on  around  the  spheroid.  This  assumption,  used 
primarily  for  lack  of  an  alternative  procedure,  needs  to  be  placed  on  a 
sound  basis  or  a  different  technique  devised.  Second,  it  can  be  seen 
from  Eqs.  (  67)  ,  (  69)  ,  and  (  70)  ,  that  the  creeping  wave  peak  predicted 
is  independent  of  the  spheroid  axial  ratio  and  equals,  in  all  cases,  that 
of  the  sphere;  a  questionable  result.  Both  of  these  questions  are 
related  to  the  nature  of  the  creeping  wave  on  smooth  convex  surfaces. 

In  Appendix  II,  some  insight  to  this  problem  is  gained  from  the  time 
domain  interpretation  of  asymptotic  estimates  of  the  creeping  wave 
on  spheres  and  circular  cylinders.  However,  in  this  study  the  results 
in  Appendix  II  have  not  been  utilized  to  improve  the  axial  model.  In 
terms  of  the  basic  goal  of  this  Chapter,  that  of  a  general  model,  two 
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results  of  this  section  are  used.  One,  the  fact  that  an  exponential 
approximation  of  the  creeping  wave  consisting  of  a  join  of  two  expo¬ 
nentials  with  finite  slopes  at  the  time  of  the  creeping  wave  peak 
yielded  reasonable  results,  and  two,  the  fact  that  the  magnitude  of 
this  creeping  wave  peak  could,  at  least  for  a  2:1  axial  ratio  spheroid, 
be  assumed  to  be  independent  of  axial  ratio. 

B.  A  General  Impulse  Response 
Waveform  Model 

For  the  general  problem,  the  prolate  spheroid  and  the  associated 
coordinates  are  shown  in  Fig.  10.  The  spheroid  has  a  semiminor  axis 
a  and  a  semimajor  axis  b.  The  incident  plane  electromagnetic  field  is 
incident  from  a  direction  defined  by  the  spherical  angles  0  and 
Because  of  the  rotational  symmetry  of  the  spheroid,  the  direction  of 
propagation  of  the  incident  field  can  be  restricted  to  lie  in  the  y,z 
plane  without  loss  of  generality.  With  this  restriction,  two  principal 
polarizations  of  the  incident  field  are  defined;  T.  E.  ,  where  the 
incident  electric  field  vector  is  normal  to  the  y,z  plane,  and  T.  M.  , 
where  the  incident  electric  field  vector  lies  in  the  y,  z  plane. 

Associated  with  each  of  the  two  principal  polarizations  is  a  path  length 
corresponding  to  line  of  sight  to  the  shadow  boundary  and  a  geodesic 
on  the  surface  beyond.  Note  that  these  paths  are  measured  from  a 
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Fig.  10--The  prolate  spheroid  and  associated  coordinates. 

plane  perpendicular  to  the  line  of  sight,  and  located  at  the  point  where 
such  a  plane  first  touches  the  spheroid. 

In  an  earlier  study  of  the  general  prolate  spheroid  target[  12]  , 
the  author  derived  an  approximate  spheroid  model  based  on  an  approxi 
mation  to  the  ramp  response  waveform  of  the  target.  The  results  of 
this  initial  study  are  of  some  interest;  both  as  an  example  of  how  a 
very  crude  approximation  in  the  time  domain  can  yield  reasonable 
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results  in  the  frequency  domain,  and  as  a  demonstration  of  the  neces¬ 
sity  for  a  well  choosen  impulse  response  model  in  order  to  achieve 
fine  accuracy  in  the  frequency  domain.  A  general  discussion  of  the 
derivation  of  transient  and  impulse  response  waveforms  for  a  finite, 
three  dimensional  target  was  given  in  Chapter  I;  it  is  reiterated  how¬ 
ever  that  it  is  reasonable  that  the  twice  integrated  function  of  a  real 
variable  be  easier  to  approximate  with  fewer  parameters  than  the 
function  itself. 

The  above  model  for  the  prolate  spheroid  target  was  based  upon 
an  earlier  study  by  Kennaugh[  39]  demonstrating  that  the  ramp  re¬ 
sponse  waveform  of  a  conducting  sphere  could  be  approximated  quite 
well  using  very  simple  functions.  In  this  approach,  the  ramp  response 
waveform  consists  of  the  physical  optics  approximation  to  the  shadow 
boundary  and  an  amplitude  proportional  to  the  circumference  of  the 
section  cut  from  the  sphere  surface  by  a  plane  moving  along  the  line 
of  sight  with  half  the  free  space  velocity  beyond  the  shadow  boundary. 
For  those  times  after  the  plane  moves  beyond  the  sphere,  this  latter 
waveform  is  extended  continuously  and  multiplied  by  an  exponentially 
damped  term  whose  magnitude  is  chosen  to  satisfy  the  second  moment 
condition.  The  amplitudes  of  the  above  two  sections  are  matched  at 
the  shadow  boundary.  The  resulting  approximation  to  the  ramp  re¬ 
sponse  waveform  of  a  conducting  sphere  of  radius  a  is 
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(71)  FR(t')  =  £-1211+^-1  [u(t')  -u(t'-t0)] 

-  cos(  t'-t0)  [  u(  t'-t0)  -u(t'-(  l  +  ir/2)t0)] 

4 

-  (to(  1  +  w/2)“t')  co8(t‘-t0)  u(  t'-(  l+w/2)t0) 

4 

Several  salient  features  of  this  waveform  should  be  noted.  The  first 
term  in  Eq.  (71)  is  simply  the  physical  optics  approximation  to  the 
ramp  response  waveform.  It  was  shown  in  Chapter  I  that  this  is 
simply  proportional  to  the  cross-sectional  area  intercepted  by  a 
transverse  plane  moving  in  the  direction  of  the  incident  wave  at  half 
the  free  space  velocity.  The  second  term  corresponds  to  the  cir¬ 
cumference  of  this  cross-sectional  area  beyond  the  shadow.  Note, 
however,  that  the  times  corresponding  to  these  amplitudes  are  meas¬ 
ured  from  the  shadow  along  the  surface  of  the  sphere.  Thus  the  cir¬ 
cumference  reaches  zero  at  t'  =  (  1  +  it/2)  t0  rather  than  2t0.  As  was 
shown  by  Kennaugh[  39]  ,  this  waveform  is  a  rather  good  approximation 
to  the  ramp  response  waveform  obtained  by  Fourier  synthesis  from 
the  exact  Mie  solution  and  shown  in  Fig.  7a.  As  will  be  seen,  how¬ 
ever,  such  a  waveform  leads  to  certain  difficulties  where  twice 
differentiated  to  obtain  an  impulse  response  waveform. 
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One  can  derive  a  similar  approximation  for  the  ramp  response 
waveform  of  a  prolate  spheroid  for  the  general  off-axis  case.  Note 
that  in  this  case  a  polarization- sensitive  approximation  is  obtained 
because  of  the  dependence  of  the  second  moment  condition  on  the 
Rayleigh  coefficient  of  the  target,  which  is  dependent  on  both  the 
incident  polarization  and  the  orientation  of  the  spheroid.  In  order 
to  secure  a  closed  form  expression  for  the  ramp  response  waveform, 
certain  simplifying  approximations  are  made  for  the  spheroid  case. 
The  first  of  these  involves  the  point  of  termination,  in  time,  of  the 
physical  optics  approximation.  For  the  sphere,  the  shadow  boundary 
is  always  normal  to  the  direction  of  propagation  of  the  incident  field. 
For  the  spheroid,  however,  this  only  occurs  at  axial  and  broadside 
(0=0  and  90  degrees  respectively  in  Fig.  10).  At  other  orientations, 
the  shadow  boundary  is  not  perpendicular  to  the  line  of  sight.  As  a 
crude  approximation,  the  shadow  boundary  was  ignored  and  the 
physical  optics  approximation  terminated  at  the  peak  transverse 
cross  section  encountered  by  the  cutting  plane.  A  second  difficulty 
encountered  is  in  the  measurement  of  path  lengths  on  the  shadowed 
side  of  the  spheroid.  This  is,  for  either  polarization,  a  distance 


*  The  Rayleigh  dependence  of  prolate  and  oblate  spheroids  for 
arbitrary  polarization  and  orientation  and  for  arbitrary  constitutive 
parameters  is  given  in  Appendix  I. 
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along  the  perimeter  of  an  ellipse,  i.  e.  ,  an  incomplete  elliptic  integral. 
The  circumlerence  of  the  area  cut  from  the  surface  by  the  cutting  plane 
is  also  an  ellipse.  Again,  in  the  interests  of  a  simple  closed  form 
expression,  the  second  portion  of  the  waveform  was  crudely  approxi¬ 
mated  by  a  shifted  cosine  wave  whose  period  was  determined  by  the 
creeping  wave  geodesic  path,  i.  e.  ,  T.  E.  or  T.  M. .  The  ramp  re¬ 


sponse  waveform  so  obtained  is  given  by 

(72)  FR(t')  =  [ -At' +  Bt'2/2]  fu(t')  -u(t’  -  1 10)  1 

-  Tb  c08[P(t,-|t°)][u(t,-tto)  -U(t'-T2t°>] 

'T5e°(Tjto’t'>cos[p(t' -f*®)  u<t'-T*to> 
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The  constants  A  and  B  are  obtained  from  the  physical  optics  ap¬ 
proximation  to  the  ramp  response  waveform.  The  constant  a  is  de¬ 
termined  from  the  second  moment  condition  and  K  is  the  Rayleigh  coef¬ 
ficient  of  the  sphe  roid  (Appendix  I).  The  time  T2  is  determined  from 
the  geometry  of  the  spheroid  and  corresponds  to  a  creeping  wave  path 
length,  following  the  line  of  sight  to  the  shadow  boundary  and  a 
geodesic  path  on  the  spheroid  beyond.  Assuming  that  the  circumfer¬ 
ence  of  the  ellipse  in  the  x,  z  plane  may  be  expressed  as 

(73)  Cir=27r 

with  negligible  error,  T2  is  given  by 


T.  M.  Polarization 


and 


T.  E.  Polarization 


(75)  _  sin2  6  cos2  0  (b2  -  a2  )2 


+  Jb2cos20+  a2sin20 


Jb2  cos2  0+a2  sin2  0  (b2  sin2  0  +  a2  cos2  0) 


+  2 


it  ci2  +  c22 
2  .  2 
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Cj2  =  a2  1 


sin2  0  cos2  0(b2  -  a2)(a2-  b2) 

sin20  co820(b2  -  aV  +  a2b2,  b2cos20+  a2sin20 
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b  2sin20  +  a2  cos2  0 
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b  cos  0+a  sin20  (b  sin  0+a  cos  0) 
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In  actual  calculations  in  this  study,  the  quarter  circumference  of  an 
ellipse  which  appears  in  the  T.  M.  case  as 

ir  r^Tb*— 

2^2 

and  in  the  T.  E.  case  as 


tr  C?  +  Cl 


was  computed  using  the  polynomial  approximation  to  the  complete 


elliptic  integral[  40] 


(76) 


Cir  =  4b  ]  1+a,  -  +a2  - 


a\  4 


where  aj  =  0.  4630151 
a2  =  0.  1077812 
bj  =  0.  2452727 
b2  =  0.  0412496 


for  an  ellipse  of  major  axis  2b  and  minor  axis  2a.  Note  that  for  a 
unity  axial  ratio  spheroid,  Eq.  (72)  reduces  to  Eq.  (71). 
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The  ramp  response  waveform  in  Eq.  (  72)  was  used  in  Reference 
12  by  the  author  as  the  basic  model  for  the  conducting  prolate  spheroid. 
Calculations  of  the  echo  area  of  the  spheroid  using  this  model  were 
compared  with  measured  data  on  a  2:1  axial  ratio  spheroid.  The  com¬ 
parison  covered  a  frequency  range  from  1.  5  to  5.  0  minor  axis  cir¬ 
cumferences  in  wavelengths  in  increments  of  0.  25  wavelenths  and  an 
aspect  range  of  0(Fig.  10)  from  0°  to  90°  in  increments  of  10°,  for 
both  the  T.  E.  and  T.  M.  polarizations.  Considering  the  crudeness  of 
the  ramp  waveform  model)  the  results  of  this  comparison  were  quite 
good.  An  example  of  the  type  of  agreement  obtained  is  shown  for  an 
aspect  angle  0  of  40°  in  Figs.  11a  and  lib  for  the  T.  E.  and  T.  M. 
polarizations  respectively.  Note  that  the  basic  deficiency  for  both 
polarizations  appears  to  be  a  constant  shift  of  the  maxima  and  minima 
of  the  calculated  data  with  respect  to  the  measured  data.  But,  taking 
into  account  this  shift,  the  calculated  and  measured  periodicity  and 
amplitudes  are  in  fair  agreement.  It  would  serve  no  purpose  here 
to  duplicate  all  of  the  results  obtained  with  this  ramp  response  model. 
These  are  reported  in  detail  in  Reference  12  and  the  reader  is 


referred  to  that  publication. 
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Fig.  1  la--Calculated  (  solid)  and  measured  echo  area 
of  a  2jl  axial  ratio  prolate  spheroid.  Ramp  response 
model,  T.  E.  polarization,  9  =  40°  . 


Aside  from  the  specific  calculations  for  the  prolate  spheroid 
targets,  the  results  obtained  with  the  crude  ramp  response  model 
yielded  two  general  conclusions.  First,  the  study  reported  in  Refer¬ 
ence  12  was  the  initial  attempt  to  extend  the  application  of  transient 
and  impulse  response  approximations  to  a  nonsymmetrical  or  off-axis 
target  orientation.  As  such  they  were  distinctly  encouraging  and  lent 
real  hope  that  other  targets  might  be  amenable  to  a  similar  attack. 
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Fig.  1  lb- -Calculated  (  solid)  and  measured  echo  area 
of  a  2:1  axial  ratio  prolate  spheroid.  Ramp  response 
model,  T.  M.  polarization,  0  =  40°  . 


Second,  it  had  been  contended  by  E.  M.  Kennaugh  and  the  author  in 
1 9 65[  13]  that,  with  experience  and  intuition,  one  might  ultimately 
be  able  to  roughly  sketch  the  ramp  response  waveform  of  a  target 
from  simply  a  consideration  of  physical  optics,  the  target  geometry 
and  the  low-frequency  derived  moment  conditions.  Noting  that  to 
this  point,  the  above  considerations  were  the  only  ones  involved  in 
deriving  the  ramp  response  waveform,  it  is  felt  that  this  premise 
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has  Leen  proved  -  and  proved  for  a  target  which  heretofore  has  stub¬ 
bornly  resided  even  an  approximate  practical  theoretical  solution. 

Two  basic  deficiencies  can  be  noted  with  the  ramp  response 
waveform  given  in  Eq.  (  72).  These  are  most  apparent  from  a  con¬ 
sideration  of  the  related  impulse  response  waveform.  Differentiating 
Eq.  (  72)  twice  with  respect  to  time, 
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A  rough  sketch  of  the  form  of  the  impulse  response  waveform  in  Eq. 

(  77)  is  shown  in  Fig.  12.  From  a  knowledge  of  the  axial  impulse  re¬ 
sponse  waveform  in  Section  A  and  the  nature  of  the  creeping  wave  in 
Appendix  II,  it  is  clear  that  the  impulse  response  waveform  in  Eq.  (  77) 
has  two  primary  faults.  First,  the  erroneous  jump  occuring  at 
t'  =  A/b  tQ,  and  second  the  form  of  the  response  waveform  at  the 
creeping  wave  peak  (  T2tQ) .  In  the  first  case,  there  is  no  evidence 
that  in  the  neighborhood  of  the  shadow  boundary  a  discontinuity  should 
occur  in  the  waveform.  The  axial  model  (Section  A)  predicts  a  smooth 


98 


Fig.  12--Rough  sketch  of  impulse  response  waveform 
predicted  by  basic  ramp  model. 


continuous  waveform  in  this  region.  The  creeping  wave  portion  of  the 
waveform  in  Fig.  12  resembles,  if  anything,  that  to  be  expected  from 
a  cylinder  (  see  Appendix  II) .  In  the  earlier  study[  12]  ,  part  of  the 
discrepancy  in  the  shifted  periodicity  of  the  calculated  data  with  re¬ 
spect  to  the  measured  data  was  attributed  to  the  approximation  in 
Eq.  (  73)  with  regard  to  the  creeping  wave  times.  For  a  2jl  axial 
ratio  spheroid,  Eq.  (  72)  is  approximately  6  percent  in  error.  How¬ 
ever,  a  later  correction  of  this  error  using  Eq.  (  76)  did  not  signifi¬ 
cantly  alter  the  positions  of  maxima  and  minima  in  the  theoretical 
curves.  It  is  now  clear  that  the  character  of  the  creeping  wave 
contribution  in  conjunction  with  the  erroneous  jump  at  the  peak 
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cross-sectional  area  produced  the  error.  This  will  be  demonstrated 
with  calculations  on  the  improved  impulse  derived  model.  One  should 
note  that,  in  principle,  it  would  be  possible  to  correct  the  ramp  re¬ 
sponse  derived  model.  That  is,  a  polynomial  section  rather  than  a 
shifted  cosine  could  be  used  to  eliminate  the  jump  at  the  shadow 
boundary  and  at  the  same  time  .orce  the  waveform  to  assume  the 
desired  behavior  in  the  neighborhood  of  the  creeping  wave  maximum. 
However,  the  additional  parameters  required  would  exceed  the  known 
conditions  on  the  waveform.  This  is  the  type  of  dilemma  one  en¬ 
counters  when  starting  with  an  analytical  approximation  of  the  ramp 
response  waveform.  Since  this  waveform  is  smoother  than  its  first 
or  second  derivatives,  it  can  be  estimated  with  fewer  parameters. 

But  the  type  of  functional  dependence  assumed  is  extremely  critical, 
when  twice  differentiated,  in  the  impulse  response  waveform.  It 
will  be  seen,  for  example,  that  th  impulse  response  waveforms  of 
the  new  model  differ  markedly  from  that  sketched  in  Fig.  12. 

Turning  directly  to  the  derivation  of  a  model  which  will  have 
the  desired  features  lacking  above,  it  is  advantageous  to  start  with 
an  assumed  approximation  for  the  impulse  response  waveform  of  the 
target,  while  at  the  same  time  retaining  those  features  seen  to  be 
desirable  in  the  ramp  model.  There  seems  no  reason  to  doubt  the 
validity  of  the  physical  optics  approximation  for  short  times,  thus  the 
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specular  delta  function  and  jump  discontinuity  will  be  given  by  the  con¬ 
stants  A  and  B  respectively  in  Eq.  (  72).  It  is  desired,  with  this 
model,  to  secure  the  secondary  peak  type  behavior  for  the  creeping 
wave  maximum  occuring  at  the  time  T*  given  in  Eqs.  (  74)  and  (  75). 
Evidence  from  the  axial  waveform  in  Section  A,  the  asymptotic  esti¬ 
mates  in  Appendix  II,  and,  indirectly,  the  ramp  response  derived 
model  discussed  in  this  section,  indicates  this  type  of  waveform  for 
the  creeping  wave.  The  basic  impulse  response  waveform  for  the 
spheroid  is  assumed  to  be  of  the  form 

(78)  Fm(t')  =  -A5(t')  +  [B-Ce"*']  [u<  f)  -  u(t'-T,t0)] 

+  De'ttt'  u(  t'  -  T2tQ) 

where  the  subscript  H  denotes  this  to  be  the  high  frequency  portion  of 
the  final  impulse  response  waveform.  Two  conditions  are  imposed  on 
this  waveform 

(79)  C  <  <  <  B  , 
and 

(80)  B  -  Ce°Tzto  =  De'aTzto  =  P  , 

where  P  is  a  constant.  Note  that  this  high  frequency  portion  of  the 
waveform  has  deliberately  been  made  as  simple  as  possible  by  chosing 
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the  same  exponential  coefficient  on  either  side  of  the  creeping  wave 
peak.  The  waveform  is  designated  as  high  frequency  because  it  is  the 
relatively  high  frequency  properties  of  the  scatterer  which  it  must 
satisfy,  the  low  frequency  conditions  will  be  added  later.  From  Eq. 
(80), 

(81)  D  =  -  [  PB-P2]  , 

C 

and 


It  is  clear  from  Eqs.  (  81)  and  (  82)  that  two  constants  in  the  waveform 
of  Eq.  (78)  ,  P  and  C,  have  not  as  yet  been  specified.  With  these  two 
parameters,  it  is  possible  to  effect  a  degree  of  control  on  the  character 
of  the  creeping  wave  contribution.  For  example,  with  the  peak  creeping 
wave  contribution,  P,  fixed,  C  controls  the  slope  of  the  waveform  prior 
to  the  peak.  The  basic  assumption  with  this  model  is  that  the  impulse 
response  waveform  of  the  spheroid,  regardless  of  orientation  or  polar¬ 
ization,  should  be  roughly  similar  in  character  to  the  waveform  of  the 
conducting  sphere  and  the  waveform  derived  in  Section  A  for  the 
spheroid  axial  case.  Note  particularly  that  no  attempt  has  been  made 
to  utilize  the  low-frequency  derived  moment  conditions  to  determine 
the  constants  P  and  C  or  to  include  additional  parameters  and  enforce 
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higher  moment  conditions.  Such  an  approach  has  been  found  to  be  inef¬ 
fective,  primarily  because  one  cannot  conveniently  maintain  simul¬ 
taneously  the  desired  character  of  the  impulse  waveform.  Even  for 
very  simple  waveforms,  the  integrations  indicated  by  the  moment  con¬ 
ditions  can  lead  to  a  very  complicated  set  of  simultaneous  nonlinear 
equations. 

Two  possibilities  exist  for  establishing,  approximately,  the 
peak,  P,  of  the  creeping  wave  contribution.  The  first  is  an  evaluation 
of  the  asymptotic  expression  obtained  in  Appendix  II.  The  second  is  to 
accept  the  assumption  in  Section  A  that  the  peak  magnitude  is  invariant 
with  axial  ratio,  and  use  the  known  magnitude  of  the  sphere  waveform. 
In  either  case,  the  derivation  of  the  remaining  low-frequency  portion 
of  the  impulse  response  waveform  for  the  spheroid  is  unchanged.  It 
is  proposed  to  superimpose  a  second  waveform  on  that  given  in  Eq. 

(78)  whose  sole  purpose  is  to  ensure  the  correct  low-frequency  be¬ 
havior  of  the  response.  This  waveform  is  an  exponentially  damped 
sinusoid  which  is  written  in  the  form 

(  83)  FjJt')  =  -  e"*2*'}  u(  t') 

where  E,<j>i  ,  and  4>2  can  be  complex.  Let  the  first  three  moments  of 
the  high-frequency  waveform  in  Eq.  (  78)  be  denoted  by  Ij  ,  I2,  and  I3 


respectively,  i.  e.  , 
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(84)  ]  FjHft’Jdt'  =1,  , 

0 

etc.  Setting  the  impulse  response  waveform  of  the  spheroid  to  be 

(85)  Fi(t')  =Fm(f)  +FIL(f) 

and  imposing  the  zero*  first,  and  second  moment  conditions  on  the 
waveform  yields  the  set  of  simultaneous,  nonlinear  equations 


(86)  E[  ]  — Ii4>i  <t>2 

E[  44-4,5]  »  -I24>i4>1 
2E[  4>!-4>i]  =  (  K-I  )4>i4>l 


Solving  Eq.  (  86)  ,  if 


(87) 


9 


then  <j>i ,  <J>2,  and  E  are  real  and  given  by 


104 


(89) 


and 


9 


(90) 


E  =  - 


ll  <t>l  4*2 


If  the  inequality  in  Eq.  (  87)  is  not  satisfied,  then  <J>j  and  <(>2 
become  complex  conjugates  and  are  given  by 


(91) 


4>ir  "  4*2 r  " 


2Ii 


(1L.V+  2Li 

Vli  )  21, 


(92) 


4>ii  =  -4>2i  =  -  T 
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ll+&+ 


2!zla 

2Ii 


1/2 


and 


(93) 


9 


where  the  subscripts  r  and  i  refer  to  the  real  and  imaginary  parts 
respectively.  To  complete  the  formulas,  the  moments  of  the  high- 
frequency  waveform  in  Eq.  (  78)  are  given  by 


(94) 


_  aT, 

I,  =  -A  +  BT2  -  — -  +  — 


-aT, 
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(95) 


Bt\  CT2eaT2  CeoTz  DT2e”oT2  De"oT2 

I2  =  — - - - + -  +  - - -  +  - 

■>  2  2 
2  a  a  a  a 


(96) 


I  =  51*  +  li  [De-°Tj  .  CeoT2]  +  ^  [CeoT*  +  De'*T«] 
3  a  or 


+  —  [De'aT2  -  CeoTz  ] 

/v3 


Thus  with  the  exception  of  the  constants  C  and  P,  all  of  the  parameters 
in  the  model  of  the  impulse  response  waveform  of  the  spheroid  in  Eq. 

(  85)  are  determined. 

Consider  first  the  expressions  for  the  creeping  wave  in  the  time 
domain  obtained  in  Appendix  II.  Using  the  approximations  valid  in  the 
neighborhood  of  the  creeping  wave  peak  and  specializing  to  the  time  at 


the  peak,  Peters'  expression  reduces  to 


(97) 
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and  Keller's  result  to 


(98) 
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A  change  of  variable  yields  for  Eq.  (97) 


(99) 


FfW(  a)  = 


A  3,/4 


(2n)jA  (P)7« 


and  for  Eq.  (  98) 


A  39,Z  T(  7/2)  A  39'2  T  (  5/2) 


,9/2 


(  100) 


Ff W(o  ) 


27te  (  tt)  3te  p4  4(ir)s,z  2,fc  p4 


Peters'  result  applies  strictly  to  the  sphere  in  the  form  given  and 
cannot  be  used  for  the  spheroid.  It  will  be  of  interest,  however,  to 
compare  Eq.  (99)  with  the  estimated  peak  obtained  by  Fourier  synthe¬ 
sis  of  the  exact  solution  for  the  sphere.  On  the  other  hand  Keller's 
formula  is  for  a  general  spheroid,  and  while  the  constants  given  apply 
to  the  sphere,  for  the  general  spheroid  only  the  constants  A  and  P 
change  in  Eq.  (  100)  ,  as  well  of  course  as  or.  It  will  also  be  of  interest 
to  see  how  these  results  compare  for  the  sphere. 

A  final  comment  is  in  order  concerning  the  use  of  Keller's  result 
in  Eq.  (  100)  to  predict  the  creeping  wave  peak,  P,  for  the  arbitrarily 
oriented  spheroid.  As  derived  by  Keller,  the  geometrical  theory  of 
diffraction  expression  is  for  axial  incidence  (0=0°)  only.  A  corre¬ 
sponding  expression  for  the  arbitrarily  oriented  spheroid  has  not  been 
given.  The  contention  here,  however,  is  that  for  estimating  the 
creeping  wave  peak  at  an  arbitrary  orientation,  Eq.  (  100)  could  be 
used  provided  the  constants  A  and  P  are  chosen  to  correspond  to  the 
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ellipse  cut  from  the  spheroid  by  the  plane  defined  by  the  incident  elec¬ 
tric  field  vector  and  the  line  of  sight.  This  obviously  is  a  gross  dis¬ 
tortion  of  the  theory  of  geometrical  diffraction,  and  the  author  intends 
in  no  way  to  implicate  Keller  or  others  in  this  assumption.  It  has 
been  seen,  however,  with  the  ramp  response  waveform  model,  that 
very  crude  approximations  in  the  time  domain  can  lead  to  good  fre¬ 
quency  domain  results.  Moreover,  it  will  be  shown  that  an  even 
cruder  assumption  concerning  the  creeping  wave  peak  leads  to  what 
must  be  considered  su-tisfactory  results  for  a  2:1  axial  ratio  spheroid. 
In  any  event,  with  reference  to  the  designated  T.  E.  and  T.  M.  polari¬ 
zations,  for  the  T.  M.  polarization  a  constant  creeping  wave  peak  is 
predicted  with  this  approximation  corresponding  to  the  actual  axial 
ratio  of  the  spheroid.  This  is  so  because  while  the  total  path  length 
changes  (  T2  in  Eq.  (  74) )  ,  the  distance  traveled  along  a  geodesic 
path  on  the  spheroid  remains  constant  for  all  aspects.  However,  for 
the  T.  E.  polarization,  the  ellipse  cut  by  the  above  defined  plane  has 
an  axial  ratio  which  decreases  with  0  (  Fig.  10)  ,  and  is  un;*y  at  0  =  90°. 
The  semiminor  and  semimajor  axis  of  this  ellipse  are  given  by  the 
constants  C\  and  Cz  respectively  in  Eq.  (75). 

In  Fig.  13,  the  creeping  wave  peak  as  calculated  from  Eq.  (  100) 
is  shown  for  a  range  of  axial  ratios  from  1  to  10.  Peters'  estimate 
of  the  creeping  wave  peak  for  a  sphere  (  Eq.  (  99)  )  is  also  shown.  In 
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SPHER0I0  AXIAL  RATIO 


Fig.  13--Peak  magnitude  of  the  creeping  wave  response 
in  the  time  domain  predicted  by  Keller's  spheroid 
diffraction  theory  (axial  incidence)  as  a  function  of  the 

sphe  roid  axial  ratio. 

comparing  these  results  for  the  sphere  with  those  obtained  from  a 
Fourier  synthesis  of  the  exact  solution*  it  must  be  remembered  that 
optical  contributions  are  also  present  in  the  synthesis  picture,  but 
are  not  included  in  the  creeping  wave  estimates.  With  this  in  mind, 
the  estimates  of  the  creeping  wave  peak  in  Fig.  13  for  a  unity  axial 
ratio  seem  reasonable,  if  not  precise. 

An  even  cruder  estimate  of  the  creeping  wave  peak  is  obtained 
by  accepting  the  result  in  Section  A  that  the  creeping  wave  peak  is 
independent  of  the  spheroid  axial  ratio.  In  this  case  P  takes  the  value 
-1.0  approximately,  and  is  independent  of  both  orientation  and 
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polarization.  The  remaining  problem,  with  either  the  above  estimate 
of  P  or  the  geometrical  result,  is  the  choice  of  the  constant  C  in  Eq. 

(  78) .  In  order  to  achieve  a  basic  waveform  with  the  essential  character 
of  the  sphere  response,  it  is  apparent  that  C  must  decrease  with  in¬ 
creasing  axial  ratio  since  the  exponential  in  the  second  term  of  Eq. 

(  78)  should  only  be  strongly  effective  as  the  creeping  wave  peak  is 
approached.  Clearly,  one  can  generate  an  infinite  variety  of  wave¬ 
forms  by  varying  C,  but  several  checks  on  the  magnitude  of  this  pa¬ 
rameter  can  be  used.  This  parameter  is  not,  therefore,  a  curve  fitting 
device  in  the  sense  that  C  is  varied  until  the  best  fit  with  independent 
measured  or  calculated  echo  area  data  is  obtained.  The  first  check  on 
C  is  the  form  of  impulse  response  waveform  obtained;  in  essence  a 
smoothly  decreasing  function  such  as  that  shown  in  Fig.  8  (dashed 
curve)  for  the  exponential  sum  approximation  is  desired.  Note  that 
the  exponential  sum  approximation  in  Reference  ! 2  gives  a  clue  as  to 
the  order  of  magnitude  required  for  C  at  least  for  a  ?.:1  spheroid. 

C  in  this  case  roughly  corresponds  to  the  coefficient  A2  in  Eq.  (  17) 
of  Reference  12,  which  is  of  the  order  of  1  X  10"14,  a  very  small 
number.  A  third  check  on  the  assumed  magnitude  of  C  is  provided 
by  the  coefficients  calculated  for  the  low-frequency  waveform  in  Eq. 

(  83).  The  low-frequency  waveform  is  a  correction  waveform  in  the 
sense  that  it  is  added  to  enforce  the  moment  conditions.  However, 
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If  the  assumptions  concerning  the  similarity  of  the  spheroid  waveform 
to  the  sphere  are  correct,  and  the  choice  of  C  reasonable,  then  very 
little  correction  should  he  needed.  Thus  if  the  magnitude  of  the  coef¬ 
ficients  in  the  low-frequency  waveform  are  such  as  to  seriously  distort 
the  waveform,  then  it  is  an  indie  ition  that  the  assumed  magnitude  of  C 
is  in  error. 

The  basic  model  proposed  for  the  impulse  response  waveform  of 
the  prolate  spheroid  target,  for  the  defined  T.  E.  and  T.  M.  polari¬ 
zations  and  any  orientation,  is  that  given  by  Eq.  (  85).  It  is  deceptive¬ 
ly  simple  in  form,  and  leads,  through  the  Laplace  transform,  to  a 
computationally  simple  frequency-dependent  function.  It  remains  to 
test  the  various  hypothesis  and  assumptions  inherent  in  the  model  both 
as  to  the  form  of  the  response  waveforms  generated  and  ultimately,  of 
course,  with  respect  to  calculated  echo  areas  for  the  target.  These 
tasks  are  undertaken  in  the  next  Section. 

C.  Calculated  and  Measured  Results 
for  a  2:1  Axial  Ratio  Spheroid 

Of  the  two  suggested  models  for  the  general  impulse  response 
waveform  of  the  spheroid,  that  one  utilizing  an  invariant  (with  axial 
ratio,  orientation,  and  polarization)  creeping  wave  peak  taken  from 
the  sphere  is  by  far  the  easiest  to  implement.  It  is  reasonable,  there¬ 
fore,  to  consider  this  model  first.  From  the  exponential  sum 
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approximation  to  the  spheroid  impulse  response  waveform  given  in 
Reference  12,  taking  into  consideration  the  fact  that  for  the  present 
model  the  jump  magnitude  (  F'  does  not  have  an  exponential  decay, 
the  constant  C  in  Eq.  (  78)  is  set  equal  to  1  X  10~12.  The  above 
choice  for  C  was  the  first  one  made.  The  best  proof  of  this,  regard¬ 
less  of  the  resultant  echo  area  calculations,  is  that  neither  the  impulse 
response  waveforms  obtained  for  the  spheroid  nor  the  magnitude  of  the 
damped  sinusoid  correction  are  by  any  means  exactly  what  was  desired 
for  all  orientations  and  polarizations.  It  is  also  clear  that  much  more 
sophisticated  techniques  for  obtaining  both  P  and  C  could  be  devised. 
For  example,  one  could  confine  these  parameters  to  a  particular  range 
of  values  and  then  minimize  the  area  beneath  the  impulse  and  step  re¬ 
sponse  wavefoi  ms  by  some  type  of  testing  procedure.  But  such  so¬ 
phistication  is  accomplished  at  the  expense  of  a  more  complicated 
model.  One  of  the  most  distinctive  features  of  the  time  domain  ap¬ 
proach  is  the  results  which  can  be  obt?  .  ed,  in  terms  of  predicting 
frequency  domain  behavior,  with  very  crude  models.  This  will  be  de¬ 
monstrated  presently  for  the  spheroid  model  using  exactly  the  gross 
approximations  described  above. 

It  is  instructive  first  to  consider  the  magnitudes  of  the  various 
parameters  in  the  impulse  response  waveform  model  described  above 
as  a  function  of  both  orientation  and  polarization.  These  parameters 
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are  given  in  Tables  I  and  II  for  the  T.  E.  and  T.  M.  polarizations  re¬ 
spectively.  The  above  tables  list*  reading  from  left  to  right;  the 
aspect  angle*  specular  magnitude  (A)  ,  jump  magnitude  (  B) ,  time  of 
creeping  wave  peak  (  T2) ,  Rayleigh  coefficient  (  K) ,  exponential  decay 
(a),  exponential  coefficients  C  and  D»  real  and  imaginary  parts  of  the 
correction  exponentials  ( and  $2)  »  and  the  coefficient  of  the  cor¬ 
rection  term  {  E) .  The  symbols  in  the  table  headings  correspond  to 
the  symbols  in  Eqs.  (  78)  and  (  83).  In  Table  IE*  the  corresponding 
moments  of  the  basic  model,  i.  e.  ,  without  the  added  damped  exponen¬ 
tial  of  Eq.  (  83)  are  given,  again  for  both  T.  E.  and  T.  M.  polarizations. 
Note  that  if  the  basic  model  were  correct.  Ij  and  ^  in  Table  III  would 
be  zero,  and  I3  equal  to  the  corresponding  Rayleigh  coefficient. 

Consider  first  the  T.  E.  polarization;  from  Table  I  and  Table  III 
it  is  seen  that  as  the  aspect  angle  increases  from  the  axial  case 
(0  =  0*),  the  basic  waveform  {  Eq.  (  78) )  becomes  first  a  better  and 
then  a  poorer  model  and  the  correction  waveform  ( Eq.  (83))  becomes 
progressively  larger  with  aspect  angle  after  0  =  20°  .  Note  that  the 
uncorrected  moments  of  the  basic  model  are  not  wholly  compatible. 

That  is,  in  terms  of  the  zero  and  first  moments  ( 1}  and  I2) ,  the  model 
is  best  at  approximately  6  =  20°,  in  terms  of  the  second  moment  ( I3 ) , 

8  =  30°  is  the  optimum  model.  In  considering  the  magnitude  of  the 
correction  coefficient  ( last  column  in  Table  I)  it  must  be  remembered 


PARAMETERS  OF  IMPULSE  WAVEFORM  MODEL 


TABLE  n 

PARAMETERS  OF  IMPULSE  WAVEFORM  MODEL 


TABLE  III 

MOMENTS  OF  UNCORRECTED  IMPULSE 
WAVEFORM  MODEL 


that  the  duration  of  the  waveforms  is  decreasing  with  increasing  $  thus 
the  correction  must  effectively  be  made  at  smaller  times.  From  the 
trend  of  the  moments  and  the  correction  coefficient,  it  is  concluded 
that  a  broadening  of  the  creeping  wave  in  the  neighborhood  of  the  peak 
for  aspect  angles  greater  than  30°  and  a  narrowing  of  this  waveform 
below  30°  is  required.  In  the  neighborhood  of  the  creeping  wave  peak, 
rather  small  changes  in  the  slope  of  the  waveform  will  have  a  sub¬ 
stantial  effect  on  the  higher  moments. 


Turning  to  the  T.  M.  polarization  in  Table  II,  a  similar  trend  is 
observed.  Here,  however,  the  zero  and  first  moments  are  smallest 
on  axis  (0=0*)  whereas  the  second  moment  is  best  at  0  =  30*.  For  the 
T.  M.  polarization,  the  moments  of  the  uncorrected  model  are  seen  to 
be  in  error  by  a  greater  amount  for  aspects  near  90°  than  were  those 
for  the  T.  E.  case.  Again,  a  required  broadening  of  the  creeping  wave 
waveform  is  indicated  with  somewhat  greater  broadening  necessary  for 
the  T.  M.  than  for  the  T.  E.  polarization.  It  is  of  interest  next  to 
examine  the  ramp  and  impulse  response  waveforms  predicted  by  this 
model,  with  the  correction  waveform  added,  and  to  compare  the  corre¬ 
sponding  frequency-dependent  function  with  experimental  data.  It  is 
desired  to  compare  the  ramp  response  waveforms  predicted  by  this 
model  with  the  ramp  waveforms  obtained  with  the  earlier  ramp  wave¬ 
form  model  given  in  Eq.  (  72).  This  comparison  is  shown  in  Figs.  14 
through  18  for  both  the  T.  E.  and  T.  M.  polarizations  for  aspect  angles 
Gfrom  0°  to  90°  in  increments  of  10°.  It  is  evident  that  except  for 
small  times,  the  ramp  waveforms  predicted  by  the  two  models  differ 
substantially.  The  differences  are  primarily  a  longer  endurance  and 
higher  magnitude  for  the  impulse  response  derived  model.  Note  also 
that  the  magnitude  of  the  ramp  derived  model  never  exceeds  the  physi¬ 
cal  optics  estimate  whereas  the  impulse  derived  model  does,  particu¬ 
larly  for  aspects  approaching  broadside.  On  the  basis  of  the  ramp 


Fig.  14--Ramp  response  waveforms  of  2:1  axial  ratio 
spheroid  obtained  from  ramp  derived  (  dashed)  and 
impulse  derived  models.  T.  E.  polarization. 

(a)  6  =  0°  ;  (b)  0=  10°  ;  (  c)  0=  20°  ;  ( d)  0=  30° 
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Fig.  15--Ramp  response  waveforms  of  2:1  axial  ratio 
spheroid  obtained  from  ramp  derived  (dashed)  and 
impulse  derived  models.  T.  E.  polarization. 

(a)  0=  40°  ;  (b)  0=  50*  ;  (c)  0=  60°  ;  (  d)  0=  70° 
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ig.  16--Ramp  response  waveforms  of  2*  1  axial  ratio 
spheroid  obtained  from  ramp  derived  (  dashed)  and 
impulse  derived  models. 

(a)  0  =  80°  ;  (  b)  0  =  90  0  ,  T.  E.  polarization 
(  c)  0  =  0°  ;  (  d)  0  =  10 °  ,  T.  M.  polarization 
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Fig.  17--Ramp  response  waveforms  of  2: 1  axial  ratio 
spheroid  obtained  from  ramp  derived  (  dashed)  and 
impulse  derived  models.  T.  M.  polarization. 

(a)  0  =  20°  ;  (  b)  9=  30°  ;  (  c)  0=  40°  ;  (d)  e=  50° 
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Fig.  18--Ramp  response  waveforms  of  2:1  axial  ratio 
sphe  roid  obtained  from  ramp  derived  (  dashed)  and 
impulse  derived  models.  T.  M.  polarization. 

(a)  0=60°  ;  (b)  0=  70°  ;  (c)  0=80°  ;  ( d)  0  =  90° 
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waveform  for  axial  incidence  derived  in  Section  A,  one  would  conclude 
that  the  ramp  derived  model  is  more  nearly  correct,  at  least  for 
aspects  near  axial  incidence.  It  is  not  possible  to  settle  this  question 
completely  with  the  limited  low-frequency  data  available  ( this  point  is 
discussed  later).  One  could,  however,  with  the  techniques  discussed 
in  Chapter  I;  utilize  additional  experimental  data  or  "point  matching" 
solutions  to  estimate  the  mean  square  error  of  the  waveforms.  For  the 
stated  purpose  of  this  study,  that  of  obtaining  an  approximate  general 
model,  it  is  not  necessary  to  resolve  this  question.  It  is  noted  however 
that  some  reservations  concerning  the  ramp  response  waveforms  must 
be  retained  and  that  the  next  investigation  of  the  prolate  spheroid  target 
should  begin  with  the  ramp  response  waveform. 

The  impulse  response  waveforms  predicted  by  the  new  (  Eq.  (  85) ) 
model  are  shown  in  Figs.  19  through  23.  Two  points  are  immediately 
evident  from  these  waveforms.  First,  the  model  has  succeeded  in 
enforcing  the  type  of  creeping  wave  contribution  felt  to  be  necessary. 
Second,  the  damped  sinusoid  correction  term  while  adjusting  the 
creeping  wave  peaks  to  some  extent  has,  in  the  process  of  enforcing 
the  moment  conditions,  distorted  the  waveforms  significantly  for 
earlier  times.  This  is  vividly  evident  for  either  polarization  as  the 
aspect  angle  approaches  90°.  Note  that  the  specular  component  has 
been  omitted  from  these  figures.  It  is  simply  a  negative  delta  function 


Fig.  19--Impulse  response  waveform  of  2:1  axial  ratio 
prolate  spheroid  obtained  from  impulse  derived  model. 

T.  E.  polarization. 

(a)  0  =  0°  ;  (b)  6=  10°  ;  (c)  0=  20°  ;  (  d)  0=  30° 
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Fig.  20--Impulse  response  waveform  of  2:1  axial  ratio 
prolate  spheroid  obtained  from  impulse  derived  model. 

T.  E.  polarization. 

(a)  9=  40°  ;  (b)  0=  50°  ;  (c)  0=  60°  ;  (  d)  6=  70° 
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Fig.  21--Impulse  response  waveform  of  2jl  axial  ratio 
prolate  spheroid  obtained  from  impulse  derived  model, 
(a)  6  =  80°  ;  (  b)  0  =  90°  ,  T.  E.  polarization 
(  c)  0  =  0°;  (  d)  0=10°  ,  T.  M.  polarization 


Fig.  22--Impul.se  response  waveform  of  2:1  axial  ratio 
prolate  spheroid  obtained  from  impulse  derived  model. 

T.  M.  polarization. 

(a)  0=  20°  ;  (b)  0=  30°  ;  (c)  0=  40°  ;  (d)  0  =  50° 
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Fig.  23- -Impulse  response  waveform  of  2;1  axial  ratio 
prolate  spheroid  obtained  from  impulse  derived  modeL 

T.  M.  polarization. 

(a)  0  =  60*  ;  (b)  0=  70°  ;  (  c)  0=  80°  ;  (d)  9  =  90° 
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of  weight  A/c  at  t'  =0  which  is  given  in  Tables  I  and  II  in  the  second 
column.  The  same  conclusions  observed  earlier  from  an  examination 
of  the  model  parameters  in  Tables  I,  II,  and  III  are  confirmed  by  the 
impulse  response  waveforms.  In  each  case  it  is  seen  that  the  cor¬ 
rection  sinusoid  has  increased  the  negative  area  beneath  the  wave¬ 
forms.  Unfortunately,  the  correction  is  not  concentrated  in  the  vi¬ 
cinity  of  the  creeping  wave  peak  and  has,  as  is  apparent,  substantially 
distorted  the  waveforms  for  earlier  times. 

The  final  test  of  the  model  is  a  comparison  of  calculated  and 
measured  results  for  the  backscattering  from  a  prolate  spheroid  target. 
Extensive  data  have  been  measured  on  a  2:1  axial  ratioprolate  spheroid 
for  both  theT.  E.  andT.  M.  polarizations  over  360°of  aspect.  The  meas¬ 
ured  patterns  corresponding  to  this  data  were  reported  in  Reference 
12  and  will  not  be  repeated  here.  In  Reference  12  the  accuracy  of  the 
measured  results  is  estimated  as  -1.  OdB  at  a  level  of  20  dB  less  than 
a  square  wavelength.  Two  comparisons  of  calculated  data  from  the 
model  and  experimental  measurements  are  desired.  The  first  is  the 
echo  area  at  a  fixed  orientation  and  polarisation  as  a  function  of  fre- 

t 

quency,  and  the  second  is  the  echo  area  as  a  function  of  orientation 
for  fixed  frequency  and  polarization.  The  variable  frequency  com¬ 
parisons  will  be  more  extensive.  In  Figs.  24  through  33,  the  predicted 


Fig.  25--Calculated  (  solid  curve)  and  measured  echo  area  of  a  2:. 
axial  ratio  prolate  spheroid  in  dB  >  \2  as  a  function  of  minor  axis 
circumference  in  wavelengths.  0=  10°  . 
a.  T.  E.  polarization. 


Fig.  2 6- -Calculated  (  solid  curve)  and  measured  echo  area  of  a  2: 
axial  ratio  prolate  spheroid  in  dB  >  \2  as  a  function  of  minor  axis 
circumference  in  wavelengths.  0=  20°  . 
a.  T.  E.  polarization. 
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Fig.  ^ 7_ -Calculated  (  solid  curve)  and  measured  echo  area  of  a  2: 
axial  ratio  prolate  spheroid  in  dB  >  \  as  a  function  of  minor  axis 
circumference  in  wavelengths.  0=  30° 
a.  T.  E.  polarization. 
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i*  ig.  c /--calculated  (  solid  curve)  and  measured  echo  area  of  a  h 
axial  ratio  prolate  spheroid  in  dB  >  X2  as  a  function  of  minor  axis 
circumference  in  wavelengths.  0=  30°  . 
b.  T.  M.  polarization. 
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Fig.  28--Calculated  (  solid  curve)  and  measured  echo  area  of  a  2: 
axial  ratio  prolate  spheroid  in  dB  >  X2  as  a  function  of  minor  axis 
circumference  in  wavelengths.  0=  40°  . 
a.  T.  E.  polarization. 


Fig.  28--Calculated  (  solid  curve)  and  measured  echo  area  of  a  2:1 
axial  ratio  prolate  spheroid  in  dB  >  X2  as  a  function  of  minor  axis 
circumference  in  wavelengths.  G=  40°  . 
b.  T.  M.  polarization. 
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Fig.  29 - -Calculated  (  solid  curve)  and  measured  echo  area  of  a  2; 
axial  ratio  prolate  spheroid  in  dB  >  X2  as  a  function  of  minor  axis 
circumference  in  wavelengths.  0=  50°  . 
a.  T.  E.  polarization. 


Fig.  29--Calculated  (  solid  curve)  and  measured  echo  area 
axial  ratio  prolate  spheroid  in  dB  >  X2  as  a  function  of  min 
circumference  in  wavelengths.  0=  50°  . 
b.  T.  M.  polarization. 


Fig.  30 - -Calculated  (  solid  curve)  and  measured  echo  area  of  a  2: 
axial  ratio  prolate  spheroid  in  dB  >  X2  as  a  function  of  minor  axis 
circumference  in  wavelengths.  0=  60°  . 
a.  T.  E.  polarization. 


Fig.  3 1  -  -Calculated  (  solid  curve)  and  measured  echo  area  of  a  2:1 
axial  ratio  prolate  spheroid  in  dB  >  as  a  function  of  minor  axis 
circumference  in  wavelengths.  0=  70°  . 
a.  T.  E.  polarization. 
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Fig.  33--Calculated  (  solid  curve)  and  measured  echo  area  of  a  2:1 
axial  ratio  prolate  spheroid  in  dB  >  X2  as  a  function  of  minor  axis 
circumference  in  wavelengths.  0=  90°  . 
a.  T.  E.  polarization. 


axial  ratio  prolate  spheroid  in  dB  >  X2  as  a  function  of  minor  axis 
circumference  in  wavelengths.  0=  90°  . 
b.  T.  M.  polarization. 
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and  measured  echo  area  of  the  2:1  prolate  spheroid  are  compared  for 
aspect  angles;  0,  from  0°  in  increments  of  10°  to  90°  as  a  function  of 
the  minor  axis  circumference  of  the  spheroid  in  wavelengths.  In  each 
figure,  the  calculated  data  is  the  solid  curve.  Part  a  of  each  figure  is 
for  the  T.  E.  polarization  and  part  b  for  the  T.  M.  polarization.  It  is 
evident  from  these  figures  that  considering  the  crudeness  of  the  theo¬ 
retical  model,  substantial  agreement  with  the  measured  data  has  been 
obtained.  Minor  disagreements  are  apparent,  as  in  the  location  of 
maxima  and  minima  for  aspects  of  0°  and  10°,  but  the  general  indi¬ 
cation  is  that  the  model  rather  accurately  predicts  the  echo  area  of  the 
2:1  spheroid  for  arbitrary  orientation  and  linear  polarization. 

The  question  concerning  the  ramp  response  waveforms  of  the 
ramp  derived  and  impulse  derived  models  is  partially  resolved  by  the 
comparison  with  measured  echo  data.  It  was  found  that  for  aspect 
angles  less  than  40°  (for  both  polarizations)  calculations  from  the 
ramp  derived  model  were  in  somewhat  better  agreement  with  the 
measured  data  for  minor  circumferences  less  than  one  wavelength. 

For  aspect  angles  greater  than  40°,  there  was  little  to  choose  between 
the  two  models.  In  the  limit  of  low  frequencies  the  models  agree,  as 
dictated  by  the  enforced  moment  conditions.  For  aspects  less  than 
40°  and  minor  axis  circumferences  less  than  one  wavelength,  the 
calculations  shown  in  Figs.  24-28  are  from  the  ramp  derived  model. 
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Turning  to  the  second  test,  in  Figs.  34  through  37,  the  measured 
and  calculated  echo  and  phase  patterns  of  a  2:1  axial  ratio  spheroid  for 
the  T.  M.  polarization  are  shown  for  minor  axis  circumferences  of  1.  5, 
2.  0,  2.  5,  and  3.  0  wavelengths  respectively.  In  these  figures  the 
measured  data  are  shown  as  the  solid  curve  (part  a  of  each  figure  is 
the  echo  data,  part  b  the  phase  data).  The  phase  data  are  referenced 
to  the  center  of  the  sphe  roid.  It  is  estimated  that  the  echo  area  data 
have  an  error  of  -  1.0  dB  at  20  dB  below  a  square  wavelength,  and 
the  phase  data  an  error  of  -  10  degrees.  The  agreement  between 
calculated  and  measured  data,  while  perhaps  not  as  striking  as  for 
the  variable  frequency  curves,  is  still  reasonably  good. 

Calculations  of  the  creeping  wave  peak  predicted  by  Keller 's 
geometrical  theory  of  diffraction  result  were  shown  in  Fig.  13  as  a 
function  of  the  spheroid  axial  ratio.  Two  points  should  be  noted  from 
this  curve.  First,  Keller's  result  for  the  degenerate  case  of  the 
sphere  is  in  error,  as  is  seen  in  Appendix  II  from  the  Fourier  synthe¬ 
sis  figures,  and  would  have  to  be  shifted  in  magnitude  to  be  used. 
Second,  for  the  2:1  spheroid  the  change  in  creeping  wave  peak  from 
the  sphere  case  is  so  slight  that  a  second  model,  for  the  case  of  a  2:1 
spheroid,  does  not  seem  to  be  warranted.  Suggestions  for  the  use  of 
the  curve  in  Fig.  13  for  higher  axial  ratio  spheroids  are  given  in  the 
next  Section.  However,  it  is  noted  here  that  with  estimates  of  the  slope 
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(b) 


Fig.  34- -Measured  (  solid)  and  calculated  echo  and  phase  patterns 
of  a  2:1  axial  ratio  prolate  spheroid.  Minor  axis  circumference 
equal  1.  5  wavelengths,  (a)  Echo  area.  (  b)  Phase. 
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Fig.  35--Measured  (  solid)  and  calculated  echo  and  phase  patterns 
of  a  2:  axial  ratio  prolate  spheroid.  Minor  axis  circumference 
equal  2.  0  wavelengths.  (  a)  Echo  area.  (  b)  Phase. 
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Fig.  36--Measured  (  solid)  and  calculated  echo  and  phase  patterns 
of  a  2:1  axial  ratio  prolate  spheroid.  Minor  axis  circumference 
equal  2.  5  wavelengths,  (a)  Echo  area.  (  b)  Phase. 
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Fig.  37--Measured  (  solid)  and  calculated  echo  and  phase  patterns 
of  3  2:1  axial  ratio  prolate  spheroid.  Minor  axis  circumference 
equal  3.0  wavelengths,  (a)  Echo  area.  (  b)  Phase. 
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of  the  creeping  wave  on  either  side  of  the  peak  from  Keller's  theory, 
additional  parameters  could  be  added  to  the  impulse  waveform  model. 
Realistically,  there  seems  little  advantage  in  this  added  complexity 
until  failures  of  the  simpler  model  are  noted  or  a  more  exact  asymp¬ 
totic  estimate  of  the  spheroid  creeping  wave  is  obtained. 

Finally  it  is  desired  to  list  some  results  obtained  from  the  im¬ 
pulse  response  waveform  model  of  the  2:1  spheroid  which  may  be  of 
use  to  other  investigators.  The  primary  concern  is  not,  as  was  noted 
in  Chapter  I,  the  details  of  the  impulse  response  waveform  of  the  de¬ 
rived  model,  but  is  instead  the  predicted  ramp  response  waveform. 

A  reasonable  approximation  to  this  ramp  waveform  can  be  generated 
by  Fourier  synthesis  from  the  first  nine  harmonics  of  the  phasox  re¬ 
sponse,  and  a  comparison  of  the  ramp  waveform  predicted  by  other 
theories  with  those  given  here  would  be  extremely  useful.  To  this  end, 
the  predicted  phasor  response  for  both  the  T.  E.  and  T.  M.  polarizations 
at  aspects  of  0°  (  10°)  90°  for  minor  axis  circumferences  from  0.  25 
(  0.  25)2.  25  wavelengths  are  tabulated  in  Table  IV.  The  phase  reference 
for  these  data  is  at  the  point  the  incident  wave  first  strikes  the  spheroid. 
Note  in  Table  IV  that  the  real  and  imaginary  parts  are  normalized  such 
that  the  cross  section  normalized  by  the  shadow  cross-sectional  area 
■rra2  is  given  by  the  sum  of  the  squares  of  the  real  and  imaginary  parts. 


TABLE  IV 

■ira2  =  1  R2  +  I2  1  T.  E.  POLARIZATION  (  Minor  Axis  Circumference  in  Wavelengths) 
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TABLE  IV  (  Contd.  ) 

T.  M.  POLARIZATION  (  Minor  Axis  Circumference  in  Wavelengths) 


156 


O' 

m 

r- 

00 

no 

O' 

00 

O' 

00 

H4 

no 

O) 

H4 

05 

05 

O' 

h- 

in 

m 

no 

ro 

* 

H4 

ro 

no 

no 

m 

05 

r- 

O' 

»H 

h- 

00 

no 

r*- 

m 

oo 

no 

o 

o 

pH 

no 

— 1 

in 

o 

05 

© 

v0 

H 

05 

no 

m 

no’ 

• 

o 

• 

o 

• 

o 

• 

o 

s 

a 

• 

H 

O 

• 

pH 

• 

o 

• 

H 

o 

• 

pH 

• 

O 

• 

pH 

• 

o 

C 

i 

i 

H 

■ 

i 

i 

1 

i 

i 

i 

i 

I 

1 

i 

m 

no 

■<* 

m 

no 

o 

05 

H4 

00 

ro 

05 

ro 

00 

o 

00 

pH 

O 

00 

o 

o 

O 

no 

■<* 

O' 

in 

in 

00 

no 

O 

03 

05 

no 

n- 

no 

o 

m 

sO 

m 

in 

ro 

ro 

no 

05 

O 

O 

H 

m 

no 

i'- 

no 

no' 

• 

o 

o 

o 

• 

o 

• 

O 

o 

• 

O 

• 

o 

o 

• 

o 

• 

o 

o 

• 

pH 

• 

o 

t 

pH 

• 

o 

• 

^H 

• 

o 

i 

i 

i 

i 

i 

1 

i 

1 

1 

1 

1 

• 

i 

i 

O' 

o 

no 

m 

m 

* 

o 

n- 

o 

ro 

00 

in 

05 

O' 

00 

in 

o 

no 

00 

no 

05 

ro 

05 

in 

05 

H4 

pH 

05 

no 

no 

no 

in 

o 

O 

© 

pH 

no 

no 

m 

ro 

O' 

ro 

no 

O' 

o 

H4 

ro 

05 

• 

• 

o 

• 

O 

o 

• 

O 

• 

O 

• 

O 

• 

o 

• 

o 

o 

• 

o 

• 

pH 

• 

O 

• 

pH 

o 

no* 

• 

O 

• 

no 

• 

o 

1 

t 

i 

1 

i 

i 

i 

i 

i 

i 

i 

i 

r- 

pH 

05 

no 

o 

pH 

05 

o 

00 

CO 

o 

05 

ro 

O' 

o 

sO 

O 

no 

o 

o 

O' 

© 

o 

ro 

H4 

m 

O' 

m 

05 

05 

H4 

O 

ro 

05 

no 

cO 

03 

in 

O' 

no 

o 

no 

pH 

O 

no 

H 

ro 

H4 

H4 

r- 

05 

O 

r- 

ro 

00 

H4 

. 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

, 

• 

• 

• 

• 

— i 

o 

o 

pH 

O 

pH 

o 

pH 

o 

pH 

o 

pH 

pH 

pH 

pH 

H 

pH 

1 

i 

i 

i 

1 

i 

1 

' 

1 

1 

i 

i 

i 

co 

vO 

ro 

H* 

05 

no 

no 

CO 

05 

ro 

pH 

O) 

o 

r- 

O' 

00 

o 

in 

O 

in 

O' 

n- 

O' 

H 

h- 

o 

no 

n- 

O' 

O' 

o 

o 

ro 

H4 

05 

no 

sO 

05 

in 

05 

vO 

r? 

05 

H4 

05 

H4 

m 

m 

05 

h- 

05 

n- 

05 

• 

"  pH 

d 

• 

o 

o 

o 

o 

o 

• 

o 

o 

• 

O 

• 

o 

d 

• 

o 

o 

o 

• 

o 

• 

O 

• 

o 

• 

O 

i 

1 

i 

i 

i 

i 

1 

i 

1 

1 

i 

i 

i 

1 

1 

l 

i 

1 

no 

o 

H* 

CO 

h4 

pH 

05 

ro 

pH 

05 

00 

no 

pH 

ro 

00 

H4 

o 

00 

pH 

no 

r- 

ro 

H4 

pH 

O' 

sO 

no 

05 

no 

O' 

ro 

no 

05 

05 

00 

° 

O 

o 

O 

pH 

^H 

no 

m 

no 

m 

ro 

00 

ro 

pH 

CO 

in 

CO 

05 

ro 

• 

• 

o 

o 

o 

• 

o 

• 

o 

• 

O 

• 

o 

• 

O 

o 

• 

O 

o 

o 

pH 

• 

o 

• 

pH 

• 

O 

• 

pH 

• 

O 

1 

i 

1 

1 

1 

i 

i 

i 

i 

O'  'O 

O' 

in 

m 

H4 

pH 

o 

H4 

H4 

ro 

ro 

© 

no 

ro 

no 

in 

no 

m 

O' 

ro 

00 

pH 

n- 

o 

pH 

00 

00 

O' 

O' 

H4 

m 

no 

Is- 

O' 

o 

O' 

o 

o 

O 

ro 

o 

m 

no 

n- 

ro 

pH 

m 

05 

00 

00 

O' 

• 

o 

• 

o 

o 

o 

o 

• 

• 

O 

• 

H 

• 

o 

• 

•H 

• 

o 

• 

H 

• 

O 

no" 

• 

o 

no 

o 

no* 

• 

O 

1 

1 

i 

i 

i 

1 

1 

1 

i 

1 

i 

i 

i 

1 

i 

1 

m 

O) 

o 

ro 

H4 

ro 

in 

o 

no 

00 

O 

no 

o> 

no 

m 

o 

o 

in 

o 

no 

H4 

00 

pH 

t'- 

^H 

H4 

H4 

no 

00 

•O 

h- 

H* 

o 

r- 

m 

00 

00 

H4 

O' 

O 

00 

in 

O' 

O' 

o 

H4 

-4 

05 

o 

o 

• 

O 

o 

• 

o 

o 

o 

• 

pH 

o 

o 

• 

pH 

o 

pH 

• 

pH 

no 

• 

pH 

• 

no 

i 

i 

I 

i 

1 

i 

i 

1 

i 

i 

i 

1 

i 

i 

i 

i 

1 

o 

^H 

in 

03 

00 

r- 

in 

O' 

ro 

no 

H4 

no 

^H 

pH 

M 

no 

O 

in 

r- 

pH 

03 

m 

03 

O'  '0 

H4 

05 

00 

^H 

00 

N 

O' 

no 

pH 

no 

no 

rj 

no 

no 

no 

no 

no 

ro 

no 

ro 

no 

H4 

nj 

* 

"VI 

• 

o 

o 

o 

o 

• 

o 

O 

o 

o 

o 

o 

O 

O 

o 

d 

O 

o 

o 

• 

o 

o 

i 

1 

i 

i 

i 

i 

I 

| 

i 

R 

I 

1 

o 

o 

o 

o 

o 

o 

o 

o 

o 

II 

o 

H 

no 

ro 

H4 

m 

05 

r- 

00 

O' 

sssjSdQ  ui  0 


Calculated  data  on  a  prolate  spheroid  with  a  10:1  axial  ratio  for 
axial  incidence  (  6  =  0°)  have  been  given  by  Siegel,  et  al.  [  38] .  These 
data  are  shown  in  Fig.  38  as  the  solid  curve.  Also  shown  in  Fig.  38 
are  calculations  for  this  target  from  the  ramp  response  derived  model. 
These  data  are  designated  by  crosses.  An  impulsive  response  wave¬ 
form  model  for  this  target  was  also  given  by  Kennaugh  and  Cosgriff[  3]  , 
and  showed  reasonable  agreement  with  the  calculated  data.  For  such 
a  large  axial  ratio  spheroid,  it  is  clear  that  the  simple  model  in  Eq. 

(  78)  is  not  a  good  approximation.  As  shown  by  the  model  in  Reference 
3,  the  creeping  wave  contribution  for  such  a  large  axial  ratio  is  well 
approximated  by  a  very  narrow  rectangular  pulse.  It  is  not  practical 
to  obtain  such  a  waveform  with  one  exponential  term  as  in  Eq.  (78). 
Therefore,  for  large  axial  ratio  spheroids,  the  ramp  response  derived 
model  should  be  used.  This  point  is  discussed  in  the  next  section. 

D.  Conclusions 

It  was  mentioned  earlier  in  this  Chapter  that  certain  of  the  ap¬ 
proaches  and  procedures  followed  in  attacking  the  prolate  spheroid 
problem  were,  with  the  experience  gained  herein,  not  necessarily 
those  recommended  either  as  an  approach  to  future  problems  or  a 
continuation  of  the  spheroid  study.  One  must  remember  that  the  use 
of  transient  and  impulse  response  waveform  estimates  in  electro¬ 
magnetic  scattering  problems  is  still  an  embryonic  art.  Significant 


Fig.  38--Axial  (  0  =  0°)  echo  area  of  a  10:1  axial 
ratio  prolate  spheroid. 
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contributions  have  been  made,  as  cited  in  Chapter  I,  both  to  the  so¬ 
lution  of  new  scattering  problems ,  and  providing  additional  insight  to 
the  scattering  mechanism.  In  a  sense,  however,  as  discussed  in 
Chapter  I,  these  contributions  have  been  made  with  a  limited  ap¬ 
proach  in  that  a  systematic  procedure  for  application  of  the  concepts 
has  not  been  developed.  Unfortunately,  these  comments  hold  true 
for  the  material  reported  in  this  Chapter.  A  general  model  has  been 
developed  for  the  electromagnetic  backscattering  by  a  perfectly  con¬ 
ducting  prolate  spheroid  target  for  arbitrary  linear  polarization  and 
arbitrary  orientation  which  has  been  shown  to  be  in  reasonable  agree¬ 
ment  with  measured  data  for  a  spheroid  of  2:1  axial  ratio.  But  despite 
this  success,  the  void  with  regard  to  a  systematic  procedure  remains. 
The  conclusions  from  this  Chapter  therefore  are  of  two  types;  those 
pertaining  strictly  to  the  spheroid  model,  its  deficiencies,  possible 
improvements,  and  applicability  to  other  untested  axial  ratios,  and 
those  pertaining  to  the  application  of  a  time  domain  approach  to  other 
problems.  In  this  Section,  the  discus sion  will  be  confined  to  the 
former  type,  i.  e.  ,  to  the  specific  model  developed  for  the  conducting 
prolate  spheroid  target.  What  insight  has  been  gained  in  developing 
a  systematic  approach  to  the  general  problem  has  been  incorporated 
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into  Chapter  l,  primarily  into  Section  C  of  Chapter  I,  on  the  analyti¬ 
cal  and  experimental  approximations  of  transient  and  impulse  re¬ 
sponse  waveforms. 

Consider  first  the  deficiencies  of  the  model,  and  admittedly  it 
has  many.  These  deficiences  are  primarily  in  terms  of  failures  to 
relate  the  geometrical  properties  of  the  spheroid  to  the  corresponding 
detail  of  the  model.  It  is  certainly  clear  that  the  different  radii  of 
curvature  in  the  two  planes,  the  changing  shadow  boundary  with  re¬ 
spect  to  the  line  of  sight,  and  the  polarization  properties  of  the 
spheroid,  must,  in  any  complete  analysis,  materially  affect  the  tran¬ 
sient  waveform.  Some  of  these  effects  may  have  been  partially 
accounted  for  by  the  conditions  imposed  on  the  model.  But  these  con¬ 
ditions  were  not  explicitly  related  to  the  geometrical  detail,  conse¬ 
quently  when  rather  large  changes  in  the  geometry  are  made,  as  with 
a  large  axial  ratio  spheroid  or  a  different  smooth  geometry,  there  is 
no  assurance  that  the  same  conditions  will  continue  to  effect  adequate 
compensation.  These  comments  apply  both  to  the  illuminated  and 
shadowed  side  of  the  sphe  roid.  On  the  illuminated  side,  corrections 
are  needed  to  the  physical  optics  approximation  which  are  related  to 
the  geometry.  For  the  conducting  sphere,  Weston{  44]  has  supplied 
such  an  additional  term  from  an  asymptotic  expansion  of  the  optics 
field,  and  asymptotic  estimates  might  be  utilized  for  the  spheroid 
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in  a  future  analysis.  Similar  comments  apply  to  the  creeping  wave 
portion  of  the  model.  Despite  the  success  of  the  comparison  with 
measured  data,  the  exact  nature  of  the  creeping  wave  contribution  is 
still  unknown.  How  the  contribution  builds  to  the  secondary  peak,  the 
slope  of  the  waveform  on  either  side  of  this  peak,  and  particularly 
how  the  geometry  of  the  spheroid  is  related  to  these  effects  is  not 
known.  In  this  regard,  the  noncausal  nature  of  the  creeping  wave 
contribution  in  the  time  domain  (  see  Appendix  II)  predicted  by  asymp¬ 
totic  estimates  and  the  attendant  questions  would  seem  to  indicate  that 
a  complete  understanding  of  the  sphere  scattering  in  all  detail  has  not 
been  attained  as  yet. 

With  regard  to  the  estimates  of  the  time  domain  waveforms,  it 
must  be  concluded  in  light  of  the  above  remarks  that  the  detailed 
features  of  the  waveform  may  bear  scant  resemblence  to  the  true  re¬ 
sponse.  This  points  out  an  extremely  useful  and  at  the  same  time 
frustrating  nature  of  the  transient  waveforms.  Very  crude  approxi¬ 
mations  in  the  time  domain  lead,  in  many  cases,  to  estimates  of  the 
phasor  response  which  are  out  of  all  proportion  to  the  exactness  of 
the  approximation.  In  establishing  the  correctness  of  the  time  domain 
waveform  then,  comparisons  of  the  calculated  and  measured  echo  area 
of  an  object  are  of  limited  usefulness.  This  is  particularly  true  if  the 
measured  data  are  confined  to  a  small  portion  of  the  spectrum.  It  is  felt 
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that  the  ramp  response  waveforms  predicted  by  the  ramp  derived 
model  are  a  reasonable  estimate  of  the  true  response.  No  such  confi¬ 
dence  is  expressed  in  the  ramp  and  related  step  and  impulse  response 
waveforms  of  the  impulse  derived  model,  however,  except  in  the 
gross  features  of  the  response. 

With  regard  to  the  utilization  of  the  impulse  derived  model  for 
other  axial  ratio  spheroids,  the  following  comments  are  offered. 
Keller's  model  for  the  creeping  wave  contribution  is  clearly  inade¬ 
quate  since  it  fails  even  in  the  degenerate  case  of  the  sphere  to  pre¬ 
dict  the  correct  detail  of  the  waveform.  However,  with  no  present 
alternative,  it  is  suggested  that  for  higher  axial  ratios  the  peak  of 
the  creeping  wave  shown  in  Fig.  1  3  be  utilized  after  first  scaling 
the  results  so  that  they  agree,  for  the  degenerate  sphere  case,  with 
the  Fourier  synthesis  result  using  the  exact  solution.  Thus  for  the 
T.  E.  polarization  a  varying  magnitude,  P,  for  the  creeping  wave 
peak  should  be  used  while  for  the  T.  M.  case  it  remains  fixed  by  the 
axial  ratio  of  the  spheroid.  The  predicted  ramp  response  waveform 
and  the  magnitude  of  the  correction  term  remain  clues  to  the  esti¬ 
mate.  Any  radical  departure  of  the  form  of  the  ramp  waveform, 
aside  from  a  general  stretching,  should  be  suspect  as  should  the 
magnitude  of  the  correction  term.  The  parameter  C  should  be  de¬ 
termined  from  tests  of  the  impulse  response  waveform  at  0=  0°  and 
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90®  .  The  basic  premise  being  to  retain  the  gross  resemblance  of  the 
waveform  to  that  of  the  sphere.  Both  Keller's  result  and  simply  an 
examination  of  the  geometry  indicate  that  as  the  axial  ratio  of  the 
spheroid  increases,  the  attenuation  of  the  creeping  wave  should  de¬ 
crease.  In  the  time  domain  this  would  indicate  that  the  waveform  in 
the  neighborhood  of  the  creeping  wave  peak  should  become  narrower 
with  increasing  axial  ratio.  As  mentioned  earlier,  the  impulse  re¬ 
sponse  waveform  of  a  10:1  axial  ratio  spheroid  given  in  Reference  3 
demonstrated  this  effect.  As  the  axial  ratio  of  the  spheroid  increases, 
the  creeping  wave  contribution  apparently  approaches  a  very  narrow, 
high  spike,  with  negligible  value  except  in  the  immediate  vicinity  of 
the  creeping  wave  peak.  It  is  clear  that  the  simple  impulse  response 
derived  model  in  Eq.  (  85)  is  not  practical  for  large  axial  ratio 
spheroids.  For  example,  a  value  of  1  X  lO"4’*  is  required  for  C  to 
achieve  a  reasonable  waveform  for  a  10:1  axial  ratio  spheroid.  For 
such  large  axial  ratios,  models  such  as  that  given  in  Reference  3  or 
the  ramp  response  derived  estimate  in  Section  B  should  be  used. 

For  axial  ratios  greater  than  2:1  or  3:1,  it  will  probably  not  be  pos- 
sible  to  use  a  single  value  for  the  parameters  C  and  P  for  all 
orientations  and  polarizations.  The  model  loses  much  of  its  attractive¬ 
ness  in  this  case  since  individual  adjustments  for  each  orientation  and 
polarization  will  be  required.  The  impulse  response  derived  model 
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has  not  been  tested  for  axial  ratios  greater  than  2:1.  It  is  evident 
however  that  utilization  of  the  model  becomes  more  difficult  with 
increasing  axial  ratio  and  is  probably  prohibitively  difficult,  with¬ 
out  modifications,  for  axial  ratios  greater  than  5:1. 

It  is  readily  apparent  that  tbit,  3tudy  in  the  time  domain  of  the 
prolate  spheroid  target  has  raised  more  questions  than  it  has 
answered.  A  novel  result  has  been  obtained;  a  model  predicting  the 
backscattering  from  the  prolate  spheroid  target  for  arbitrary  orien¬ 
tation,  and  linear  polarization  of  the  incident  field  has  been  derived; 
and  the  general  validity  of  the  model  demonstrated  for  the  case  of  a 
2:1  spheroid.  In  the  opinion  of  the  author,  this  model,  with  judicious 
application,  may  continue  to  demonstrate  reasonable  agreement  with 
other  measured  or  calculated  data  for  axial  ratios  less  than  5:1.  In 
fairness  both  to  this  study  and  to  the  time  domain  approach  in  general, 
one  should  be  cognizant  of  the  numerous  frequency  domain  studies  of 
the  spheroid  target  which  have  had  far  less  success.  However,  it  is 
clear  than  an  understanding  of  the  scattering  by  the  spheroid,  and 
consequently  an  ability  to  project  these  results  to  other  shapes,  is 
meager.  The  most  pressing  requirement  at  present  is  a  detailed 
model  for  the  creeping  wave  contribution  from  smooth  surfaces  in 
the  time  domain.  The  asymptotic  high  frequency  estimates  of  this 
response,  when  transformed  to  the  time  domain,  offer  clues  in  this 
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direction,  particularly  with  regard  to  the  behavior  of  this  reponse  in 
the  neighborhood  of  its  peak.  A  feeling  persists  however  that  the 
ultimately  most  fruitful  attack  on  the  problem  will  be  one  undertaken 
in  the  time  domain.  Considering  the  relative  simplicity  of  the  real, 
time-dependent  waveforms  corresponding  to  complex  frequency-de- 
pendent  functions;  it  does  not  seem  unrealistic  to  suggest  that  such 
time-dependent  functions  may  be  more  intimately  and  simply  related 
to  the  geometry  of  the  scattering  object.  This  comment  is  felt  to  be 
pertinent  in  general,  and  not  restricted  to  the  creeping  wave  problem. 

It  is  apparent  that  one  could  make  a  number  of  adjustments  on 
the  impulse  response  derived  model  to  minimize  the  low  frequency 
corrections  required.  By  examining  the  ramp  and  impulse  response 
waveforms  and  the  uncorrected  moments  as  the  parameters  P  and  C 
are  varied,  an  improved  model  for  each  aspect  and  polarization 
could  be  obtained.  However,  as  indicated  in  Chapter  I,  such  an  ap¬ 
proach  does  not  lead  to  an  "understanding"  of  the  scattering  by  a 
prolate  spheroid  target.  It  is  recommended  that  future  research  on 
this  problem  begin  with  improved  estimates  of  the  ramp  response 
waveform  using  measured  and  calculated  low  frequency  data.  The 
"point  matching"  solutions  mentioned  previously  appear  to  offer  a 
convenient  and  accurate  method  for  securing  the  necessary  phasor 


response  data. 


APPENDIX  I 

RAYLEIGH  SCATTERING  BY  SPHEROIDS 


It  is  well  known  that  for  objects  sufficiently  small  in  terms  of 
wavelength,  the  scattering  cross  section  varies  as  the  sixth  power 
of  dimension  and  the  fourth  power  of  frequency.  The  purpose  of 
these  notes  is  to  develop  explicit  expressions  for  the  factor  of  pro¬ 
portionality  for  scattering  from  spheroids.  Stevenson(  42 ]  ,  in  his 
excellent  paper,  has  provided  a  general  solution  for  the  coefficients 
in  a  power  series  expansion  in  the  wavenumber  of  the  scattered  field. 
We  are  concerned  here  with  the  coefficient  of  the  wavenumber  squared 
or  the  so-called  Rayleigh  term  for  the  case  of  spheroids.  We  wish  to 
couch  our  development  in  terms  of  electric  and  magnetic  dipoles 
aligned  with  the  principal  axes  of  the  body  whose  relative  strengths 
are  determined  by  the  aspect  and  polarization  of  the  incident  field. 

Following  Stevenson,  let  f,m,n  be  the  direction  cosines  of  the 
direction  of  propagation  of  the  incident  field;  l\ ,  mj  ,  ni  be  the  di¬ 
rection  cosines  of  the  incident  electric  vector;  and  f2,m2,  n2  be  the 

loot 

direction  cosines  of  the  incident  magnetic  vector.  With  e” 
understood  the  incident  field  is  given  by[  42] 
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(101)  El  =  (  ij  ,  mi  ,  n)  e 


ik(  lx  +  my  +  nz) 


and 


H1  =  (  i2,  m2,  n2)  e 


ik(  lx  +  my  +  nz) 


and  the  scattered  fields  by[  42] 


(  102) 


'0 


s  P 
=  =  [ ] 


8P  l 

+ 


_ 8P_ 

d0s  sin  0g  dtp 


S 


ikR 

k 

0 _ 

kR 


and 


E? 


=  -  H8  =  — - 
6  sin 


8P 


8P 


in  08  8<t>s  88 


sj 


ikR 


kR 


where  R,  0S,  <j)s  are  the  sphe  rical  coordinates  of  the  field  point 


(103)  P  =  (ka)S  [Kilja  +  K2mi  0  +  Ksn1Y] 


and 


r\j 

P  =  (ka)3[Kji2  a  +  K2m2  0+  Kj  n2  y] 


where  a,  0,  y  are  the  direction  cosines  of  the  radius  vector  to  the 
field  point.  The  constants  Ki .  .  ,  K*.  .  ,  are,  respectively,  the 
electric  and  magnetic  coefficients  which  are  a  function  only  of  the 
body  and  independent  of  the  incident  field.  Note  that  in  the  "a"  in 
Eq.  (  103)  is  the  principal  axis  along  the  x  and  y  coordinates 
(  Fig.  10). 
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(  104)  Kj  = 


Z(€-l) 


3a5  [(«  - 1 )I.  + 

L  a  a*bj 


»  = 


2(H-1) 


3a*  fC|a-l)I  +  *  1 

L  a  a2bJ 


K*  = 


2(<-l) 


3a>  f(€-l)Ia  +-4 
.  a  d. 


K,  = 


2m- 1) 


and 


K,  = 


2(« -1) 


3aJ  f{« -l)Ib  +4- 

L  a'bj 


-i-J 


where  €  and  y.  are  the  constitutive  parameters  of  the  spheroid.  Note 
that  we  have  altered  Stevenson's  expressions  slightly  to  conform  to 
our  desired  dipole  development.  The  expressions  Ia,  Ib  are  related 
by 

2 

(105)  2Ia  +  Ib=  — 

azb 


and 


2a*  Ia  +  b*  Ib  =  I  , 


where  now  we  must  differentiate  between  two  cases.  For  the  prolate 
spheroid  (  a  <  b) 


9 


where  t]  =  n/  1  -(  a/b)  ^  . 


For  the  oblate  spheroid  ( a  >  b) 


where 


(  108) 


and 


r|  -  \J  l-(  b/a)  *  .  Thus  for  the  prolate  spheroid  we  obtain 

2(«-l)  _ 


K,  =  K2  = 


3^[2^  (ttS 

»  L 


+  z/TTr?" 


K,  = 


2(«-J) 


<«-U 

— 


Ki  =  K2  - 


2(p-l) 


3  L2^  |20  -i2  -  ( 1  -  i2)S/2  * n  (tt^) 


M, 

+  2/U?- 

\l-h/ 

4 

4 

*.  = 


2(u-l) 


(k-i) 

“I — 


P  _ 

zj  i-ti2  (n2- 

u 


1)  +  iiilL — t 


.2)3/  2 


-  (£)]  *  ^ 
«  « 


For  the  perfect  conducting  case  (|  e  |  —  ®  ,  p  —  0),  the  above  reduce 
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l 


For  the  oblate  spheroid 


For  a  perfect  conductor  the  above  become 


Now  for  the  general  case, 


(  H2) 


a  =  sin  9g  cos  <J>S  , 

(3  =  sin  0S  sin  4>s  ,  and 
■y  =  cos  0g 


Let  the  coordinates  of  the  source  point  be  R,  0^,  4>i  >  then 
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(  113)  I  =  -sin  0j  cos  4^  » 

m  =  -sin  0i  sin  ,  and 

rj  =  -cos  0^ 

There  is  no  loss  of  generality  if  we  restrict  the  direction  of 
incidence  to  the  y,z  plane  (4>j  =  tt/2)  whence 

(114)  /  =  0, 

m  =  -sin  0j  ,  and 

T|  -  -COS  0- 

We  now  consider  two  particular  polarizations  of  the  incident 
field  from  which,  by  superposition,  the  solution  for  any  arbitrary 
polarization  can  be  obtained. 

T.  E.  Incident  electric  vector  perpendicular  to  yz  plane. 


*1  - 

1. 

0 

*2  - 

o 

O 

mi  = 

0. 

0 

m2  = 

-COS  0} 

Til  = 

0. 

0 

Tl2  = 

sin  0j 

T.  M.  Incident  electric  vector  parallel  to  yz  plane. 


(116) 


=  0.  0 

mi  =  -cos  0. 
ni  =  sin  0^ 


i2  =  -1.0 
m2  =  0.  0 
Tl2  =  o.  0 


Let  us  define  a  normalization  cuch  that 


(117) 


Then  the  normalized  scattered  fields  are  given  by 
T.  E. 


(  118)  E0ti=  =  2(ka)*[Kicos  08  cob  <J>8  -  cob  0{  Coa  4>8J 

and  /v 

=  2(ka)2[-Kiflin  4>8  +  K2  coa  cos  08  sin  4>s 

+  Kj  sin  0^  sin  0flJ ; 


T.  M. 


(119)  Een  =  =  2(ka)2  [-K2  cos  cos  0,  sin  <|>8  -  Kj  sin  0t  sin  0 

+  sin  4>SJ 


and 


JE^ri  =  -  Hen  =  2(ka)2  [-K2  cob  0*  cob  <J>8  +  Ki  cos  08  cob  <t>8]  . 


Let  us  define  the  Rayleigh  coefficient  K  as 
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Then 

T.  E. 

(  121)  Ke(8i#ea,<l>8)  =  2[Ki  cos  e#  cos  <J>8  -  Kz  cob  0,  COS  <bfl  1 

and 

K<j)(8t#  ®84>8)  =  2[-Kj  sin  <J>8  +  K2  cos  0^  cos  08  sin  <t>8 

+  1^3  s  in  0j  s  in  8_  J  ; 

T.  M. 

(  122)  K0(0i#  0a#  4>s)  =  2[-K2  cos  0*  cos  08  sin  <J>8  -  Kj  sin  0|  sin  0fi 

+  'ki  sin  <j>8  J 

and 

K<J>(0i»0s»4>8)  =  2[-Ka  cos  0i  cos  4>s  +  cos  08  cos  <|>8]  . 
For  backscatter  (0g  =  0^,  <£s  =  tt/2) 

T.  E. 

(123)  1^(9)  =  2[-Ki  +  k2  cos2 9  +  Kj  ai„26] 

and 

T.  M. 


Ke(0)  =  2[-K2  cos2t>  -  sin20+  £,]  . 


175 


It  is  clear  from  Eqs.  (  121-123)  that  a  knowledge  of  the  dipole  coef- 
ficients  Ki  ...»  Kj  . .  .  ,  suffices  to  determine  the  Rayleigh  coefficient 
for  any  polarization  and  for  any  arbitrary  orientation  of  source  and 
receiver. 

For  the  general  spheroid,  the  dipole  coefficients  may  all  be 
obtained  from  the  quantities  Ia  and  1^,  these  are  plotted  in  Figs.  39. 
40,  41,  and  42  as  a  function  of  the  axial  ratio  for  the  prolate  and 
oblate  spheroids,  respectively.  For  the  perfectly  conducting  case 
we  note  that 

( 124)  k3  -  -  Ik, 

and 


Thus  all  the  dipole  coefficients  may  be  obtained  from  Kj  and  K3  . 
These  quantities  are  plotted  in  Figs.  43,  44,  45,  and  46  as  a 
function  of  the  axial  ratio  for  the  prolate  and  oblate  spheroids, 
respectively.  From  the  curves  in  Figs.  39-46,  the  Rayleigh  coef- 
licient  for  the  prolate  or  oblate  spheroid  with  any  combination  of 
constitutive  parameters  for  any  polarization  and  orientation  can  be 


quickly  obtained. 
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As  an  example  of  the  use  of  these  curves,  consider  a  perfectly 
conducting  prolate  spheroid  with  a  2:1  axial  ratio.  From  Figs.  4  3 
and  45,  Kj  ~  &2  and  K3  ~  3.  8.  Then  from  Eq.  (  124)  ,  Kj  =  -1.  13 

and  K3  =  -0.  8.  Then  from  Eq.  (  123) 

(125)  Kte(8)  =  2[-l  .62  -  1.13  cos20  -  0.8  sin20] 

► 

and 

Ktm(*^  =  2[-  1 .62  cos2  8-  3.8  sin2 0-1. 13]  . 

j 

In  Table  V  the  Rayleigh  coefficients  are  given  for  perfectly  con¬ 
ducting  spheroids  with  axial  ratios  from  1(  1)10  and  incidence  angles 
Ofrom  0(  10°)  90°  for  both  the  T.  E.  and  T.  M.  polarizations. 


RAYLEIGH  COEFFICIENTS  OF 
PERFECTLY  CONDUCTING  PROLATE  SPHEROIDS 
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Fig.  ‘11- -The  quantity  K3  for  a  perfectly  conducting  prolate 
spheroid  as  a  function  of  axial  ratio. 
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PERFECTLY  CONDUCTING 
OBLATE  SPHEROID 


Fig.  46--The  quantity  K3  for  a  perfectly  conducting  oblate 
spheroid  as  a  function  of  axial  ratio. 


APPENDIX  II 

ON  THE  NATURE  OF  THE  CREEPING  WAVE 


In  the  derivation  of  transient  and  impulse  response  approxi¬ 
mations  for  smootl  curved  objects,  a  model  for  the  creeping  wave 
contribution  is  an  important  consideration.  High  frequency  asymp¬ 
totic  solutions  have  been  employed  with  some  success  in  the  pre¬ 
diction  of  the  diffraction  from  smooth  convex  curved  surfaces.  It 
is  of  interest  to  examine  certain  of  these  solutions  in  the  time  do¬ 
main,  where  additional  insight  both  to  the  scattering  mechanism 
and  to  the  nature  of  the  approximation  may  possibly  be  gained.  In 
this  Appendix,  attention  will  be  given  to  the  perfectly  conducting 
circular  cylinder  and  the  perfectly  conducting  sphere.  The  geo¬ 
metrical  theory  of  diffraction  solutions  given  by  Kouyoumjian[  43] 
and  Peters[  44]  for  the  cylinder  and  sphere  and  by  Keller[  45]  for 
the  sphere  '  will  be  examined  in  addition  to  Senior's[  46]  asymp¬ 
totic  expression  for  the  conducting  sphere  creeping  wave. 


*  Keller's  result  is  actually  for  axial  incidence  on  a  spheroid, 
but  considered  first  is  the  specialization  to  the  sphere. 
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A  general  observation  concerning  the  geometrical  diffraction 
theory  can  be  made.  For  conducting  objects,  each  surface  ray  mode 
according  to  the  theory  has  associated  with  it  a  phase  delay  of  at  least 
e  J  ,  where  k  is  the  free  space  wave  number  and  i  the  geodesic 
path  length  along  the  surface  of  the  diffracted  ray.  Thus  the  geo¬ 
metrical  theory  expression  for  the  diffracted  field,  when  transformed 
to  the  time  domain,  must  have  a  time  delay  of  at  least  l / c  provided 
the  complete  expression  is  one  which  leads  to  a  causal  time  function. 

As  will  be  seen,  there  is  an  essential  difference  between  the  geometri¬ 
cal  theory  expressions  for  the  cylinder  and  sphere.  The  geometrical 
expression  for  the  diffracted  field  of  a  cylinder  (  normal  incidence) 
leads  to  a  causal  time  function  whereas  the  corresponding  expression 
for  the  sphere  does  not.  This  difference  is  caused  by  the  presence  of 
a  cylindrical  caustic  for  the  sphere  case,  for  which  a  phase  jump  must 
be  added.  If  a  total  caustic  correction  of  other  than  -  nn  radians  is 
dictated  by  the  theory,  the  resultant  expression  for  the  diffracted  field 
apparently  leads  to  a  noncausal  time  function.  This  is  true  for  edge  - 
as  well  as  surface-diffracted  fields.  In  the  frequency  domain,  the  non¬ 
causal  nature  of  the  inverse  transformed  diffracted  field  is  of  no  conse¬ 
quence  and  in  fact  the  caustic  correction  is  essential  for  the  success 
of  the  theory.  In  the  time  domain,  however,  when  one  attempts  to 
construct  a  model  for  the  transient  or  impulse  response  waveform, 
the  presence  of  a  noncausal  time  function  is  not  convenient. 
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Consider  the  far-field  backscattcred  impulse  response  waveforms 
of  the  conducting  sphere  and  the  infinite  conducting  cylinder  for  broad¬ 
side  incidence.  The  sphere  waveform  is  shown  in  Fig.  4  7,  the  cylin¬ 
der  waveforms  for  parallel  polarization  (  solid)  and  for  perpendicular 
polarization  (  dashed)  are  shown  in  Fig.  48.  These  approximate  wave¬ 
forms  were  obtained  from  the  Fourier  synthesis  procedure[  13]  using 
the  exact  solutions  for  the  respective  scattering  problems.  Fo:.*  the 
sphere,  475  harmonics  with  a  fundamental  of  0.  2512  sphere  circumfer¬ 
ences  in  wavelengths  were  used  For  the  cylinder,  79  harmonics  with 


Fig.  4  7  -  -  A  pp  r  ox  ima  te  impulse  response  waveform  of 
conducting  sphere,  backscatter. 
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Fig.  48- -Approximate  impulse  response  waveform  of  conducting 
infinite  circular  cylinder,  backscatter,  broadside  incidence. 
Parallel  (  solid)  and  perpendicular  (  dashed)  polarization. 


a  fundamental  of  0.  157  cylinder  circumferences  in  wavelengths  were 
used.  In  both  cases,  the  pulse  width  and  pulse  separation  of  the 
synthesized  incident  cosine  pulse[  13]  are  sufficient  to  ensure  a 
reasonably  accuracte  approximate  to  the  impulse  response  waveform. 
Assuming  the  progression  of  creeping  wave  contributions  along 


#  An  adequate  representation  of  the  sphere  impulse  waveform  can 
be  obtained  from  as  few  as  64  harmonics. 
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geodesic  paths,  the  creeping  wave  contributions  for  both  the  sphere 
and  cylinder  should  appear  first  at  a  time  t1  =  (1  +  ir/2)tQ.  For  the 
sphere  (  Fig.  47)  ,  the  creeping  wave  contribution  is  a  secondary  peak 
occurring  at  the  prescribed  time.  For  parallel  polarization  on  the 
cylinder  (  Fig.  48)  ,  where  the  incident  electric  vector  is  effectively 
shorted  out  at  the  shadow  boundary,  there  is  no  evidence  of  a  creeping 
wave  contribution.  For  perpendicular  polarization  the  contribution 
begins  at  the  above  prescribed  time.  The  distinct  difference  in  the 
creeping  wave  contributions  for  the  sphere  and  cylinder  is  vividly 
evident  in  Figs.  47  and  48. 

One  should  note  particularly  that  the  creeping  wave  contribution 
for  the  sphere  peaks  at  t'  “  (  1  +  ir/2)tQ,  whereas  the  contribution  for 
the  cylinder  (perpendicular  polarization)  begins  at  this  same  time. 
With  reference  to  either  the  sphere  or  cylinder  waveforms,  the  wave¬ 
form  for  t'  >  (  1  +  ir/2)  tQ  might  properly  be  explained  by  a  creeping 
wave  theory,  since  the  creeping  wave  cannot  arrive  before  this  time. 
Thus  one  would  conclude  that  a  creeping  wave  theory  would  be  inade¬ 
quate  for  the  sphere  since  it  could  not  explain  the  leading  edge  of  the 
peak,  but  adequate  for  the  cylinder.  It  will  be  shown,  however,  that 
the  referenced  sphere  solutions  predict  "precursors"  in  the  time  do¬ 
main  which  agree  approximately,  in  the  neighborhood  of  the  peak,  with 


the  waveform  in  Fig.  47. 


191 


Consider  first  the  geometrical  theory  of  diffraction  solution  for 
the  conducting  cylinder  of  radius  a  for  broadside  incidence.  Both 
Kouyoumjianf  43]  and  Peters[  44]  give  the  result  (  retaining  only  one 
term  in  the  diffraction  and  attenuation  coefficients  and  neglecting  all 
but  the  first  surface  ray  mode) 

.  ir 

(126)  F.c.w  =  0.  81  a1/3  \l/6  e“J  12 

( J<*> ) 

.  IT 

-  [  jkira  +  1.  65e-i £  a“2^ \-1^3  -rra] 

•  e  . 

sTr 


Using  an  s  notation  (  s  =  jut0  =  j2ka)  ,  shifting  the  phase  reference  to 
the  specular  point  of  the  cylinder,  and  normalizing  such  that  the  echo 
width  predicted  for  the  creeping  wave  contribution  is 


(  127) 

cr(  oj)  =  \Trr  |  GCW(  jw  )  | 

(  128) 

„  a  e~[  06  +  P  s 

~cw,  .  ,  A  e 

G  (ju,)  --  s„6 

where 

A  =0.  8lN/Ta  (  4ti)  1,6 

a  =  (  1  +  it/2)  tQ, 
_  1.  65tt^3 


and 
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Friedlander[  47]  has  obtained  an  approximate  general  result  for  the 
inverse  Laplace  transform  of  functions  of  the  type  in  Eq.  (  128)  using 


a  saddle  point  technique. 


the  result  in  Eq.  (  129)  holding  for  T  small  and  n  2.  0.  Applying  Eq. 

(  129)  to  Eq.  (  128) ,  the  cylinder  creeping  wave  impulse  response 
waveform  predicted  by  the  geometrical  theory  of  diffraction  is  approxi¬ 


mately 


(130) 


F 


t'  -  (  l  +  ir/2)  tQ 


•  u(  t'  -  (  1  +  ir/2)  tQ) 


where  Eq.  (  130)  holds  for  t'/t0  -1  -  ir/2  small.  Note  particularly 
that  the  geometrical  theory  has  led  to  a  causal  time  function,  i.  e.  , 
it  is  identically  zero  for  t'  <  (  1  +  ir/2)tQ.  Equation  (  130)  yields  a 
good  estimate  of  the  creeping  wave  contribution  for  t'  -  ( 1  +  ir/2)t0 


small. 
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The  impulse  response  waveforms  obtained  by  Fourier  synthesis 
for  the  conducting  circular  cylinder  shown  in  Fig.  48  are  felt  to  be  of 
some  interest  in  their  own  right*  without  regard  to  the  comparison  to 
the  geometrical  theory.  To  the  best  of  the  author's  knowledge*  such 
waveforms  have  not  appeared  previously  in  the  literature.  For  paral¬ 
lel  polarization  it  is  evident  that  except  at  the  lowest  frequencies*  the 
geometrical  optics  specular  contribution  is  an  excellent  approximation. 
For  perpendicular  polarization  Peters[44]  has  suggested  the  combi¬ 
nation  of  the  creeping  wave  component  in  Eq.  (  126)  and  the  specular 
term  as  an  approximation;  and  demonstrated  the  validity  of  this  model 
for  cylinder  radii  greater  than  0.  3  wavelengths.  Again,  it  is  evident 
from  Fig.  48  that  such  a  model  must  fail  at  frequencies  somewhat 
higher  than  for  the  specular  model  for  parallel  polarization  and  that 
the  failure  of  the  model  is  not  due  to  inadequacies  of  the  creeping  wave 
component*  but  rather  to  the  omission  of  additional  terms  for  the 
optical  contribution.  Finally  it  is  noted  that  for  perpendicular  polari¬ 
zation  in  the  neighborhood  of  t'  =  0+,  the  details  of  the  response  wave-, 
form  are  somewhat  obscured  because  of  the  finite  width  of  the  incident 
quasi-impulsive  wave  and  the  occurrence  of  both  the  specular  jump 
and  a  second  opposite  sign  discontinuity  at  the  same  instant  Details 
of  the  response  waveform  in  this  region  are  not  germain  to  the  present 
study,  and  will  not  be  pursued  except  to  note  that  the  synthesis 
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procedure  could  be  used  to  resolve  the  question  by  the  use  of  additional 
harmonics  to  narrow  the  incident  pulse  and  by  subtracting  out  the  well 
known  specular  contribution. 

Turning  now  to  the  conducting  sphere  of  radius  a,  it  can  be  shown 
that  the  various  high  frequency  asymptotic  estimates  for  the  creeping 
wave  reduce  to  the  following,  where  in  each  case  the  phase  reference 
has  been  shifted  to  the  specular  point  of  the  sphere  and  the  expression 
normalized  such  that  the  echo  area  predicted  for  the  creeping  wave  is 
given  by 

(131)  T<“>  =  TT  I  GCW(  B)  (  * 


From  Peters[  44] 

eXAe-(°'  +  ?'U’l 
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From  Keller[  45]  ,  specializing  the  spheroid  solution  to  the 


sphere  case 


(  133) 


;CW,  v  _  ;T  Ifl  ’  [  «  8  +P  SWS] 
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(  \Z)xt>  q,D2 


p  =  -3^—  ,  q,  -  1.  469  3  54,  and 
(  12)  "* 


Dz  =  (  1.  1 6680) * 


It  can  be  shown[  48]  that  when  the  same  number  of  terms  are 
retained,  Kouyoumjian1  s[ 43]  geometrical  theory  solution  for  the  creep¬ 
ing  wave  and  Senior's[  46]  asymptotic  high  frequency  expansion  for 
the  creeping  wave  are  identical.  Using  only  the  first  term,  and  the 
acoustically  hard  component  which  dominates,  there  is  obtained 


(  134) 
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196 


q0  =  1.  01879 


Noting  the  distinct  differences  in  the  various  formulas  given 
above,  it  should  be  pointed  out  that  each  of  the  results  has  been  shown 
to  yield  some  degree  of  agreement,  when  combined  with  the  optical 
components,  with  the  well  known  exact  solution  for  the  sphere  back** 
scattering  cross  section.  Keller's  expression  lacks  the  correction 
terms  to  the  diffraction  and  attenuation  coefficients  which  were 
included  by  Kouyoumjian.  It  was  demonstrated  indirectly  by  the 
author[  49]  ,  in  a  study  of  the  cone- sphere  axial  cross  section,  that 
Keller:s  result  for  the  sphere  creeping  wave  was  incomplete.  The 
purpose  here  however  is  not  a  comparison  or  criticism  of  the  various 
results  in  the  frequency  domain  but  an  examination  of  the  general  type 
of  waveform  predicted  in  the  time  domain. 

Note  in  the  expressions  in  Eqs.  (  132) ,  (  133) ,  and  (  134) ,  that 
they  are  of  the  general  form 

GCW(  s)  =  t  jF(  s) 


(  135) 


f 
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where  F(  s)  differs  for  the  three  formulas.  In  Peter's  case  this  is 

v 

one  of  the  two  terms  obtained  when  e^  is  expanded.  Assuming  that 
F(  s)  yields  a  causal  time  function,  which  will  be  shown  later,  it  is 
seen  that  the  asymptotic  estimates  of  the  creeping  wave  add  a  con¬ 
stant  tr/2  radians  phase  shift.  If  F(  s)  yields  a  causal  time  function 

then  Gcw(  s)  is  noncausal.  To  show  this, 

A+j® 

(  136)  Fy  w(  t * )  =  —  C  Gcw(  s)  est'  ds 

1  2irj  J 

A-j« 


GCW(jw )  e^wt  du> 


f 


where  GCW(-s)  =  Gcw(  s)*  has  been  utilized.  To  obtain  a  real  time- 

dependent  waveform,  the  constant  phase  shift  is  assumed  to  change 

sign  for  negative  frequencies.  From  Papoulis[  50]  , 

00 

(137)  ]ejo,t'  d«=  it  6(  t')  +  j-i- 

0 


thus 


(  138) 


If  the  inverse  transform  of  F{  s)  is  fQ(  t * )  ,  i.  e.  , 
(i39)  f0<  *’) 


then  by  convolution 
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and  Eq.  (  143)  reduces  to 
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(  145)  jf 


mi 


A  e“t°  8  +  P8^] 

sfl  1^3 


r  2p>fe  | 

AJ- 


z./57s'“  it'..)'" 


u(t'-a)  , 


which  is  exactly  the  result  obtained  by  using  Friedlander's  approxi¬ 
mate  formula  (  Eq.  (  129) )  in  Eq.  (  141) .  From  Eqs.  (  142)  and  (  143) , 


the  exact  expression  for  Peters'  creeping  wave  in  the  time  domain  is 
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and  for  t'-a  small  this  reduces  to 
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For  Keller's  expression  in  Eq.  (133) 


(148) 
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From  Eq.  (  142)  and  the  differentiation  formula 
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where  the  prime  indicates  differentiation  with  respect  to  t'  .  The  exact 
expression  in  the  time  domain  for  the  creeping  wave  predicted  by  Keller 


is 


(  150) 
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and  for  t'-o  small, 
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Unfortunately,  the  inverse  transform  of  these  functions  can  only  be 
written  in  terms  of  an  infinite  series.  Using  Eq.  (  140)  ,  the  creeping 
wave  time  domain  response  of  the  Senior  result  ( Eq.  (  144))  can  be 
written,  but  the  explicit  results  are  of  little  interest  because  of  their 
complexity. 

The  noncausal  nature  of  the  impulse  response  waveforms  pre¬ 
dicted  for  the  creeping  wave  on  the  sphere  is  demonstrated  in  Eqs. 

(  146)  and  (  150)  ,  and  implicitly  from  a  combination  of  Eqs.  (  152) 
and  (  153).  This  is  also  evident  from  Eq.  (  135)  alone,  once  the 
causal  nature  of  the  inverse  transform  of  the  frequency-dependent 
function  without  the  constant  phase  shift  is  demonstrated.  It  is  further 
interesting  to  note  that  in  the  applications  of  the  geometrical  theory, 
i.  e.  ,  Keller,  and  Kouyoumjian,  this  constant  phase  shift  is  included 
because  of  a  caustic  correction,  and  is  not  a  part  of  the  diffraction  coef 
ficient  obtained  from  the  solution  of  a  canonical  problem.  However, 
Senior's  expression  for  the  creeping  wave  also  includes  the  constant 
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phase  shift  and  in  this  case  results  from  simply  an  asymptotic  high  fre¬ 
quency  expansion  of  the  exact  solution.  In  time  domain  terms,  the  non- 
causal  nature  is  that  of  a  network  required  to  yield  a  tt/ 2  radians  phase 
shift  at  all  frequencies,  ^he  impulse  response  of  such  a  network  would 
have  a  "precursor",  and  the  network  itself  is  nonphysical. 

Before  further  comment,  the  time  domain  response  waveforms 
predicted  by  certain  of  the  asymptotic  approximations  will  be  calculated. 
This  is  done  most  simply  for  the  results  of  Senior  and  Keller,  using  the 
Fourier  synthesis  procedure,  since  calculations  of  the  frequency- 

jjc 

dependent  expressions  at  *he  requisite  harmonics  were  available.  In 
both  cases,  two  results  are  obtained;  the  causal  waveform  corresponding 
to  removal  of  the  tt/2  radians  phase  shift,  and  the  noncausal  waveform 
corresponding  to  the  complete  expression.  The  waveforms  for  the 
causal  case  are  shown  in  Fig.  49.  Included  in  Fig.  49  are  calculated 
results  from  the  analytical  expression  in  Eq.  (  149)  with  the  simplifi¬ 
cation  of  the  modified  Hankel  function  for  large  arguments  (  Eq.  (  144) ) 
employed.  In  Figs.  50  and  51,  the  corresponding  noncausal  results 
for  the  actual  asymptotic  high  frequency  expansions  of  Keller  and 


#  Professor  L.  Feters,  Jr.  of  The  Ohio  State  University  supplied 
calculated  data  from  Senior's  results.  Professor  E.  M.  Kennaugh 
had  previously  obtained  calculated  data  from  Keller's  result  and 
had  also  obtained  a  portion  of  the  time  domain  results  via  Fourier 
synthesis. 
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Fig.  49--Time  dependent  response  waveforms  predicted 
by  Keller-Levy  (  solid)  and  Senior  (  dashed)  asymptotic 
estimates  of  creeping  wave  on  conducting  sphere  when 
g-jtr/2  phase  factor  removed.  Backscatter. 
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Fig.  50--Time-dependent  response  waveform  predicted 
by  Keller- Levy  asymptotic  estimate  of  creeping  wave 
on  conducting  sphere,  backscatter. 

Senior  respectively  are  shown.  Aside  from  the  demonstrated  agree¬ 
ment  of  the  analytical  result  for  short  time  (t'-(  1  +  ir/2)t0)  small  in 
Fig.  49)  ,  the  significant  result  in  these  figures  is  that  the  constant 

phase  shift  of  ir/2  radians  (  which  amounts  to  convolution  with  l/t' 

* 

in  the  time  domain)  has  transformed  the  causal  waveforms  into  non- 
causal  ones  in  approximate  agreement  in  character  with  the  creeping 
wave  peak  shown  in  Fig.  47.  It  is  apparent  that  the  addition  of  the 
optical  contributions  to  the  noncausal  response  waveform  in  Fig.  51 
will  yield  approximate  agreement  in  form  with  the  sphere  impulse  re¬ 
sponse  waveform.  This  comparison  is  made  in  Fig.  52.  Shown  are 
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Fig.  52- -Impulse  response  waveform  of  conducting  sphere 
backscatter.  Solid,  Fourier  synthesis  from  Mie  Series. 
Dashed,  Fourier  synthesis  from  Senior's  creeping 
wave  estimate  plus  optical  terms. 


207 


the  Fourier  synthesis  derived  waveforms  obtained  from  the  exact 
solution,  and  from  Senior's  expression  plus  the  optical  terms.  It 
is  easily  seen  from  Fig.  52  that  the  combination  of  Senior's  creep¬ 
ing  wave  expression  with  the  optical  contributions  yields  a  reason¬ 
able  approximation  in  the  neighborhood  of  the  creeping  wave  peak 
to  the  sphere  impulse  response  waveform. 

In  one  sense,  the  results  in  this  Appendix  undoubtedly  are  a 
reiteration  of  what  is  already  known,  namely  that  asymptotic  esti¬ 
mates  of  the  creeping  wave,  when  combined  with  the  optical  contri¬ 
butions,  yield  a  fair  estimate  of  the  echo  area  of  spheres  and 
cylinders  over  a  range  of  higher  frequencies.  From  the  time  do¬ 
main  viewpoint,  however,  the  results  explain  a  paradox  which  had 
existed  as  to  how  the  asymptotic  estimates  of  the  creeping  wave  with 
the  addition  of  optical  terms  could  possibly  yield  good  cross  section 
data  in  the  resonance  region  for  the  sphere.  The  explanation  lies  in 
the  noncausal  nature  of  the  predicted  creeping  wave  waveform.  That 
is,  it  was  clear  that  if  a  response  waveform  beginning  at 
f  =(  1  +  ir/2)t  were  added  to  the  optical  contributions  (  a  negative 
delta  function  and  a  positive  step  function) ,  the  result  would  not  be 
a  good  approximation  for  the  impulse  response  waveform  of  the 
sphere.  But  the  addition  of  a  noncausal  waveform  can,  as  has  been 
shown,  yield  a  reasonable  approximation. 
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The  time  dependent  response  waveforms  predicted  by  high  fre¬ 
quency  asymptotic  estimates  of  the  creeping  waves  on  spheres  and 
circular  cylinders  have  been  demonstrated.  For  the  circular  cylin¬ 
der  a  closed  form  expression  for  the  waveform,  holding  for  short 
times  after  the  arrival  of  the  diffracted  front,  was  obtained.  The 
complexity  of  the  asymptotic  expressions  precludes  such  a  result  for 
the  sphere.  In  this  case  a  unique  tool,  the  Fourier  synthesis  tech¬ 
nique,  was  used  to  demonstrate  the  character  of  the  response.  This 
same  approach  can  be  used  to  resolve  the  question  of  the  exact  nature 
of  the  sphere  response  in  the  neighborhood  of  the  secondary  peak. 

That  is,  with  the  Fourier  synthesis  approach  it  is  necessary  to  use 
many  harmonics  to  secure  a  sufficiently  narrow  pulse  width  to  ap¬ 
proximately  resolve  a  jump  discontinuity  or  a  cusp  type  behavior 
of  the  waveform.  Using  the  Mie  series,  it  is  not  economically  feasi¬ 
ble  to  secure  the  requisite  harmonics.  However,  the  asymptotic 
estimates  are  easily  calculated  and  it  is  a  relatively  simple  matter 
to  secure  the  necessary  data.  Some  preliminary  work  on  this 
question  has  been  done,  utilizing  up  to  1800  harmonic*  with  a 
fundamental  corresponding  to  a  sphere  circumference  of  0.  1  wave¬ 
lengths.  These  results  represent  an  order  of  magnitude  decrease 
in  the  narrowest  pulse  width  obtained  using  Mie  series  calculations. 
These  results  are  incomplete  and  will  not  be  detailed.  However,  it 
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was  concluded  that  the  secondary  peak  is  finite  (  a  previously  unsolved 
question)  and  that  an  additional  order  of  magnitude  reduction  of  the 
pulse  width  would  permit  the  magnitude  of  the  peak  to  be  determined 
as  well  as  the  slope  of  the  waveform  on  either  side  of  the  peak.  It  is 
clear  therefore  that  asymptotic  high  frequency  estimates  of  diffracted 
waves  can  be  utilized  in  the  time  domain  to  resolve  questions  concern¬ 
ing  the  character  of  the  response  waveform  at  times  corresponding  to 
the  arrival  of  the  diffracted  front.  The  Fourier  synthesis  procedure 
is  an  indispensible  tool  in  this  case  since  many  of  the  asymptotic  esti¬ 
mates  are  not  amenable  to  an  analytical  treatment.  However,  as  a 
device  for  "understanding"  the  scattering  mechanism  of  an  object, 
those  asymptotic  estimates  which  employ  a  nonphysical  model  and 
predict  a  noncausal  waveform  leave  a  great  deal  to  be  desired.  For 
the  specific  case  of  the  sphere,  both  time  and  frequency  domain 
results  demonstrate  that  the  estimate  of  the  energy  arriving  via 
creeping  wave  paths  is  correct.  But  the  interpretation  that  these 
paths  are  confined  to  great  circles  is  incorrect.  Unless  one  is  willing 
to  accept  group  velocities  for  the  diffracted  fields  which  are  greater 
than  the  speed  of  light,  the  near  time  build-up  of  the  creeping  wave 
peak  can  only  be  explained  on  the  basis  of  nongreat  circle  or  non¬ 
geodesic  paths.  Note  that  in  terms  of  geometry,  the  difference 
between  the  cylinder  and  sphere  is  that  there  are  no  possible  shorter 
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paths  on  the  cylinder  (  i.  e.  ,  shorter  than  the  cylinder  circumference) 
whereas  such  shorter  paths  do  exist  on  the  sphere.  It  is  concluded 
therefore  that  the  physical  mechanism  for  creeping  waves  actually 
consists  of  propagation  along  nongc^desic  as  well  as  geodesic  paths 
with  the  attenuation  along  the  paths  increasing  with  departure  from 
the  geodesic.  It  is  suggested  that  an  examination  of  the  time-de- 
pendent  currents  on  the  shadow  side  of  the  sphere  when  it  is  illumi¬ 
nated  by  an  incident  impulsive  wave  will  reveal  the  location  and 
relative  attenuation  of  these  paths.  Some  evidence  to  support  this 
has  been  found  using  calculations  of  the  transient  currents  obtained 
by  Kennaugh[  53]  .  It  is  estimated  that  the  additional  calculations  of 
Mie  series  and  corresponding  Fourier  synthesis  procedure  required 
to  produce  conclusive  evidence  of  the  nongeodesic  paths  as  well  as 
estimates  of  the  attenuation  would  cost  approximately  $300.  00  on  an 
IBM  7094  computer. 

An  alternative  approach  confined  entirely  to  the  frequency  do¬ 
main  may  also  be  possible.  The  asymptotic  estimate  of  the  creeping 
wave  derived  by  Senior  is  not  unique.  This  suggests  that  a  different 
asymptotic  expansion  corresponding  to  a  causal  response  waveform 
may  be  possible.  Such  an  expansion,  while  of  no  greater  utility  than 
Senior's  result  in  the  frequency  domain,  would  be  of  much  greater 
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use  in  the  time  domain.  From  a  philosophical  viewpoint,  the  high- 
frequency  asymptotic  estimates  of  diffracted  fields  must  be  consider¬ 
ed  incomplete,  insofar  as  the  time  domain  picture  is  concerned, 
until  such  estimates  can  be  interpreted  in  the  time  domain  with  a 
strictly  causal  relationship  between  the  geometry  of  the  object  and 
specific  locations  in  time  of  the  response  waveform.  It  is  hoped 
that  the  results  in  this  Appendix  may  provide  the  stimulus  for  future 
investigations  in  this  direction. 
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